1/f noise: implications for solid-state quantum information 
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The efficiency of the future devices for quantum information processing wih be hmited 
mostly by the finite decoherence rates of the individual qubits and quantum gates. Re- 
cently, substantial progress was achieved in enhancing the time within which a solid-state 
qubit demonstrates coherent dynamics. This progress is based mostly on a successful 
isolation of the qubits from external decoherence sources obtained by clever engineering. 
Under these conditions, the material-inherent sources of noise start to play a crucial 
role. In most cases, quantum devices are affected by noise decreasing with frequency 
/ approximately as 1/f. According to the present point of view, such noise is due to 
material- and device-specific microscopic degrees of freedom interacting with quantum 
variables of the nanodevice. The simplest picture is that the environment that destroys 
the phase coherence of the device can be thought of as a system of two-state fluctuators, 
which experience random hops between their states. If the hopping times are distributed 
in a exponentially broad domain, the resulting fluctuations have a spectrum close to 1/f 
in a large frequency range. 

In this paper we review the current state of the theory of decoherence due to degrees 
of freedom producing 1/f noise. We discuss basic mechanisms of such noises in vari- 
ous nanodevices and then review several models describing the interaction of the noise 
sources with quantum devices. The main focus of the review is to analyze how the 1/f 
noise destroys their coherent operation. We start from individual qubits concentrating 
mostly on the devices based on superconductor circuits, and then discuss some special 
issues related to more complicated architectures. Finally, we consider several strategies 
for minimizing the noise-induced decoherence. 
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I. INTRODUCTION 

Evidence of properties that fluctuate with spectral den- 
sities varying approximately as 1/f over a large range of 
frequencies, /, has been reported in an astonishing va- 
riety of systems. In condensed matter physics, the dif- 
ficulties in reasonably explaining the shape of the spec- 
trum and in ascribing a physical origin to the noise in 
the diversity of system where it has been observed, have 
kept 1/f noise in the forefront of unsolved problems for 
a long time. The large theoretical and experimental ef- 
fort in this direction up to the late '80s, with emphasis 
on 1/f conductance fluctuations in conducting materials, 
has been reported in the excellent reviews (Dutta and 



Horn 1981 



Weissman, 1988), book (Kogan 1996) and 



others cited therein. 

With the progressive reduction of systems size, fluctu- 
ations having l//-like spectra have been frequently ob- 
served in various mesoscopic systems. The importance of 
magnetic flux noise in Superconducting Quantum Inter- 
ference Devices (SQUIDs) was recognized already in the 
1980s (iKoch et al\ 119831 iWeissmanl 119881) thus opening 



the debate about its physical origin - noise from the sub- 
strate/mount or noise from trapped flux in the SQUID 



frequency polarization fluctuations, see, e. g., (Krupenin 



et 



1998, 2001 Wolf et a/., 1997 Zorin et a/., 1996). 



Nanodevices are the subject of intense research at 
present also because of their long-term potential for 
quantum information. Similarly to atomic systems, the 
quantum nature of nanocircuits, despite being hundreds 
of nanometers wide and containing at huge number of 
electrons, is observable. Because of these characteris- 
tics solid state quantum bits (qubits) can be relatively 
easily addressed to perform desired quantum operations. 
The drawback of tunability is sensitivity to fluctuations 
of control parameters. Fluctuations are partly extrinsic, 
like those due to the local electromagnetic environment. 
This source of noise has been greatly reduced via clever 
engineering. In almost all quantum computing nanode- 
vices fluctuations with 1/f spectral density of different 
variables and of different physical origin have been ob- 
served. There is clear evidence that 1/f noise is detri- 
mental to the required maintenance of quantum coherent 
dynamics and represents the main source of decoherence. 
This fact has stimulated a large effort of both the exper- 
imental and theoretical communities aimed on one side 
at a characterization of the noise, on the other at un- 
derstanding and eventually reducing noise effects. On 
a complementary perspective, nanodevices are sensitive 
probes of the noise characteristics and therefore may pro- 
vide important insights into its microscopic origin. 

In this review we describe the current state of theo- 
retical work on 1/f noise in nanodevices with emphasis 
on implications for solid state quantum information. We 
will focus on superconducting systems and refer to other 
implementations, in particular, those based on semicon- 
ductors, whenever physical analogies and/or formal sim- 
ilarities are envisaged. 

Previous reviews on 1/f noise mainly focused on re- 
sistivity fluctuations of conducting materials ( Dutta and 



Horn 1981 



Kogan[ 


1996 



Weissman 



1988). Important 



questions have been addressed, like universality of the 
mechanisms leading to conductivity fluctuations and of 
the kinetic patterns leading to the 1/f spectral form. 
Detailed investigations in different materials (metals and 
semiconductors) failed to confirm the appealing impres- 
sion of universality and lead instead to the conclusion of 
the existence of a variety of origin of 1/f conductance 
noise in diverse materials. 

Recent experiments with superconducting circuits ev- 
idenced 1/f low- frequency fluctuations of physically dif- 
ferent observables, thus providing new important in- 
sights into noise microscopic sources. This is due to 
the fact that three fundamental types of superconducting 
qubits exist: flux, charge and phase, for a recent review 



and temperature dependence ( |Savo et al\ |1987[ |Well-| see ([Clarke and Wilhelm[ |2008[ |Ladd et al\ |2010[ [Steffen 



stood et at. 1987b). Single-electron and other tunneling et al. 2011 You and Nori, 2011). The main difference 



devices have provided a compelling evidence that fluc- 
tuating background charges, either within the junctions 



between them is the physical observable where informa- 
tion is encoded: superconducting current, excess charge 



in a superconducting island or the superconducting phase 
difference across a Josephson junction. Different observ- 
ables couple more strongly to environmental variables of 
different nature and therefore are sensitive probes of dif- 
ferent noise sources. As a result, magnetic flux noise, 
polarization or "charge" noise and critical current noise 
with 1// spectrum in some frequency range are presently 
routinely measured in the three implementations. 

One scope of this review is to present the current state 
of understanding of the microscopic sources of 1// noise 
in superconducting nanocircuits. Despite relevant mech- 
anisms have been largely identified, in most solid state 
nanodevices this problem cannot be considered as totally 
settled (Section In]). In some cases, available experiments 
do not allow drawing solid conclusions and further in- 
vestigation is needed. A number of basic features of the 
phenomenon are however agreed upon. According to pre- 
vious reviews, 1// noise results from a superposition of 
fiuctuators (whose nature has to be specified case by case) 
having switching times distributed in a very broad do- 
main (Dutta and Horn, 1981[ Kogan[ 1996 Weissman 



1988). Statistical properties of 1// noise have been dis- 
cussed in depth in the book by iKogan (1996). These 



properties are at the origin of 1// noise-induced loss of 
coherence of solid-state qubits. For the sake of clarity, 
here we recall the basic definitions and specify their use 
in the context of the present review. 

For a Gaussian random process all nonzero n-th order 
moments can be expressed in terms of the second-order 
moments, i. e., pair correlations. In general, statistical 
processes producing 1// noise are non-Gaussian. This 
fact has several important implications. On one side, 
statistical correlations higher than the power spectrum 
should be considered in order to characterize the process, 
see (Kogan 1996), paragraph 8.2.2. On the other side. 



deviations from Gaussian behavior are also expected to 
show up in the coherent quantum dynamics of solid state 
qubits. For Gaussian noise, random noise-induced phases 
acquired by a qubit obey the Gaussian distribution. We 
will refer to a process as Gaussian whenever this situation 
occurs. As we will see, even in cases where the noise can 
be considered as a sum of many statistically-independent 
contributions, the distribution of the phases can be essen- 
tially non-Gaussian. A number of investigations aimed 
at predicting decoherence due to non- Gaussian 1// noise: 



([Bergli et 



2006, 2009 



Burkard, 2009 Galperin et al. 



2006[ |2007| IGrishin et 0/112005, .Paladino et al\ |2002 



Yurkevich et al\ 2010) (more references can be found in 
Section 11111). 



The other crucial property of 1// noise is that it cannot 
be considered a Markovian random process. A statistical 
process is Markovian if one can make predictions for the 
future of the process based solely on its present state, just 
as well as one could do knowing the process's full history. 
We will see that even if the noise can be considered as a 
sum of Markovian contributions, the overall phase fluctu- 



ations of a qubit can be essentially non-Markovian. The 
consequence of 1// noise being non-Markovian is that the 
effects of 1// noise on the system evolution depend on the 
specific "quantum operation" and/or measurement pro- 
tocol. In this review we will illustrate various approaches 
developed in recent years to deal with the non-Markovian 
nature of 1// noise starting both from microscopic quan- 
tum models and from semi-classical theories. We will 
discuss the applicability range of the Gaussian approxi- 
mation as well as deviations from the Gaussian behavior 
in connection with the problem of qubit dephasing. 

A statistical process is stationary if all joint probability 
distributions remain invariable under identical shifts of 
all time points. To the best of our knowledge, at present 
time there is no clear evidence of non-stationarity of the 
processes leading to 1// noise (see the discussion in Sec- 
tion [ml). 



Low-frequency noise is particularly harmful since it is 
difficult to filter it out by finite band filters. In recent 
years different techniques have been proposed, and some- 
times experimentally tested, in order to limit the effect 
of low frequency fluctuations. One successful strategy 
to increase phase-coherence times is to operate qubits 
at working points where low-frequency noise effects van- 
ish to the lowest-order; such operating conditions are 
called "optimal point" or "magic point" (|Vion et al, 
2002). Further substantial improvement resulted from 



the use techniques inspired to Nuclear Magnetic Reso- 



nance (NMR) ( |Schlichter , 1992). In a quantum informa- 
tion perspective, any approach aimed at limiting deco- 
herence should be naturally integrated with other func- 
tionality, such as quantum gates. In addition, achieved 
fidelities should be sufficiently high to allow for the suc- 
cessful application of Quantum Error Correction codes. 
The question about the best strategy to limit 1// noise 
effects via passive or active stabilization is still open. We 
will review the current status of the ongoing research 



along this direction in Section III.D 



Building scalable multi-qubit systems is presently 
the main challenge towards the implementation of a 



solid-state quantum information processor (Nielsen and 



Chuang 1996) 



The effect of 1// noise in solid state 
complex architectures is a subject of current investiga- 
tion. Considerable improvement in minimizing sensitiv- 
ity to charge noise has been reached via clever engineer- 
ing. Within the newly framework named circuit quan- 
tum electrodynamics (cQED), based on the analogy of 
superconducting circuits and the fields of atomic physics 
and quantum optics, important steps further have been 
done. Among the newest we mention the achievement 
of three-qubit entanglement with superconducting nano- 



circuits (Di Carlo et a/., 2010 Neeley et al , 2010) which. 



in combination with longer qubit coherence, illustrate a 



potentially viable approach to factoring numbers ( Lucero 



et al. 



et al. 



2012 



2012[) , implementing quantum algorithms (Fedorov 
Mariantoni et al\ [2Q11|) and simple quantum 



error correction codes (Chow et a/., 2012 Reed et al. 



2012 Rigetti et a/., 2012). 



A. General features and open issues 

Low frequency noise is commonly attributed to so- 
called fluctuators - dynamic defects, which randomly 
switch between two metastable states (1 and 2), see, 
e. g., (Dutta and Horn 1981 Kogan[ p.996[ Weissman] 
19881). Such a switching produces Random Telegraph 



(RT) noise. The process is characterized by the switch- 
ing rates 71^2 and 72^1 for the transitions 1^2 and 
2^1. Only the fluctuators with energy splitting E less 
than a few ksT (T is temperature) contribute to the de- 
phasing of a qubit, since the fluctuators with large level 
splitting are frozen in their ground states. As long as 
E < ksT the rates 71^2 and 72^1 are close in magni- 
tude, and to describe general features of decoherence one 
can assume that 71^2 = 72^1 = 7, i- e., the fluctuations 
can be described as a RT process, see (Buckingham, 
1989| [Kirton and Urenj |1989| [Koganl |1996D for reviews. 



A set of random telegraph fluctuators with exponentially 
broad distribution of relaxation rates, 7, produces noise 
with a 1// power spectrum at 7min <C cj = 27r/ <C 7max- 
Here, 7min is the switching rate of the "slowest" fluc- 
tuator, whereas 7max is the maximal switching rate for 
fluctuators with energy difference E = ksT. Random 
telegraph noise has been observed in numerous nano- 
devices based on semiconductors, normal metals and su- 



perconductors (Duty et al.\ 2004 



Eroms et ( 



Rogers and Buhrmanl 



1984 



2006 



Par- 



man et al\ |1991| [Peters et al\ fl999| [Ralls et al\ fl984| 



1985). Multi-state fluctu- 



ators with number of states greater than two were also 
observed (Bloom et a/., 1993a|b ). 

Various microscopic sources can produce classical ran- 
dom telegraph noise. Here we briefly summarize some of 
them in connection with charge, flux and critical current 
noise and mention some recent key references; a detailed 
discussion will be presented in Section [TT[ 

The obvious source of RT charge-noise is a charge 
which jumps between two different locations in space. 
Various hypotheses about the actual location of these 
charges and the nature of the two states are still under 
investigation. The first attempt at constructing such a 
model in relation to qubit decoherence appeared in ( [Pal-[ 
adino et al. , 2002 ) , where electron tunneling between a lo- 
calized state in the insulator and a metallic gate was stud- 
ied. The quantitative importance (in explaining observed 
spectra) of quantum effects of hybridization between elec- 
tronic states localized at the trap and electrodes extended 



states was pointed out by Abel and Marquardt (2008) 



and Grishin et al (2005). Models considering the ac- 



tual superconducting state of the electrodes have been 
studied in (Faoro et al. 2005 Faoro and loffe 2006), 



tal observations of charge noise based on measurements 
of relaxation rates in charge qubits reported by jAstafiev] 
et a/.[([2004[). 



Studies of flux noise have a long history. Already in 



(Koch et al., 1983) it was demonstrated that it is the 



flux and not the critical current noise that limits the 
sensitivity of dc SQUIDs. The interest in this prob- 
lem was recently renewed when it was realized that flux 
noise can limit the coherence in flux and phase super- 



conducting qubits ( Harris et al. , 2008 Yoshihara et al. 



2006 ) . Two recent models for fluctuators producing low- 



frequency noise were suggested. In (Koch et al.\ 2007) 



flux noise is attributed to the electron hopping between 
traps in which their spins have fixed, random orienta- 



tions. In (de Sousa, 2007) it is proposed that electrons 



flip their spins due to interaction with tunneling two- 
level systems (TLSs) and phonons. A novel mechanism 
- based on spin diffusion along the surface of a super- 
conductor - was suggested by [Faoro and loffe (2008). It 
seems to agree with recent experiments on measurements 



of the 1// flux noise reported by Bialczak et al. (2007) 
and[Sendelbach et a/.[([2008[). 



The microscopic mechanism and the source of the fluc- 
tuations of the critical current in a Josephson junction are 
long-standing open problems. These fluctuations were 
initially attributed to the charges tunneling or hopping 
between different localized states in the barrier forming 
glass-like TLSs. However, a more detailed comparison 
with experiments revealed an important problem: the 
noise spectrum experimentally observed by [Van Harlin 



[gen et a/.[ (|2004) and [Wellstood et fl/.[ ( [2004[ ) was propor- 
tional to T^, which is incompatible with the assumption 
of constant TLS density of states. The recent experi- 
ments by Eroms et al. (2006) on fluctuations in the nor- 
mal state in small Al junctions - similar to those used in 



et 



several types of qubits ( Chiorescu et al. 2003 Martinis 



2002 Vion et al. , 2002) - brought a new puzzle. It 



turned out that the temperature dependence of the noise 
power spectrum in the normal state is linear, and the 
noise power is much less than that reported for large su- 
perconducting contacts. A plausible explanation of such 



behavior was proposed by Faoro and loffe (2007), who 



leading to predictions in agreement with the experimen- 



suggested that the critical current noise is due to electron 
trapping in shallow subgap states that might be formed 
at the superconductor-insulator boundary. 

The dramatic effect of 1// charge noise was pointed 
out already in the breakthrough experiment performed 
at NEC and reported by [Nakamura et al\ (p^999 ) , where 
time-domain coherent oscillations of a superconducting 
charge qubit were observed for the first time. This exper- 
iment has renewed the interest for charge noise and has 
stimulated a number of investigations aimed at explain- 
ing decoherence due to noise sources characterized by a 
1// power spectrum in some frequency range. The the- 
oretical understanding of decoherence in solid-state sin- 
gle qubit gates is presently quite well established. Quite 



often in the literature the statistics of the fluctuations 
of the qubit parameters displaying 1// spectrum is as- 
sumed to be Gaussian. This assumption is not a priori 
obvious. In order to discuss the applicability range of 
the Gaussian approximation as well as deviations from 
the Gaussian behavior in connection with the problem 
of qubit dephasing we will consider in details some of 
the above mentioned models where the qubit response to 
typical manipulation protocols can be solved exactly. 

Peculiar features originated from fluctuations with 1// 
spectrum show up in the qubit dynamics superposed to 
effects due to high frequency fluctuations. The intrin- 
sic high-frequency cut-off of 1// noise is in fact hardly 
detectable, measurements typically extending to 100 Hz. 
Recently, charge noise up to 10 MHz has been detected 



et a/., 2002) and it is now applied to all types of su- 



perconducting qubits (except for phase qubits). Further 
substantial improvement resulted from the use of charge- 



or flux-echo techniques (Bertet et al. 2005 



et al. 



Nakamura 



2002 ) . In NMR the spin-echo removes the inhomo- 



geneous broadening that is associated with, for example, 
variations of static local magnetic fleld over the sample, 
changing the NMR frequency. In the case of qubits, the 
variation is in their energy-level splitting frequency from 
measurement to measurement. For some qubits defocus- 
ing is strongly suppressed by combining optimal point 
and echo techniques, thus providing further evidence of 
the fundamental role of 1// noise. Dynamical decou- 
pling, which uses sequences of spin flips to average out 
effectively the coupling to the environment, is another 



in a single-electron transistor (SET) by Kafanov et 

(|2008|) and flux noise in the 0.2 - 20 MHz range by |2004[|Faoro and Viola[|2004[). Optimized sequences limit 



proper pulse sequences has been measured by |Bylander| 
et al. (2011). Incoherent energy exchanges between sys- 



promising strategy (Bylander et al. 2011 



Falci et al. 



tem and environment, leading to relaxation and deco- 
herence, occur at typical operating frequencies (about 
10 GHz). Indirect measurements of noise spectrum in 
this frequency range (quantum noise) often suggest a 



"white" or Ohmic behavior (Astaflev et al., 2004 Ithier 



et al. [2005|. In addition, narrow resonances at selected 
frequencies (sometimes resonant with the nanodevice- 
relevant energy scales) have been observed ([Cooper et al. 



2004||Eroms et a/.[|2006]|Simmonds et a/.||2004D In cer- 



tain devices they originate from the circuitry (Van der 



Wal et al. 2000) and may eventually be reduced by cir- 



cuit design. More often, resonances are signatures of the 
presence of spurious fluctuators which also show up in 
the time resolved evolution, unambiguously proving the 



discrete nature of the noise sources (Duty et al., 2004). 



Such fluctuators may severely limit the reliability of nan- 
odevices (iFalci et al.\ |2005l iGalperin et al.\ |2006|. Ex- 



planation of this rich physics is beyond phenomenologi- 
cal theories describing the environment as a set of har- 
monic oscillators. On the other side, an accurate char- 
acterization of the noise sources might be a priori in- 
efficient, since a huge number of parameters would be 
required for a microscopic description. A roadmap to 
treat broadband noise which allows to obtain reasonable 
approximations by systematically including only the rel- 
evant information on the environment, out of the huge 
parametrization needed to specify it has been proposed 



in ( ,Falci et a/.[|2005| ). The obtained predictions for the 
decay of the coherent signal are in agreement with obser- 
vations in various superconducting implementations and 
in different protocols, like the decay of Ramsey fringes in 



charge-phase qubits (iVion et al., 2002). 



One successful strategy to increase phase-coherence 
times in the presence of 1// noise is to operate close to 
qubits' "optimal points" where low- frequency noise ef- 
fects vanish to the lowest-order. This strategy was flrst 
implemented in a device named "quantronium" (iVion, 



ing the blowup in resources involved in 1// noise suppres- 



sion ([Biercuk et a/.[[2QTT||Du et a/.|[2009||Lee et a/.||2008 



Uhrig, 2007) as well as "optimal control theory" design 



of quantum gates (Montangero et al. 2007 Rebentrost 



et al. , 2009 ) have been recently explored. The exper- 



imental realization of optimal dynamical decoupling in 
solid-state systems, and the implementation of quantum 
gates with integrated decoupling in a scalable and/or hy- 
brid architecture is an open problem. 



B. Why l//-noise is important for qubits? 

Let us briefly discuss in which way noise influences the 
operation of an elementary part of a quantum computer 
- a quantum bit (qubit). This topic will be addressed in 
details in Section |III| The qubit can be described as a 
two-level system with effective Hamiltonian 



Ha 



{cGz + Acr^) . 



(1) 



Here he represents diagonal splitting of the individual lev- 
els, h/S. represents their tunneling coupling, cjx^z are the 
Pauli matrices. The physical meaning of the quantities 
e and A depends on the speciflc implementation of the 
qubit. In its diagonalized form the qubit Hamiltonian ([l]) 
reads 



Ha 



hft 



hft 



■Cfz 



(cosO a^ + sin^cTa^) 



(2) 



where Q = \Je^ + A^ and the quantization axis a^' forms 
an angle d with a^. The Hamiltonian (IT]) corresponds to 
a pseudospin 1/2 in a "magnetic fleld" B, which can be 
time-dependent : 

H^ = {h/2)B a, B, = e(t), B, = A{t) . (3) 

Any state vector, |^), of the qubit determines the Bloch 
vector M through the density matrix 



|^)(^| = (l + M-t7)/2. 



(4) 



The Schrodinger equation turns out to be equivalent to 
the precession equation for the Bloch vector: 



the phase (j){t) accumulated by m+ during the time t as 
a sum of a regular, ^q, and stochastic, (p{t), parts: 



M = B xM. 



(5) 



^{t) = Mt) + ^W, Mt) = Bot, ^{t) 



The problem of decoherence arises when the "magnetic 
field" is a sum of a controlled part Bq and a fluctuating 
part h{t) which represents the noise, i. e., the field is 
a stochastic process, B(t) = Bq + b(t), determined by 
its statistical properties. The controlled part Bq is not 
purely static - to manipulate the qubit one has to apply 
certain high-frequency pulses of Bq in addition to the 
static fields applied between manipulation steps. In this 
language, the role of the environment is that it creates a 
stochastic field b(t), i. e., stochastic components of e{t) 
and A(t). These contributions destroy coherent evolution 
of the qubit making its coherence time finite. 

Let us consider first the simple case of "longitudi- 
nal noise", where b || Bq, and let the z-axis lie along 
the common direction of Bq and b, see Fig. [l] In the 
physics of magnetic resonance, this situation is called 
pure dephasing because z-component of the Bloch vec- 
tor, M^, is conserved during the process. As long as the 
time evolution of M is governed by Eq. ([s]), the length 
|M| = ^/M^^^M^ = M is also conserved, while the 
length I (M) | of the vector M averaged over the stochas- 
tic process h{t) decays. Description of this decay is the 
main objective of the decoherence theory. In the case of 
pure dephasing, this will be the decay of the components 
Mx and My. It is convenient to introduce a complex 
combination r?i+ = {Mx -\-iMy)/M. Equation ([5| can be 
written in terms of m+ as m+ = iBruj^ with solution 



/ h(t')dt' 
Jo 



and obtain: (m+(t)) = e^^o (e^^)m+(0). The stochas- 
tic phase (p is the integral of the random process b{t). 
The Bloch vector precesses around the 2:-axis with the 
angular velocity that has random modulation b{t). In 
the Gaussian approximation the only relevant statis- 
tical characteristics of b{t) is the correlation function 
{b{ti)b{t2)) = Sb{\ti — t2\) (we assume that b{t) is a sta- 
tionary random process). The function Sb{t) decays at 
\t\ ^ oo and the scale of this decay is the correlation time. 
If the integration time, t, strongly exceeds the correlation 
time the random phase (p is a sum of many uncorrelated 
contributions. According to the central limiting theorem 
such a sum has Gaussian distribution. 



^^"^ = 7^'^K-^ 



(7) 



independently of the details of the process. Therefore, 
the Gaussian distribution should be valid as soon as t 
exceeds the correlation time of the noise. In Sec. lIII.A.l.al 
we will further discuss this conclusion. As follows from 
Eq. 0, 



m+{t) = e'*(*)m+(0), (l){t) 



(e^^) = I p{<^)e'^ dip = e-^'^')/^ , (8) 

{^^)= f dh I dt2Sb{\h-t2\). (9) 

Jo Jo 



^ ^ ' ^ ^ Representing Sb{r) by its Fourier transform. 



Here we have assumed that the variable B{t) is a clas- 
sical one, i. e., the order of times in the products 
B{ti)B{t2) ■ ■ ■ B{tn) is not important. The solution has 
to be averaged over the stochastic process b{t). We define 



M, 



X 



Sb{uj) 



dt Sb{t) cos ujt , 



(10) 



that is just the noise spectrum^ and using Eq. ([9| we 
obtain 





M 


/ 


M:, 


\( 


/ 



(:/) = 2 



iy{ 



sm.uot/2 
uj/2 



Si,{uj). (11) 



Therefore, the signal decay given by Eq. ([8| is de- 
termined (in the Gaussian approximation) only by 
the noise spectrum, Sb{uj). For large t, the iden- 
tity liiRa^ooisiii^ ax/nax'^) = S{x) implies that {(p'^) = 
27rtSb{0) and thus 



'± 



FIG. 1 Bloch vector representing the state of a qubit in the 
rotating (with angular frequency Bq) frame of reference. In 
the laboratory frame of reference it precesses around the z- 
axis with (time-dependent) angular velocity B. Fluctuations 
of this velocity is b(t) 



{e'n 



= e-t/T^ 



nSbiO) . 



(12) 



Thus, the Gaussian approximation leads to exponential 
decay of the signal at large times, the decrement being 
given by the noise power at zero frequency. The above 
expression shows that the pure dephasing is determined 
by the noise spectrum at low frequencies. In particular. 



at Sb{uj) ex 1/uj the integral in Eq. ( pT] ) diverges, cf. with 
further discussion of spin echo. 

The time dependence of (m+) ex (e*"^) characterizes 



decay of the so-cahed free induction signal ( Schlichterj 
1992). The free induction decay (FID) is the observable 



NMR signal generated by non-equilibrium nuclear spin 
magnetization precessing about the magnetic field. This 
non-equilibrium magnetization can be induced, generally 
by applying a pulse of resonant radio-frequency close to 
the Larmor frequency of the nuclear spins. In order to ex- 
tract it in qubit experiments, one usually has to average 
over many repetitions of the same qubit operation. Even 



in setups that allow single-shot measurements (Astafiev 



et al 



2004 ) each repetition gives one of the two qubit 
states as the outcome. Only by averaging over many 
repeated runs can one see the decay of the average as de- 
scribed by the free induction signal. The problem with 
this is that the environment has time to change its state 
between the repetitions, and thus we average not only 
over the stochastic dynamics of the environment during 
the time evolution of the qubit, but over the initial states 
of the environment as well. As a result, the free induc- 
tion signal decays even if the environment is too slow 
to rearrange during the operation time. This is an ana- 
logue of the inhomogeneous broadening of spectral lines 
in magnetic resonance experiments. This analogy also 
suggests ways to eliminate the suppression of the signal 
by the dispersion of the initial conditions. One can use 
the well-known echo technique [see, e. g., (Mims, 1972)] nential decay with time scales 



Along the simplified model discussed above, the com- 
ponent M^ of the magnetic moment which is parallel to 
the "magnetic field" B does not decay because longitu- 
dinal fiuctuations of the magnetic field do not influence 
its dynamics. In a realistic situation it also decays in 
time, the decay time being referred to as relaxation time 
and conventionally denoted Ti. This occurs when the 
stochastic field b(t) has a component perpendicular to 
the controlled part Bq. Relaxation processes also induce 
another decay channel for the M^ and My components 
on a time scale denoted T2. Let us suppose that stochas- 
tic field is h{t) = b{t)z, while Bq = Ax + ez. Due to 
the non-isotropic interaction term, the effect of noise on 
the qubit phase-coherent dynamics depends on the angle 
0. When ^ = the interaction is longitudinal and we 
obtain the already discussed pure dephasing condition. 
When 0^0 the stochastic field also induces transitions 
between the qubit eigenstates. As a result, if the noise 
has spectral component at the qubit energy splitting, 
this interaction induces inelastic transitions between the 
qubit eigenstates, i. e., incoherent emission and absorp- 
tion processes. This is easily illustrated for a weak ampli- 
tude noise treated in the Markovian approximation. Ap- 
proaches developped in different areas of physics, as the 
Bloch-Redfield theory (|Bloch[ [T957| |Redfield| [T957t, the 



Born-Markov Master Equation ( Cohen- Tannoudji et al 



1992), and the systematic weak-damping approximation 



in a path-integral approach (Weiss, 2008), lead to expo- 



when the system is subjected to a short manipulation 
pulse (the so-called 7r-pulse) with duration ri at the time 
ri2. The duration ri of the pulse is chosen to be such that 
it switches the two states of the qubit. This is equivalent 
to reversing the direction of the Bloch vector and thus 
effectively reversing the time evolution after the pulse 
as compared with the initial one. As a result, the ef- 
fect of any static field is canceled and decay of the echo 
signal is determined only by the dynamics of the envi- 
ronment. The d ecay of two-p ulse echo can be expressed 



1/Ti = TT sin^ Sb{n) , 

I/T2 = I/2T1 + 1/T* 



(15) 
(16) 



as (m!f^(2ri2)) ( |Mims 



1972), where 



(mi^)(t))^(e^^(^)), m 



/ rTi2 pt 
\Jo Jti: 



h{t')dt'. (13) 



where the adiabatic or pure dephasing term in the gen- 
eral case reads I/T2* = ttcos^ ^£'5(0). Even if the above 
formulas do not hold for 1// noise, which would lead to 
a singular dephasing time, they indicate that the diverg- 
ing adiabatic term containing Sh{^) may be eliminated 
(in lowest order) if the qubit operates at ^ = 7r/2, or 
equivalently when the noise is "transverse", b ± Bq. In 
the context of quantum computing with superconducting 
systems, this condition - of reduced sensitivity to 1// 
noise - is usually referred to as "optimal point" (|Vion| 



et al. , 2002 ) . In this case T2 = 2Ti , which is the upper 



Finite correlation time of h(t) again leads to the Gaussian 
distribution of '^(t) at large enough t with 



(V''(2ti2)) 



Jo 



dw 



uj/2 ) 



Sb{uj) . (14) 



This variance can be much smaller than {^p^) given by 
Eq. (11) \i Shioo) is singular at a; ^ 0. Though the in- 
tegral (14) is not divergent in the case of l//-noise, the 
time dependence of the echo signal is sensitive to the low- 
frequency behavior of Sh- Therefore, the low- frequency 
noise strongly affects coherent properties of qubit s. 



limit to the dephasing time. We will come again to this 
issue in Sec. IIII.Bi 



II. PHYSICAL ORIGIN OF 1// NOISE IN NANODEVICES 
A. Basic models for the 1// noise 

Studies of the noise with the spectral density (x. uj~^ 



have a long history, see, e. g., (Kogan, 1996) for a review. 



This behavior at low frequencies is typical for many phys- 
ical systems, such as bulk semiconductors, normal metals 
and superconductors, strongly disordered conductors as 



well as devices based on these materials. One observes, 
in practically all cases, an increase of the spectral density 
of the noise with decreasing frequency / = cj/27r approx- 
imately proportional to 1// down to the lowest experi- 
mentally achievable frequencies. Therefore, the noise of 
this type is referred to as 1// noise (the term 'flicker 
noise' proposed by Schottky is now more rarely used). 

1// noise poses many puzzles. Is the noise spectrum 
decreasing infinitely as / -^ 0, or does it saturate at 
small frequencies? What are the sources of the 1// noise? 
Why many systems have very similar noise spectra at 
low frequencies? Is this phenomenon universal and does 
a unified theory exist? Those and many other questions 
stimulated interest to the 1// noise from a fundamental 
point of view. This interest is also supported by the 
crucial importance of such a noise for all the applications 
based on dc and low- frequency response. 

Usually, the observed 1// noise of the electrical cur- 
rent is a quadratic function of the applied voltage in uni- 
form Ohmic conductors. This indicates that the noise is 
caused by fiuctuation in the sample resistance, which are 
independent of the mean current. Therefore, the current 
just 'reveals' the fluctuations. Though the low-frequency 
noise spectrum seems to be rather universal the noise in- 
tensity differs very much in different systems depending 
not only on the material, but also on the preparation 
technology, heat treatment, etc. These facts lead to the 
conclusion that in many cases the 1// noise is extrinsic^ 
i. e., caused by some dynamic defects. 

In the simplest case, when the kinetics of fluctuations 
is characterized by a single relaxation rate, 7, the corre- 
lation function of a fluctuating quantity x{t) is propor- 
tional to e~^l^l . Then the spectral density is a Lorentzian 
function of frequency. 



Sx{uj) oc C^{uj) 



1 



7 



TT CJ^ + 7' 



2 • 



(17) 



In more general cases, the kinetics of x{t) is a superpo- 
sition of several, or even many, relaxation processes with 
different rates. In general, a continuous distribution of 
the relaxation rates, Vxi'y) may exist. Then the noise 
spectral density can be expressed as 



Sx{uj) (X / djVx{7)jC^{oj). 



(18) 



Equation (18) describes a non-exponential kinetics. 



Vx (7) dj being the contribution of the processes with re- 
laxation rates between 7 and 7 + ^7 to the variance 



{Sxy 



dujSx{uj)(x / djVxil). (19) 

Jo 



If T^x{l) ex 7 ^ in some window 70 ^ 7min, but very 



small outside this interval, then, according to Eq. (18), 
one gets ^'^^(ci;) oc uj~^ in the frequency domain 70 ^ 
^ ^ 7min (Surdin 1939). Following a model of this 




FIG. 2 Schematic diagram of the potential energy of a two- 
level tunneling system versus configuration coordinate. U is 
the distance between the single- well energy levels, character- 
izing asymmetry of the potential. E is the difference between 
the lowest energy levels in a two-well potential with the ac- 
count of tunneling below the barrier. 



type, one has to specify the processes responsible for the 
noise, which depend on the speciflc properties of the sys- 
tem under consideration. Several types of kinetic pro- 
cesses were indicated as able to produce the 1// noise. 
Among them are activated processes with different relax- 
ation rates exponentially dependent on the inverse tem- 
perature: 7 = 7oe~^/^^^, where 70 is some attempt fre- 
quency and the distribution of activation energies, J-{E), 
is smooth in a sufficiently broad domain. Then the distri- 
bution of the relaxation rates has the required fo rm since 
V.h) = :Fx{E)\d^/dE\-^ = Tx{E){kBT/j) dvan der 



Ziel 1950 



If the kinetics of the ffuctuations is controlled by tun- 
neling processes then the relaxation rates depend approx- 
imately exponentially on the width and the height of the 
tunneling barrier. If the distribution of these parame- 
ters is almost constant in a wide interval, than again 
the distribution of relaxation rates will be proportional 
to 1/7. In particular, iMcWhorter ( 1957| ) suggested that 
fluctuation of number of carriers in a surface layer of a 
semiconductor arise from exchange of electrons between 
this layer and the traps lying in the oxide layer covering 
the surface, or on the outer surface of the oxide. Since 
the electron transfer takes place via tunneling, the char- 
acteristic relaxation rate exponentially depends on the 
distance x between the surface and trap: 7 = 706"^/^^. 
Since the distances x vary with a scatter :^ A, the dis- 
tribution of relaxation times is exponentially wide. This 
model has been extensively used to interpret 1// noise 
in field-effect transistors. However, the model failed to 
explain observed very small rates 7 in the devices with 
relatively thin insulator layer. 

There exist special low-energy excitations in amor- 
phous materials, in particular, in all dielectric and metal- 
lic glasses. They result in anomalous temperature depen- 
dences of the heat capacity and the thermal conductivity 
at low temperatures, as well as specific features of the 



sound absorption, see, e. g., (Black 1981 Galperin et al. 



1989 Hunklinger and von Schickfus, 1981) for a review. 



According to this model, atoms or groups of atoms exists 
that can occupy two positions. Therefore, their energy as 
a function of some configuration coordinate can be rep- 
resented as a double-well potential, as shown in Fig. [2J 



The model of two-level systems (TLSs) formulated by 
Anderson et al. ( 1972, ) and |Pliillips| ( |1972| ) contains two 
parameters - the asymmetry, [/, of the potential (which is 
approximately equal to the difference between the min- 
ima) and the tunnel matrix element, A, characterizing 
the strength of the barrier which can be estimated as 



A = huj^ 



'oe 



(20) 



where ujq is the frequency of the intra- well vibrations. 
Thus the inter-level spacing £^ of a TLS (which is its 
excitation energy) is given by 



E = VU^ + A^ . 



(21) 



The transitions of atoms or group of atoms between 
these levels, and the change in the levels' relative occu- 
pancy with varying temperature or under acoustic vibra- 
tions are responsible for the low-temperature properties 
of structural glasses. The same type of low-energy exci- 
tations were also found in amorphous metals and ionic 
conductors, see (Black , 1981) for a review. 

Due to disorder, the TLSs have different values of U 
and A. Physical considerations lead one to assume that 
the distribution of U and A, 7^(/7, A), is almost constant 
in the region A ^ 1, /7 <C ficjo important for the effects 
observed in the experiments. Therefore it is assumed that 
in this region ViU^X) = Pq. Here Pq is a constant, which 



can be found by comparison with experiments (Black 



1978 Halperin 1976) 



The rate of transitions between the two levels of a TLS 
is determined by interactions with phonons (in insulat- 
ing solids), or with electrons (in metals). Assuming that 
fluctuations of the diagonal splitting U are most impor- 
tant, we can describe the interaction between the TLS 
and the environment as 



'WxLS-env = 9 CTz 



(22) 



where c is an operator in the Hilbert space of the environ- 
ment depending on the specific interaction mechanism. It 
is convenient to diagonalize the TLS Hamiltonian, 

HTLS = ^(/7r,+Ar,) , (23) 



where r^ are the Pauli matrices, by rotating the TLS 
Hilbert space. Then 

Htls = {E/2)t, , 



-env — 9 (^ 



E 



E 



(24) 
(25) 



The interlevel transitions are described by the second 



item in the interaction Hamiltonian (25). Therefore, the 



relaxation rate for the deviation of the occupancy num- 
bers of the level s from the equilibr i um ones is propor - 
tionalto {K/Ef ([Black and Gyorffy[ [19781 |Jackle[|l972[): 



7 = 7o(£^) 



7o(^)^^^coth(^^) . (26) 



The quantity ^q{E) has the meaning of a maximal relax- 
ation rate for the TLSs with given inter-level spacing, E. 
The exponent a depends on the details of the interaction 
mechanism, its typical values are 3 (for the interaction 
with phonons) and 1 (for the interaction with electrons) . 
Using Eq. (26) one can obtain the distribution of the 



relaxation rates as 



V{l,E) 



E 



nu,x) 



Po 



27\/l - 7/70 



27 



(27) 



(here we have taken into account that small relaxation 
rates require small tunnel coupling of the wells, A <C 
U c^ E). Therefore, owing to the exponential depen- 
dence of 7 on the tunneling parameter A (7 oc e"^'^) and 
approximately uniform distribution in A, the distribution 
with respect to 7 is inversely proportional to 7 in an ex- 
ponentially broad interval, as is characteristic of systems 
showing 1// noise. 

Spontaneous transitions between the levels of the TLSs 
can lead to fluctuations of macroscopic properties, such 



as resistance of disordered metals ( ]Kogan and Nagaev 



1984b Ludviksson et al , 1984), density of electron states 



in semiconductors and metal-oxide-semiconductor struc- 



tures (Kogan and Nagaev 1984a), etc. These fluctua- 



tions have 1// spectrum. 

Two-level systems can be formed in hopping insulators 
by different electronic configurations having close ener- 
gies (|Burin£raL 2006). Giant relaxation times necessary 
for 1// noise are provided by a slow rate of simultaneous 
tunneling of many localized electrons and by large activa- 
tion barriers for their consecutive rearrangements. The 
model qualitatively agrees with the low-temperature ob- 
servations of 1// noise in p-type silicon and GaAs. Sev- 
eral other models of the 1// noise are briefly reviewed in 
the book by |Kogan| ( [19961 ) . 



B. 1// noise and random telegraph noise 

In many systems comprising such semiconductor de- 
vices as p — n junctions, metal-oxide-semiconductor field- 
effect transistors, point contacts and small tunnel junc- 
tions between the metals, small semiconductor resistors 
or small metallic samples, the resistance switches at ran- 
dom between two (or several) discrete values. The time 
intervals between switchings are random, but the two 
values of the fluctuating quantity are time-independent. 
This kind of noise is now usually called random telegraph 



noise (RTN), see (Kogan 1996[ ) for a review. 

Statistical properties of RTNs in different physical sys- 
tems are rather common. Firstly, the times spent by the 
device in each of the states are much longer than the mi- 
croscopic relaxation times. Therefore, the memory of the 
previous state of the system is erased, and the random 
process can be considered as a discrete Markovian pro- 
cess. Such a process (for a two-state system) is character- 
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ized by the equilibrium probabilities pi and P2 = ^ — Pi 
of finding the system respectively in the first or in the 
second state, as well as transition probabilities per unit 
time, 71^2 and 72^1. The spectral density of a random 
quantity, which switches between the two states xi and 
X2 can be easily derived as, cf. with ( Kogan[| 1996 1), 
(xi -^2)^ 



Sx{uj) 



-C^{oj) 



4cosh^{E/2kBT) 

where 7 = 71^2 + 72^1 • Therefore each RT process 
contributes a Lorentzian line to the noise spectrum. 

RTN was observed in numerous small-size devices. At 
low temperatures, usually only one or few telegraph pro- 
cesses were observed. However, at higher temperatures 
(or at higher voltages applied to the device) the number 
of contributing telegraph processes increased. Typically, 
at high enough temperatures discrete resistance switch- 
ing is not observed. Instead a continuous 1// noise is 
measured. This behavior was interpreted by [Rogers and 



[Buhrman] (1984, 4985) as a superposition of many un- 
correlated two-state telegraph fluctuators with various 
relaxation rates 7. The increase of temperature leads to 
an increase of the number of contributing fluctuators and 
discrete switchings become indistinguishable. 

A different interpretation was suggested by [Ralls and 



Buhrman 



(1991), see also references therein. This in- 



terpretation is based on interaction between the fluctu- 
ators leading to a deviation from the simple Lorentzian 
spectrum. Moreover, the system of interacting defects 
may pass to another metastable state where different 
defects play role of active fluctuators. This interpre- 
tation is based on the observation of RTN in metallic 
nanobridges where the record of resistance at room tem- 
perature is still composed by one or two telegraph pro- 
cesses, but their amplitudes and characteristic switching 
rates change randomly in time. According to this in- 
terpretation, the systems of interacting dynamic defects 
is similar to a glass, particularly, to a spin glass with a 
great number of metastable states between which it is 
incessantly wandering. 

Despite great progress of the 1// noise physics, for 
the major part of systems showing 1// or/and random 
telegraph noise the actual sources of the low-frequency 
fluctuations remain unknown: this is the main unsolved 
problem. Below we will discuss some simple models in 
connection with devices for quantum computation. 

C. Superconducting qubits and relevant noise mechanism^ 

Circuits presently being explored combine in variable 
ratios the Josephson effect and single Cooper-pair charg- 



^ More detailed descri p tion can be foun d in th e revie ws by [Clarke 



I and Willielm[ ( |2008| >; [Makhlin et al\ ( [2001| >; and [Xiang et cd. 
( |2Q13| >, 



ing effects. When the Coulomb energy is dominant, the 
"charge circuits" can decohere from charge noise gen- 
erated by the random motion of offset charges (jNaka- 



mura et al\ [20021 [Zimmerli eFlI] [19921 [Zorin et al 



19961). Conversely, when the Josephson energy is domi- 
nant, these "flux circuits" are sensitive to external flux 



and its noise (Friedman et al. 2000 Mooij et al, 1999 



(28) Wellstood et a/., 1987b). For the intermediate-energy 



regime, a circuit designed to be insensitive to both the 
charge and flux bias has recently achieved long coher- 
ence times (< 500 ns), demonstrating the potential of 



superconducting circuits ( Cottet et al. , 2002 Vion et al. 



2002). The third type is the phase qubit, which consists 
of a single Josephson junction current biased in the zero 

In 



voltage state (Martinis et al. 2002 



Yu et al. 



2002) 



this case, the two quantum states are the ground and 
first excited states of the tilted potential well, between 
which Rabi oscillations have been observed. 

In the case of charge qubits, the coherence times have 
been limited by low- frequency fluctuations of background 
charges in the substrate which couple capacitively to the 
island, thus dephasing the quantum state (Nakamura 



et al. 2002). Flux and phase qubits are essentially im- 



mune to fluctuations of charge in the substrate, and, by 
careful design and shielding, can also be made insensitive 
to flux noise generated by either the motion of vortices in 
the superconducting films or by external magnetic noise. 
The flux-charge hybrid, operating at a proper working 
point, is intrinsically immune to both charge and flux 
fluctuations. However, all of these qubits remain sensi- 
tive to fluctuations in the Josephson coupling energy and 
hence in the critical current of the tunnel junctions at low 
frequency /. These fluctuations lead to variations in the 
level splitting frequency over the course of the measure- 
ment and hence to dephasing. 



1. Charge noise in Josephson qubits 

The importance of the charge noise was recognized af- 
ter careful spin-echo-type experiments applied to an arti- 
ficial TLS utilizing a charge degree of freedom of a small 
superconducting electrode - so-called single- Cooper-pair 



box (CPB) ( [Bouchiat e^aq[1998| [Nakamura et aq|1997 ). 
To explain the main principle behind this device, let 
us consider a small superconductor grain located close 
a (gate) metallic electrode. The ground state energy of 
such a grain depends in an essential way on the number of 
electrons on it. Two contributions to such a dependence 
are given by the electrostatic Coulomb energy Ec{n) de- 
termined by the extra charge accumulated on the super- 
conducting grain and by the so-called parity term A^ 



([Files et al.\ |1993t jGlazman et al. 


|1994t [Hekking et al.\ 


19931 [Hergenrother et al.\ [19941 J 


oyez et al.\ [1994| [La- 


farge et al.\ 1993; Matveev et al.^ 1993 1994 Tuominen 


et al.\ [1992, p93,). The latter or 


iginates from the fact 
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FIG. 3 Left panel: The energy diagram for the ground state 
of a superconducting grain with respect to charge for the case 
Ao > Ec- For a certain bias voltage when q = (2n + l)e, 
ground states differing by only one single Cooper pair become 
degenerate. Right panel: Energy of a Cooper-pair box with 
account of the Josephson tunneling. 



that only an even number of electrons can form a BCS 
ground state of a superconductor (which is a condensate 
of paired electrons) and therefore in the case of an odd 
number of electrons n one unpaired electron should oc- 



cupy one of the quasiparticle states (Averin and Nazarov 
[1992 J p] 

The energy cost of occupying a quasiparticle state, 
which is equal to the superconducting gap Aq, brings 
a new scale to bear on the number of electrons that a 
small superconducting grain can hold. Taking the above 
into account, one presents the ground state energy Uo{n) 
in the form, see (Glazman et al.^ ^1994^ Hekking et al\ 



1993 Matveev et al 1993 1994), 



[/o(n) = Ec(n- ^)'= A,, ^n = [\^ 



even n. 



odd 



n. 



(29) 



Here q is the charge induced on the grain by the gate 
electrode. One can see from Eq. (29) that if Aq > Ec 



only an even number of electrons can be accumulated in 
the ground state of the superconducting grain. Moreover, 
for special values of the gate voltage corresponding to 
q = (2n + l)e a degeneracy of the ground state occurs 
with respect to changing the total number of electrons by 
one single Cooper pair. An energy diagram illustrating 
this case is presented in Fig. [3) The occurrence of such 
a degeneracy brings about an important opportunity to 
create a quantum hybrid state at low temperatures which 
will be a coherent mixture of two ground states, differing 
by a single Cooper pair: 

|^)=ai|2n)+a2|2(n + l)). (30) 

The idea of the device is presented in Fig. |4] where the 



Reference 7 in JAverin and Nazarov[ |l992| indicates that the 
original idea belongs to K. A. Matveev. 
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FIG. 4 A schematic diagram of a single-Cooper-pair box. An 
island of superconducting material is connected to a larger 
superconducting lead via a weak link. This allows coherent 
tunneling of Cooper pairs between them. For a nanoscale sys- 
tem, such quantum fluctuations of the charge on the island are 
generally suppressed due to the strong charging energy asso- 
ciated with a small grain capacitance. However, by appropri- 
ate biasing of the gate electrode it is possible to make the two 
states, |2n) and |2(n+ 1)), differing by one Cooper pair, have 
the same energy (degeneracy of the ground state). This allows 
the creation of a hybrid state |^) = ai\2n) + a2|2(n + 1)). 



superconducting dot is shown to be in tunneling con- 
tact with a bulk superconductor. A gate electrode is re- 
sponsible for lifting the Coulomb blockade of Cooper-pair 
tunnelling (by creating the ground state degeneracy dis- 
cussed above) . This allows the delocalization of a single 
Cooper pair between two superconductors. Such a hy- 
bridization results in a certain charge transfer between 
the bulk superconductor and the grain. At the charge 
degeneracy point, q = (2n + l)e, the Josephson tunnel- 
ing produces an avoided crossing between the degenerate 
lavels corresponding to the symmetrical and antisymmet- 
rical superpositions, |2n) ± |2(n + 1)). As a result, the 
terms are split by an energy Ej <C Ec- Far from this 
point the eigenstates are very close to being charge states. 
To summarize, in a single-Cooper-pair box all electrons 
form Cooper pairs and condense in a single macroscopic 
ground state separated from the quasiparticle states by 
the superconducting gap, Aq. The only low-energy exci- 
tations are transitions between the charge number states, 
|2n), which are the states with excess number of Cooper 
pairs in the box due to Cooper-pairs tunneling if Aq is 
larger than the single-electron charging energy of the box, 
Ec' The fluctuations of n are strongly suppressed if Ec 
exceeds both the Josephson energy, Ej, and the thermal 
energy, ksT. Then we come back to the Hamiltonian 
^ with he = Ec{n — q/e^ (where q is the induced 
charge) and fiA = Ej^ i. e., the Josephson energy of 
the split Josephson junction between the box and the su- 
perconducting reservoir, see Fig. [4J Therefore, e can be 
tuned through the gate voltage determining the induced 
charge. The Josephson junction is usually replaced by a 
dc SQUID with low inductance. Ej (and consequently, 
A) is then adjusted by applying the appropriate mag- 
netic flux. A realistic device is shown Fig. |5] adapted 



from ([Nakamura et al\ 1999p. 



Nakamura et al. (2002) have compared the decay of 
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FIG. 5 Single-Cooper-pair box with a probe junction - Mi- 
crograph of the sample. The electrodes were fabricated by 
electron-beam lithography and shadow evaporation of Al on 
a SiNcc insulating layer (400-nm thick) above a gold ground 
plane (100-nm thick) on the oxidized Si substrate. The 'box' 
electrode is a 700 x 50 x 15 nm Al strip containing ^10^ con- 
duction electrons . Adapted by permiss ion from Macmillan 
Publishers Ltd.: ([Nakamura et a/.|p^99t, copyright (1999). 



the normalized echo signal, Fig. [6J with the expression 
by|Cottet et a/.|(|200lD 



(e^^) = exp 



1 
"2^2 



dujSe{oj) 



sin^(cL;r/4) 



(31) 



where Se{uj) is the spectrum of the noise in the interlevel 
spacing, fie, of the qubit. The latter is expressed through 
the charge noise, Se{uj) as Se{uj) = {4Ec/e)'^Se{uj). The 
charge noise spectrum was determined by a standard 
noise measurement on the same device used as a single- 
electron transistor. It can be expressed as Se{oj) = a/uj 
with a = (1.3 X lO^^e)^. The estimate following from 
Eq. ( [31] ), with the mentioned value of a and cjmin = 
27r/tmax being the low-frequency cutoff due to the finite 
data-acquisition time tmax (20 ms), is shown in Fig. Ig] 
(solid line 1). Solid line 2 in the same figure corresponds 
toa = (3.0 X 10"'^e)^ Note that Eq. ^ is based on the 
assumption that the fluctuations are Gaussian; it pre- 
dicts that at small delay time, r, the echo signal decays 




FIG. 6 Decay of the normalized amplitude of the echo signal 
(filled circles) and the free induction decay signal (open cir- 
cles) compared with estimated decoherence factors (e*'^) due 
to charge noise with the spectrum a/u. Here ^/ a x 10^e~ ^ 
is 1.3 for line 1 and 2 — 0.3 for 2. Adapted from (Nakamura 



et al. 



2002). 




FIG. 7 Energy relaxation rate, Fi (closed circles and open 
squares), and phase decoherence rate F2 = T^ ^ (open circles) 
versus gate induced charge q. Adapted from ( Astafiev et al.\ 
20041 



as ln(e*'^) ex — r^. In Section III this assumption and the 



ensuing prediction will be further discussed. 

(2004) studied decoherence of the 



Astafiev et al. 



Josephson charge qubit by measuring energy relaxation 
and dephasing with help of single-shot readout. Both 
quantities were measured at different charges induced at 
the single-Cooper-pair box by the gate electrode. The 
decoherence was determined from decay of the coherent 
oscillations related to the noise spectrum as 



ln(e^^) 



2[{he)^^EJ] 



duj Se{oj) 



sm{ujt/2) 
uj/2 



n2 



(32) 

Based on the dependence of the decoherence rate on the 
induced charge shown in Fig. [7| and on estimates of the 
noise, authors conclude that the source of decoherence is 
charge noise having 1// spectrum. Another conclusion 
is that the energy relaxation rate, Fi, is also determined 
by low-frequency noise of the same origin. This conclu- 
sion is drawn from the observed Fi oc Ej/[{he)'^ + Ej] 
dependence. The importance of 1// charge noise for de- 
coherence in the so-called quantronium quantum bit cir- 



cuit ( Cottet et al. 2002 ) was emphasized by Ithier et al. 
([2005^ 



To verify the hypothesis about the common origin of 
the low- frequency 1// noise and the quantum / noise re- 
cently measured in the Josephson charge qub its, [Astafiev 
\et a/.| (g2006) studied the temperature dependence of the 
1// noise amplitude and decay of coherent oscillations. 
The 1// noise was measured in the single-electron tunnel- 
ing (SET) regime. In the temperature domain 50 mK - 
1 K it demonstrated oc T^ dependence, see Fig. [s] adapted 
from (Astafiev et a/., 2006). The measurements of the 



noise were accompanied by measurements of the decay 
rate of the coherent oscillations away from the degen- 
eracy point {he ^ Ej). The decay of the oscillations 



was fitted according to Eq. (32) yielding at small times 
the dependence (e"*'^) ex e~* /^^2 . The results turned 
out to be consistent with the strength of the 1// noise 
observed in transport measurements. To explain the 
quadratic temperature dependence of the 1// noise the 
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FIG. 8 Temperature dependence of the amplitudes a^^^, 
which can be ap proximated as o}''^ — (1.0 x 10~^ e/K)T. 
Adapted from ([Astafiev ei al. 2006). 



authors assumed that this dependence originates from 
two-level fluctuators with the density of states linearly 
dependent on the inter-level spacing. We will discuss this 
assumption later while considering models for the noise- 
induced decoherence. Recently, measurements of charge 
noise in a SET shown a linear increase with tempera- 
ture (between 50 mK and 1.5 or 4 K) above a voltage- 
dependent threshold, with a low temperature saturation 



below 0.2 K (Gustafsson et a/., 2012). Authors conclude 



that this result is consistent with thermal interaction be- 
tween SET electrons and TLSs residing in the immedi- 
ate vicinity of the device. Such possible defects inlude 
residue of Al grains formed around the perimeter of the 
SET island and leads during the two-angle evaporation 
(Kafanov et al. , 2008), as well as interface states between 



metal of the SET and its surrounding oxides ( Choi et al 



2009). 



Though the obvious source of RT charge noise is a 
charge, which jumps between two different locations in 
space. Less clear is where these charges are actually lo- 
cated and what are the two states. The first attempt 
of constructing such a model in relation to qubit deco- 



herence appeared in (Paladino et al. 2002), where elec- 



trons tunneling between a localized state in the insulator 
and a metallic gate was studied. This model has been 



further studied in (Abel and Marquardt 2008 Grishin 



tal results (Astafiev et a/., 2004) indicated a linear de- 



et al. 20051 Yurkevich et a/., 2010). Later, experimen- 



pendence of the relaxation rate on the energy splitting 
of the two qubit states. One also has to take into ac- 
count that in the experimental setup there is no normal 
metal in the vicinity of the qubit: all gates and leads 
should be in the superconducting state at the temper- 
atures of experiment. These two facts suggest that the 




Model I 



Modern 



Model III 



FIG. 9 Three possible models for the fluctuating charges 
Model I, electrons jumping bet ween a localized state and a 
normal metal, as discussed in (Abel and Marquardt 2008 



Grishin et al. [2005 Paladino et al 



Model II, elec 
trons jumping between localized states. Model III, electrons 
jumping bet ween localized stat es and a superconductor, as 
discussed in (Faoro et al. 2005). 



A constant density of these localized states gives a lin- 
early increasing density of occupied pairs, in agreement 



with experiments (Astafiev et a/., 2004). This model was 



criticized (Faoro and loffe 2006) because it required an 



unreasonably high concentration of localized states, and 
a more elaborate model involving Kondo-like traps was 



proposed. However, it was shown (Abel and Marquardt 



2008 Grishin et a/., 2005) that allowance for quantum 



effects of hybridization between the electronic states lo- 
calized at the traps and extended states in the electrodes 
relaxes the above requirement. At present it seems that 
no solid conclusions can be drawn based on the available 
experiments. 

As we have seen, the "standard" Cooper-pair boxes 
are rather sensitive to low-frequency noise from electrons 
moving among defects. This problem can be partly re- 
laxed in more advanced charge qubit s, such as trans- 



mon (Koch et a/., 2007) and quantronium (Vion et al. 



2002). The transmon is a small Cooper-pair box where 



the Josephson junction is shunted by a large external 
capacitor to increase Ec and by increasing the gate ca- 
pacitor to the same size. The role of this shunt is played 
by a transmisison line, and therefore the qubit is called 
the transmon. The main idea is to increase the ratio 
Ej/Ec making the energy bands shown in the Fig. [s] 
(right panel) almost flat. For this reason, the transmon 
is weakly sensitive to low-frequency charge noise at all 
operating points. This eliminates the need for individ- 
ual electrostatic gate and tuning to a charge degeneracy 
point. A complementary proposal for using a capacitor to 
modify the Ej/Ec ratio in superconducting flux qubits 



model (Paladino et al. 2002) was not directly applicable is put forward by, You et al. (2007) 



to explain the decoherence in charge qubits (Nakamura 



et a/., 2002) and favored a model with superconduct- 
ing electrodes (Faoro et a/., 2005). In this model, the 
two electrons of a Cooper pair are split and tunnel sepa- 
rately to some localized states in the insulator (see Fig. |9] 
for an illustration of this (Model III) and other models) . 



At the same time, the large gate capacitor provides 
strong coupling to external microwaves even at the level 
of a single photon, greatly increasing coupling for cir- 
cuit quantum electrodynamics (cQED). The schematics 
of the suggested device is shown in Fig. [lO] reprinted from 
(Blais et al. 2004) where a detailed discussion of cQED 
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FIG. 10 (Color online) Schematic layout and equivalent 
lumped circuit representation of proposed implementation of 
cavity QED using superconducting circuits. The ID transmis- 
sion line resonator consists of a full- wave section of supercon- 
ducting coplanar waveguide, which may be lithographically 
fabricated using conventional optical lithography. A Cooper- 
pair box qubit is placed between the superconducting lines 
and is capacitively coupled to the center trace at a maxi- 
mum of the voltage standing wave, yielding a strong electric 
dipole interaction between the qubit and a single photon in 
the cavity. The box consists of two small (~ 100 x 100 nm^) 
Josephson junctions, configured in a ?^ 1 /xm loop to permit 
tuning of the effective Josephson energy by an external flux 
$ext- Input and output signals are coupled to the resonator, 
via the capacitive gaps in the center line, from 50 Q transmis- 
sion lines, which allow measurements of the amplitude and 
phase of the cavity transmission, and the introduction of dc 
and rf pulses to manipulate the qubit states. Reprinted from 
dBlais et a/._, 2004j . 



devices is presented. [Koch et at. (2007) have considered 



possible sources of dephasing and energy relaxation in 
transmon devices. In particular, the contributions of the 
charge, flux, and critical current to the dephasing rate, 
T2"^, were estimated. These (quite favorable) estimates 



are based on Eqs. (11) and (12), which follow from the 
assumption of the Gaussian statistics of the noise. In 
some cases, the decoherence is predicted to be limited 
by the energy relaxation, T2 ~ 2Ti, which occurs due 



to dielectric losses (Martinis et al. 2005) and quasipar- 



tide tunneling (Lutchin et al. 2005 2006). Catelani 



et al. (2011) developed a general theory for the qubit 
decay rate induced by quasiparticles. They studied its 
dependence on the magnetic flux used to tune the qubit 
properties in devices such as the phase and flux qubit s, 
the split transmon, and the fluxonium. Recently [Hassler 
\et al.\ (l2011| proposed to use a transmon qubit to perform 
parity-protected rotations and read-out of a topological 
qubit. The advantage over an earlier proposal using a 
flux qubit is that the coupling can be switched on and 
off with exponential accuracy, promising a reduced sen- 
sitivity to charge noise. 



Barends et al. (2013) have demonstrated a planar, tun- 



able superconducting qubit with energy relaxation times 
up to 44 /is. This was achieved by using a geometry based 



on a planar transmon (Houck et al. 2007 Koch et al. 



^ /aE}\ (m) 



Lcl 



^ 



/\^J 



EjX 



$ 



-^ 



$1 I—: 



1999). 



FIG. 11 (Color online) A superconducting loop with three 
Josephson junctions (indicated with crosses) encloses a flux 
$ that is supplied by an external magnet. Two junctions 
have a Josephson coupling energy Ej, and the third junction 
has aEj, where a = 0.75. This system has two (meta) stable 
states, |0) and |1), with opposite circulating persistent cur- 
rent. The level splitting is determined by the offset of the 
flux from $0/2. The barrier between the states depends on 
the value of a. The qubit is operated by resonant microwave 
modulation of the enclosed magnetic flux by a superconduct- 
ing control line (indicated in red). From (Mooij et al.\ 
Reprinted with permission from A A AS. 



reduce coupling to materials-related defects. The authors 
report a fine structure in the qubit energy lifetime as a 
function of frequency, indicating the presence of a sparse 
population of incoherent, weakly coupled two- level de- 
fects. The suggested qubit (called 'Xmon' because of its 
special geometry) combines facile fabrication, straight- 
forward connectivity, fast control, and long coherence. 

Finally, we mention an alternative single Cooper-pair 
circuit based on a superconducting loop coupled to an LC 
resonator used for dispersive measurement analogously 
to cQED qubit s and insensitive to offset charges. The 
circuit, named fluxonium, consists of a small junction 
shunted with the Josephson kinetic inductance of a se- 
ries array of large-capacitance tunnel junctions, thereby 
ensuring that all superconducting islands are connected 
to the circuit by at least one large junction (Manucharyan 



et al. , 2009). The array of Josephson junctions with ap- 



propriately chosen parameters can perform two functions 
simultaneously: short-circuit the offset charge variations 
of a small junction and protect the strong nonlinearity of 
its Josephson inductance from quantum fluctuations. 



2. Flux and phase qubits 



A flux qubit, see, e. g.. Fig. 11 adapted from (Mooij 



et al.^ 1999) 



2007) and designed to both minimize radiative loss and 



consists of a micrometer-sized loop with 
three or four Josephson junctions.. The energy of each 
Josephson junction can be expressed as Ej{l — cos 6i) 
where Ej is the Josephson energy of i-th junction while 
Si is the phase drop on the junction. The phase drops 
on the junctions are related as ^- di + 27r($/^o) = 27rn 
where n is an integer number. Here $ is the magnetic flux 
embedded in the loop while ^0 = Trhc/e is the magnetic 
flux quantum. To get the total energy one should add the 
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FIG. 12 (Color onine) Plot u{S) for a = 0.75. The inset 
show the potential profile along the line A for different /3 = 
(2$ - $o)/$o 



charging energy of each junction, (5|/2(7i, and magnetic 
energy, {^ — <I>5)^/2L where ^5 is the bias flux created 
by external sources. Neglecting charging energies and 
magnetic energy of the loop end assuming that Ej = 



E 



(2) 



Ej^E^ 



(3) 



aEj one can express the total energy 



of the qubit as 



u{5) = 2 — cos ^1 — cos ^2 + <^ [1 + cos {(^ — 5i — 82)] (33) 

where ix(5) = [/(^i, (52)/£^j, (/:) = 7r(2^-^o)/^o. The en- 
ergy is a periodic function of 5i and 82 with a period of 
2tt. At (^ ^ 1 and a > 1/2 each elementary cell contains 
two close minima, the barrier between them being small. 
In the eigen basis of the states representing these minima 
the qubit is then described by the effective Hamiltonan 
([1]) where the asymmetry he can be tuned by the embed- 
ded magnetic flux, ^. Shown in Fig. [12] is the plot u{5) 
for a = 0.75. The inset shows the potential profile along 
the line A for different (3 = (2$ - ^o)/*^o. To tune the 
tunneling parameter. A, one can replace the third junc- 
tion by two Josephson junctions connected in parallel, as 
it is shown in Fig. [IT] (inset). The Josephson energy of 
this circuit can be tuned by the magnetic flux through 
the second loop, which in turn can be tuned by an ex- 
ternal current-carrying line. Quantum superpositions of 
these states are obtained by pulsed microwave modula- 
tion of the enclosed magnetic flux by currents in control 
lines. Such a superposition has been demonstrated by 



Friedman et al (2000) and Van der Wal et al (2000). 



Though fabricated Josephson circuits exhibit a high 
level of static and dynamic charge noise due to charged 
impurities, while the magnetic background is much more 
clean and stable. The flux qubits can be driven individu- 
ally by magnetic microwave pulses; measurements can be 




qubit 
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FIG. 13 Cubic potential U sho wing qubit states and measure- 
ment scheme. Adapted from (Martinis et al. 2002). 



made with superconducting magnetometers (SQUIDs). 
They are decoupled from charges and electrical signals, 
and the known sources of decoherence allow for a de- 
coherence time of more than 1 ms. Entanglement is 
achieved by coupling the flux, which is generated by 
the persistent current, to a second qubit. The qubits 
are small (of order 1 mm) , they can be individually ad- 
dressed and integrated into large circuits. However, they 
are slower than the charge qubits. As it follows from 
Eq. (33), fluctuations of two parameters - magnetic flux. 



(0 



^, and Josephson energies Ey (or Josephson critical cur- 
rents J^ ) - are important for operation of flux qubits. 



Phase qubits, see, e. g., (Martinis et al. 2002), are 



designed around a 10 /im scale Josephson junction in 
which the charging energy is very small, thus provid- 
ing immunity to charge noise. Although still sensitive to 
flux, the circuit retains the quality of being tunable, and 
calculations indicate that decoherence from flux noise is 
small. The main part of a phase qubit is a current-biased 
Josephson junction, which can be characterized by the 
potential energy 



U{5) = -Ej[cos5^{I/Ij)5]. 



(34) 



Because the junction bias current / is typically driven 
close to the critical current Ij the tilted washboard po- 



tential ( 34 ) can be well approximated by a cubic potential 



with the barrier height 



A/7(/) = (4V2/3)^j [1 - I/Ij] 



3/2 



Therefore, the barrier can be tuned by the bias current; 
at / ^ /j it vanishes. The bound quantum states \n) 
with energy E^^ see Fig.[T3J can be observed spectroscop- 
ically by resonantly inducing transitions with microwaves 
at frequencies uornn = {Em — En)/h. The qubit state 
can be manipulated with dc and microwave pulses at fre- 
quency ujio of bias current. The measurement of the 
qubit state utilizes the escape from the cubic potential 
via tunneling. To measure the occupation probability pi 
of state |1), microwave pulses at frequency UJ21 driving 
a 1 ^ 2 transition were used. The large tunneling rate 
then causes state |2) to rapidly tunnel. Since the poten- 
tial profile depends on Ej, i. e., on the Josephson critical 
current, as well as on the bias current through the junc- 
tion, their fluctuations are important. Unlike the other 
qubits, the phase qubit does not have a degeneracy point. 
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Below we will briefly discuss main noise sources in flux 
and phase qubits. 



a. Flux noise - The origin of magnetic flux noise in 
SQUIDs with a power spectrum of the 1// type has been 
a puzzle for over 20 years. The noise magnitude, a few 
/i^o Hz^/^ at 1 Hz, scales slowly with the SQUID area 
and does not depend signiflcantly on the nature of the 
thin fllm superconductor or the substrate on which it 
is deposited. The substrate is typically silicon or sap- 
phire, which are insulators at low temperature (Well- 



stood et al. 



1987b). Flux noise of similar magnitude 



is observed in flux (Kakuyanagi et al., 2007 Yoshihara 



et a/.||2006[ ) and phase jBialczak et a/.[|2007[ ) qubits. The 



near-insensitivity to device area of the noise magnitude 
(normalized by the device area) ([Bialczak et al. 2007 



Lanting et~al\ |2009[ [Wellstood et al.\ |1987bD suggests 



that the origin of the noise is local. 



Koch et al. (2007) proposed a model in which electrons 



hop stochastically between traps with different preferen- 
tial spin orientations. They found that the major noise 
contribution arises from electrons above and below the 
superconducting loop of the SQUID or qubit, and that 
an areal density of about 5 x 10^^ cm~^ unpaired spins 
is required to account for the observed noise magnitude. 



de Sousa (2007) proposed that the noise arises from spin 



flips of paramagnetic dangling bonds at the Si-Si02 in- 



terface. Assuming an array of localized electrons, |Faoro 
et al. (2008) suggested that the noise results from elec- 
tron spin diffusion. The model was extended in (Faoro 



et al. 



2012) where it was shown that in a typical random 



configuration some fraction of spins form strongly cou- 
pled pairs behaving as two-level systems. Their switching 
dynamics is driven by the high-frequency noise from the 
surrounding spins, resulting in low- frequency 1// noise in 
the magnetic susceptibility and other physical quantities. 



Sendelbach et al. (2008) showed that thin-fllm SQUIDs 



are paramagnetic, with a Curie (ex T~^) susceptibil- 
ity. Assuming the paramagnetic moments arise from 
localized electrons, they de duced an areal den sity of 
"^. Subsequently, IB luhm et al. (2009) used a 



5 X 10^^ ' 

scanning SQUID microscope to measure the low-T para- 
magnetic response of (nonsuperconducting) Au rings de- 
posited on Si substrates, and reported an areal density 
of 4 X 10^^ cm~^ for localized electrons. Paramagnetism 
was not observed on the bare Si substrate. 



[Choi et al.\ ( |2009| ) have proposed that the local mag- 
netic moments originate in metal-induced gap states 
(MIGS) ( [Louie and Cohen 1976) localized by potential 
disorder at the metal-insulator interface. At an ideal in- 
terface, MIGS are states in the band gap that are evanes- 
cent in the insulator and extended in the metal, see 
Fig. 14 adapted from (Choi et al., 2009). At a metal- 



insulator interface there are inevitably random fluctu- 
ations in the electronic potential. The MIGS are par- 
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FIG. 14 (Color onine) (a) Schematic density of states, (b) 
MIGS at a perfect interface with energy in the band gap 
are extended in the metal and evanescent in the insulator. 
Adapted from ( |Choi et al.\ |2009| ) 



ticularly sensitive to these potential fluctuations, and a 
signiflcant fraction of them with single occupancy be- 
comes strongly localized near the interface, producing 
the observed paramagnetic spins. The local moments in- 
teract via mechanisms such as direct superexchange and 
the RKKY interaction between themselves, and Kondo 
exchange with the quasiparticles in the superconductor. 
This system, in principle, can exhibit a spin-glass transi- 
tion. However, experiments ( [Harris et a/. , 2008) suggest 
that at T > 55 mK the spins are in thermal equilibrium 
and exhibit a 1/T (Curie law) static susceptibility. To ex- 
plain the observed 1//" (0.6 < a < 1) noise spectrum of 
the magnetic flux the authors suggest that in this region 
one can use the fluctuation-dissipation theorem leading 
to the conclusion that fluctuations of the electronic mo- 
menta have also 1//" spectrum. Unfortunately, without 
knowing the form of the interaction between the spins, 
one cannot derive this behavior theoretically - this is still 
an open question. 



Recently Sank et al. (2012) measured the dependence 



of qubit phase coherence and flux noise on inductor loop 
geometry. They concluded that while wider inductor 
traces change neither the flux noise power spectrum nor 
the qubit dephasing time, increased inductance leads to a 
simultaneous increase in both. Another important result 
is the absence of scaling with the trace aspect ratio. 



b. Critical current noise - Noise of the Josephson crit- 
ical current in various superconducting qubits incorpo- 
rating Josephson junctions has been investigated in detail 



by Van Harlingen et al. (2004). They consider critical- 



current fluctuations caused by charge trapping at defect 
sites in the tunneling barrier and compare their contri- 
bution to the dephasing time with that of the flux noise 
due to hopping of the vortices through the SQUID loop. 



see Fig. 15 adapted from (Van Harlingen et al. 2004). 



This mechanism can usually be made negligible in de- 
vices fabricated with linewidths less than approximately 
^y^o/B for which vortex trapping in the line is sup- 
pressed ( Dantsker et al.\ |T996); here B is the fleld in 
which the device is cooled. The trapped charges block 
tunneling through a region of the junction due to the 
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FIG. 15 Flux modulation from vortices hopping into and 
out of a loop, and critical-current modulation from electrons 
temporarily trapped at defect sites in the junction barrier. 



Adapted from ( Van Harlingen et al. 2004 ) 



Coulomb repulsion, effectively modulating the junction 
area. In general, a single-charge fluctuator produces a 
two-level, telegraph noise in the critical current of a junc- 
tion characterized by the lifetimes r^ in the untrapped 
state (high critical-current), and Tt in the trapped state 
(low critical-current). This produces a Lorentzian peak 
in the power spectral density with a characteristic rate 
7 = r~^ +r^~^. Experiments on dynamics of such fluctu- 
ators and their life times (Rogers and Buhrman, 1984 



1985 Wakai and Harlingen, 1986) provide strong evi 



dence that the dominant charges enter the barrier from 
one electrode and exit to the other (voltage dependence), 
and that the fiuctuators exhibit a crossover from thermal 
activation to tunneling behavior at about 15 K. In the 
tunneling regime, the fluctuating entity has been shown 
to involve an atomic mass, suggesting that ionic reconfig- 
uration plays an important role in the tunneling process, 
see also discussion in (Galperin et al., 1989). Though 



interactions between traps resulting in multiple-level hi- 



erarchical kinetics have been observed (Wakai and Har- 



lingen 1987), usually the traps can be considered to be 
local and noninter acting. In this limit, the coexisting 
traps produce a distribution of Lorentzian features that 
superimpose to give a l//-like spectrum. Careful analy- 
sis of the influence of the noise in the Josephson critical 
current on different qubit designs and various data ac- 



quisition schemes, has led Van Harlingen et al (2004) 



to the conclusion that although there is strong evidence 
that the noise derives from a superposition of random 
telegraph signals produced by charge trapping and un- 
trapping processes, the origin of 1// noise in the criti- 
cal current of Josephson junctions is still not fully un- 
derstood. In particular, the origin of the oc T^ depen- 
dence of the noise power observed by |Wellstood et al\ 
( 1987a|b , 2004) remains puzzling. To account for this be- 



havior IShnirman et al. ( 2005 ) suggested that the density 



of states for two-level fiuctuators is proportional to their 
inter- level spacing, E. This assumption was supported 
by a microscopic model by Faoro et a/.|(|2005 ). However, 
detailed measurements by Eroms et al.\ (|2QQ6) showed 



utilized in many superconductor-based qubits, the noise 
power is (X T between 150 and 1 K rather than ex T^. 
The observed spectral density saturates below 0.8 K due 
to individual strong two-level fiuctuators. Interestingly, 
the noise spectral density at 4.2 K is two orders of mag- 
nitude lower than expected from the literature survey of 



(Van Harlingen et al, 2004). 



The role of the critical current noise is convention- 
ally allowed for by introducing a coupling term into the 
Hamiltonian 



Ku and Yu 



Faoro et al. 



2005 



Faoro and loffe 



2006 



20051 IShnirman et al.\ |2005| [Simmonds 



et al.\ 2004). A simple microscopic model relevant to 



Al/AlO^^/Al was developed by |Constantin and Yu] ( |2007 ). 
This model leads to the scaling of the 1// noise with the 
junction thickness as ex L^. The results are in a reason- 
able agreement with corresponding experimental values 



of (Eroms et al. 2006 Zimmerli et al., 1992 Zorin et al. 



1996). However, to the best of our knowledge, the pre- 



dicted scaling has not yet been verified. 



c. Decoherence in Josephson qubits from dielectric losses - 
As it was mentioned in Sec. |II.A[ amorphous ma- 
terials contain low-energy excitations behaving as two- 
level tunneling systems. These states are responsible for 
low-temperature thermal and kinetic properties of struc- 
tural glasses, in particular, for specific heat and dielec- 
tric losses. Due to interaction with their environment, 
the TLSs switch between their states producing low- 
frequency noise. These noises act on a qubit reducing 
its coherence time. 



Martinis et al. ( 2005 ) pointed out that the noises pro- 



duced by TLSs in amorphous parts of qubit devices are of 
primary importance. The reason is that crossover wiring 
in complex superconducting devices requires an insulat- 
ing spacer that is typically made from amorphous Si02 
deposited by chemical vapor deposition (CVD). They 
performed a variety of microwave qubit measurements 
and showed that the results are well modeled by loss from 
resonant absorption of two- level defects. Dielectric loss 
(loss tangent) in a system formed by a superconductor 
lead and a 300 nm thick CVD Si02 layer was measured 
at / - 6 GHz and T = 25 mK <C hw/kB^ Generally, 
two mechanisms contribute to the dissipation induced by 
the two- level defects. The first - resonant - is due to 
direct microwave- induced transitions between the TLS's 
level with subsequent emission of phonons. The second 
one is due to the relaxation-induced lag in phase between 
the nonequilibrium level populations and the driving ac 
electric field. One can expect that at hw ^ ksT the 
first mechanism should dominate, see, e. g., ( jHunklinger 
and von Schickfus 1981). A hallmark of the resonant 



that in shadow-evaporated Al/AlO^^/Al tunnel junctions 



absorption is its strong dependence on the amplitude of 
the applied ac electric field. This dependence is due to 
decrease of the difference between the occupancies of the 
upper and lower level with the field amplitude increase. 
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The theoretical prediction for the loss tangent f3 is, see, 
e. g., ( von Schickfus and Hunklinger[ |1977| ) 



insensitive both to the charge and flux noise was imple- 



(3 



7rP(ed)2 tanh(hLu/2kBT) 



(35) 



Here P is the TLS density of states each having a fluctu- 
ating dipole moment ed and relaxation times Ti and T2, 
X is the dielectric constant, while ujr = eSd/h is the TLS 
Rabi frequency corresponding to the ac fleld amplitude, 
£. This expression is derived under the assumption that 
distributions of inter- level spacings and logarithms of the 
relaxation times of the TLSs are smooth. The theory flts 
the experimental data well with parameters compatible 



with previous measurements of bulk Si02 (von Schickfus 



and Hunklinger 1977). The above results are consis- 



tent with previous measurements of an AlO^^ capacitor 



by |Chiorescu et a/.| ( [2QQ4| . 

A key difference between tunnel junctions and bulk 
materials is that tunnel junctions have small volume, and 
the assumption of a continuous distribution of defects 
is incorrect. Instead, dielectric loss must be described 
by a sparse bath of discrete defects. Indeed, individual 
defects were measured spectroscopically with the phase 



qubit ([Cooper et a/T| [2004 [ [Lisenf eld et"aH|2QlQa|b[|Sim- 



monds et al. , 2004 ) . They are observed as avoided cross- 
ings in the plots of the occupation probability versus 
qubit bias. A qualitative trend is that small- area qubits 
show fewer splittings than do large-area qubits, although 
larger splittings are observed in the smaller junctions. It 
follows from quantitative analysis of the number of res- 
onances that couple to the qubit that at large area the 
decoherence rate is compatible with the loss tangent of a 
bulk material. 



Based on the obtained results. Martinis et al. (2005) 



have formulated the following trends for making devices 
with long coherence times (i) usage small-area junctions 
where number of two-level defects is small; (ii) usage of 
simple designs with no lossy dielectrics directly connected 
to the qubit junction; (iii) trying to flnd insulating ma- 
terials with low dielectric losses. 

A way to eliminate the effects of low-frequency charge 
(and flux) noise is to operate the system close to the 
degeneracy point [e = in the Hamiltonian (fl])]. In- 
deed, since the distance between the energy levels is 
hVt = h\/e^ + A^, the fluctuation part hSe{t) of the po- 
tential energy he{t) creates the additional shift 

1A2 

n^^'^' ' 2 0^ 

Therefore, at e = the charge noise vanishes in the lin- 
ear approximation and only second-order contributions 
are important. For a charge qubit, the degeneracy point 
corresponds to the induced charge q = e while for a flux 
qubit this point corresponds to an integer number of the 
half-flux quanta in the loop. This idea of tuning the de- 
vice to the double degeneracy point where the qubit is 



Sn{t) = - Se{t) + 



Mt)? 



(36) 



mented in the quantronium device (Cottet et al. 2002 



Ithier et al. 2005 Vion et al. 2002). 
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FIG. 16 (Color online) Idealized circuit diagram of the 
quantronium, a quantum-coherent circuit with its tuning, 
preparation, and readout blocks. The circuit consists of a 
CPB island (black node) delimited by two small Josephson 
junctions (crossed boxes) in a superconducting loop. The loop 
also includes a third, much larger Josephson junction shunted 
by a capacitance C. The Josephson energies of the box and 
the large junction are Ej and Ejq. The Cooper pair number 
A^ and the phases 6 and 7 are the degrees of freedom of the 
circuit. A dc voltage U applied to the gate capacitance Cg 
and a dc current /$ applied to a coil producing a flux $ in 
the circuit loop tune the quantum energy levels. Microwave 
pulses u{t) applied to the gate prepare arbitrary quantum 
states of the circuit. The states are read out by applying a 
current pulse Ib{t) to the large junction and by mo nitoring the 
voltage V{t) across it. From (Vion et al. 2002J). Reprinted 
with permission from AAAS. 

The main principle behind the device is shown in 
Fig. [16] reprinted from (|Vion et al] [2002]). The CPB 



involves two Josephson junctions with a capacitance Cg 
connected to the island separating them. The junc- 
tions are connected to a third, larger junction, with 
a larger Josephson energy, to form a superconducting 
loop threaded by a magnetic flux <I>. To achieve in- 
sensitivity to the charge noise the qubit is operated at 
Ng = CgU/{2e) =0.5, where the energy levels have zero 
slope and the energy-level splitting is Ej. Insensitivity 
to the flux noise is achieved by applying an integer num- 
ber of half-flux quanta to the loop. Therefore, a double 
degeneracy point can be achieved by tuning dc gate volt- 
age U and the current /«^. To measure the qubit state 
one has to shift the qubit from this point. This was 
achieved by the current pulse Ib{t) applied to the loop, 
which produces a clockwise or counterclockwise current 
in the loop, depending on the state of the qubit. The 
direction of the current is determined by the third (read- 
out) junction since the circulating current either adds or 
subtracts from the applied current pulse. As a result, the 
read-out junction switches out of the zero- volt age state 
at different values of the bias current. Thus, the state 
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of the qubit was determined by measuring the switching 
currents. In the quantronium, much longer relaxation 
and decoherence times can be achieved compared with 
conventional CPB. 



D. Semiconductor-based qubits 

Here we will briefly discuss main trends in designing 
the qubits based on semiconductors and semiconductor 
devices. One can find more detailed discussion in the 



1 ^im 



reviews ( Chirolli and Burkard 2008 Hanson et al. 2007 



Liu et a/., 2010 Zak et a/., 2010) 



1. Spin qubits 

It is natural to choose the electron spin as the two- 
state system that encodes the qubit. In modern semi- 
conductor structures the spin of the electron can have a 
much longer coherence time than the charge degrees of 
freedom. However, it is not easy to isolate, control, and 
manipulate the spin degree of freedom of an electron to a 
degree required for quantum computation. A successful 
and promising device for the physical implementation of 
electron spin-based qubits is the semiconductor quantum 



dot (Loss and DiVincenzo 1998). 



The quantum dots (QDs) considered for implementa- 
tion of quantum algorithms are confined regions of semi- 
conductor materials coupled with reservoirs by tunable 
tunnel barriers. The height of the barriers, and con- 
sequently the rates for tunneling through the barriers 
on and off the dot, can be controlled via the applica- 
tion of gate electrodes. The dots are actually quantum 
boxes having discrete energy levels, their positions with 
respect to the chemical potential of the reservoir can be 
also tuned by electrostatic potentials. Therefore, QDs 
can be considered as tunable artificial atoms. Coulomb 
interaction between the electrons (or holes) occupying 
the dot's levels determine the energy cost for adding an 
extra electron. Because of this cost, the electron trans- 
port through the dot can be strongly suppressed at low 
temperatures (the so-called Coulomb blockade). Since 
the energy cost can be tuned by gate electrodes, the QD 
devices are promising for many applications. 

Among many types of quantum dots, devices based 
on lateral III-V semiconductor QDs are of special inter- 
est. Such devices are usually fabricated from heterostruc- 
tures of GaAs and AlGaAs grown by molecular beam epi- 
taxy. In such heterostructures the electron motion can 
be confined to a thin layer along the interface forming a 
two-dimensional electron gas (2DEG) with high mobility 
(- 10^ - 10^ cmVV- s) and low density (- 10^^ cm-^). 
The low density results in a relatively long Fermi wave- 
length (~40 nm) and a large screening length. Therefore, 
the 2DEG can be locally depleted by an electrostatic field 




FIG. 17 SEM image of a double quantum dot device . A 
quantum point contact with conduct ance, gs^ senses charge 
on the double dot. Reprinted from (|Petta et al. 2008). 



applied to a metal gate electrode allowing to design quan- 



tum devices similar to that shown in Fig. 17 When the 
lateral size of the dot is comparable with the Fermi wave- 
length the distance between the discretized energy levels 
becomes larger than the temperature (at temperatures of 
tens millikelvin) , and quantum phenomena become im- 
portant. 



Spin relaxation and dephasing in quantum dots - Two 
kinds of environment turn out to mainly affect the dy- 
namics of an electron spin in a quantum dot, the phonons 
in the lattice, and the spins of atomic nuclei in the quan- 
tum dot. 



Starting from (Khaetskii and Nazarov 2000, 2001 ), the 
phonon-induced relaxation was extensively studied. It 
turns out that the lattice phonons do not couple directly 
to the spin degree of freedom. However, even without 
the application of external electric fields, the breaking 
of inversion symmetry in GaAs gives rise to spin-orbit 
(SO) interaction, which couples the spin and the orbital 
degrees of freedom. These orbital degrees of freedom, 
being coupled to the phonons, provide an indirect cou- 
pling between the electron spin and the phonons, which 
constitute a large dissipative bosonic reservoir and pro- 
vide a source of decoherence and relaxation. Short time 
correlations in the phonon bath induce a Markovian dy- 
namics of the electron spin, with well-defined relaxation 
and decoherence times Ti and T2. 



As we discussed in Sec. LB in the Bloch picture, pure 
dephasing arises from longitudinal fluctuations of the 
magnetic field, while a perturbative treatment of the SO 
interaction gives rise, within first order, to a fluctuat- 
ing magnetic field perpendicular to the applied magnetic 
field. As a consequence the decoherence time T2 is lim- 
ited only by its upper bound Ti, T2 = 2Ti. 

Hyperfine interaction was first taken into consideration 
as a source of decoherence for an electron spin confined in 



a quantum dot by Burkard et al (1999). This interaction 



is very important since in a ^40 nm GaAs quantum dot 
the wavefunction of an electron overlaps with approxi- 
mately 105 nuclei. The electron spin and the nuclear 
spins in the dot couple via the Fermi contact hyperfine 
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FIG. 18 Electric circuit model of a DQD containing m and n 
electrons in the left and right dot, respectively. The two QDs 
are coupled to each other via a tunnel junction, and each dot 
is connected to an electron reservoir via a tunnel junction. 
The electrostatic potential of the left (right) quantum dot is 
controlled by the gate vo ltage Vi {Vr ) through t he capacitance 
Ci (Cr). Adapted from (Fujisawa et al. 2006). 




FIG. 19 Fragment fro m the stability diagram. Adapted from 
( Fujisawa et al.\ 2006j ). 



interaction, which creates entanglement between them 
and strongly influences the electron spin dynamics. It 
turns out that long-time correlations in the nuclear spin 
system induce a non-Markovian dynamics of the electron 
spin, with non-exponential decay in time of the expecta- 
tion values of the electron spin components. 

Relative importance of the above mechanisms of de- 
coherence depends on the external magnetic field: the 
phonon-induced relaxation rate of the electron spin is 
enhanced by an applied magnetic field, whereas the in- 
fluence of the hyperfine interaction is reduced by a large 
Zeeman splitting. 

In this review, we will not focus on decoherence in 
semiconductor spin qubits, which has been extensively 
reviewed, e. g., by Chirolli and Burkard ( 2QQ8| ). More re- 
cent papers related to random telegraph or 1// noise ad- 
dress charge traps near the interface of a Si heterostruc- 
ture (Culcer et a/., 2009), as well as nearby two-level 



charge fluctuators in a double dot spin qubit (Ramon 



andHu||2009J . 



2. Charge qubits 

In semiconductor systems, a charge qubit can be 
formed by isolating an electron in a tunnel-coupled dou- 
ble quantum dot (DQD) ([Fujisawa et 



2006 Hayashi 



et al 2003). Here we will discuss a DQD consisting of 



two lateral QDs, which are coupled to each other through 
a tunnel barrier. Each QD is also connected to an elec- 
tron reservoir via a tunnel junction. The operation of 
such device can be analyzed from its electric circuit model 
shown in Fig. [ISJ The transport properties of a semicon- 
ductor DQD have been studied extensively, see reviews 



(Grabert and Devoret, 1991 Kouwenhoven et a/., 1997 



van der Wiel et al. , 2003 ) . Each tunnel barrier has a 



small coupling capacitance, Ci, as well as a finite tunnel- 
ing coupling, Tc, and single-electron transport through 
the DQD can be measured. The tunneling rates for the 
left and right tunnel barriers are denoted as Tl and F^^, 



respectively. In addition, the DQD is connected to gate 
voltages Vi and Vr via capacitors Ci and Cr, respectively, 
so that the local electrostatic potential of each dot can be 
controlled independently. The energy difference between 
the electrochemical potentials fiL and fiR of the left and 
right reservoirs corresponds to the applied source - drain 
voltage Vsd- The stable charge configuration (m, n), with 
m electrons on the left QD and n electrons on the right 
one, minimizes the total energy in all capacitors minus 
the work that has been done by the voltage sources. 
This energy can be estimated from the equivalent cir- 
cuit Fig. flS) see, e. g., analysis in the book by jHeinzeT 



(2007), Sec. 9.3. A fragment of the stability diagram in 
the Vi — Vr plane is shown in Fig. [19] When the tun- 



neling coupling, Tc, is negligibly small, the boundaries 
of the stable charge states appear as a honeycomb pat- 
tern, a part of which is shown by dashed lines. The triple 
points, E and H^ of three charge states are separated by 
a length corresponding to the inter-dot Coulomb energy 
U. Electrons pass through three tunnel barriers sequen- 
tially in the vicinity of triple points. On the other hand, 
the tunneling process at H can be viewed as hole trans- 
port, as the unoccupied state (hole) moves from the right 
to the left (not shown in the diagram). 

When the tunneling coupling is significantly large, the 
charging diagram deviates from the honeycomb pattern 
as shown by solid lines in Fig. 19 Due to quantum repul- 
sion of the levels, the minimum distance between A and 
B is increased by the coupling energy HA = hTc from its 
original value U. 

The idealized dynamics of the qubit can be described 
assuming that each dot has a single energy level, e^. Then 
the effective Hamiltonian can be expressed in the form 
(IT]) with he = cl — cr and A = Tg. The parameters of 
the Hamiltonian can be tuned by the gate voltages. 

Coherent control of a GaAs charge qubit has 
been demonstrated (Fujisawa et al 2004; Hayashi 



et 



2003), along with correlated two-qubit interac- 
tions ([Shinkai et al. 2009). In these experiments each 



dot contained a few tens of electrons, potentially com- 
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FIG. 20 Scanning electron microphotograp h of a device sim- 
ilar to the one measured. Adapted from ( [Petersson et al] 
[2Q1QJ . 



plicating the qubit level structure. In addition, the dots 
were strongly coupled with the leads typically limiting 
coherence times to ~ 1 ns (due to quantum co-tunneling) . 
A coherent control of tunable GaAs charge qubit con- 
taining a single electron was demonstrated by IPetersson 



et aL\ (2Q1Q[), see Fig. 20 The gate electrodes are ar- 



ranged in a triple QD geometry and deplete 2DEG in the 



GaAs/AlGaAs heterostructure (Petta et al 2010). The 



DQD was formed using the left and middle dots of the 
structure while the right side of the device was configured 
as a non-invasive quantum point contact (QPC) charge 
detector. The device operated near the (1,0) — (0,1) 
charge transition, where {riL^nR) denote the absolute 
number of electrons in the left and right dots. The level 
detuning, he, was adjusted by the voltage Vr on the gate 
R while the tunnel splitting HA was adjsuted by the volt- 
age Vm on the gate M. The coherence time was extracted 
as a function of the detuning (from the charge degeneracy 
point) voltage, the maximal value being ~ 7 ns. The re- 
sult is ascribed to 1// noise, whose influence is analyzed 
along the Gaussian assumption. 

In spite of the successful manipulation of a single- 
charge qubit, the qubit is actually influenced by uncon- 
trolled decoherence, which is present even in the Coulomb 
blockade regime. Several possible decoherence mecha- 
nisms were discussed. 

First, background charge (1//) noise in the sample and 
electrical noise in the gate voltages cause fluctuation of 
the qubit parameters e and A, which gives rise to deco- 



herence of the system ( Hayashi et al] |2003[ [Itakura and| 



Tokura 2003 Paladino et al 2002). The amplitude of 



low-frequency fluctuation in he is estimated to be about 
1.6 /ieV, which is obtained from low- frequency noise in 
the single-electron current, or 3 /ieV, which is estimated 
from the minimum line width of an elastic current peak 



at the weak coupling limit (Fujisawa et a/., 1998). Low- 



frequency fluctuations in HA are relatively small and es- 
timated to be about 0.1 /ieV for HA = 10 /ieV, assuming 



ence rate observed at the off-resonant condition (e > A). 
1// is usually attributed to a set of bistable fluctuators, 
each of which produced Lorentzian spectrum. However, 
the microscopic origin of the charge fluctuators is not 
fully understood, and their magnitude differs from sam- 
ple to sample, even when samples are fabricated in the 



same batch. It was shown ( Jung et al , 2004 ) that fluc- 



tuations in e can be reduced by decreasing the temper- 
ature as suggested by a simple phenomenological model 
where the activation energy of the traps is uniformly dis- 
tributed in the energy range of interest. Cooling samples 
very slowly with positive gate voltage is sometimes ef- 
fective in reducing charge fluctuation at low temperature 
( Pioro-Ladriere et al 2005|. Authors suggest that the 
noise originates from a leakage current of electrons that 
tunnel through the Schottky barrier under the gate into 
the conduction band and become trapped near the ac- 
tive region of the device. According to [Buizert et al 
( 2008| ), an insulated electrostatic gate can strongly sup- 
press ubiquitous background charge noise in Schottky- 
gated GaAs/AlGaAs devices. This effect is explained by 
reduced leakage of electrons from the Schottky gates into 
the semiconductor through the Schottky barrier. A sim- 
ilar result has been recently reported in Schottky gate- 
deflned QPCs an DQDs in Si/SiGe heterostructures with 
a global top gate voltage ( [Takeda et al. 2013). By neg- 



atively biasing the top gate, 1//^ switching noise is sup- 
pressed in a homogeneous 1// charge noise background. 
It is suggested that this technique may be useful to elim- 
inate dephasing of qubits due to charge noise via the 
exchange interaction (Culcer et al. 20091). 



In contrast, the decoherence at the resonant condi- 
tion (e = 0) is dominated by other mechanisms. Al- 
though the flrst-order tunneling processes are forbidden 
in the Coulomb blockade regime, higher-order tunneling, 
namely co-tunneling, processes can take place and deco- 
here the system (lEtoj 120011) . Actually, the co-tunnelling 
rate estimated from the tunneling rates is close to the 
observed decoherence rate and may thus be a dominant 



mechanism in the present experiment (Eto, 2001). How- 



ever, since one can reduce the co-tunneling effect by mak- 
ing the tunneling barrier less transparent, it is possible 
to eventually eliminate it. 

The electron-phonon interaction is an intrinsic deco- 
herence mechanism in semiconductor QDs. Spontaneous 
emission of an acoustic phonon persists even at zero tem- 
perature and causes an inelastic transition between the 



two states (Fujisawa et a/., 1998). The phonon emission 
rate at the resonant condition e = cannot be directly 
estimated from the experimental data on the FID-type 
protocols, but it may be comparable to the observed de- 
coherence rate. Strong electron-phonon coupling is re- 
lated to the fact that the corresponding phonon wave- 



length is comparable to the size of the QD (Brandes 



local potential fluctuations in the device (Jung et al. , and Kramer 1999 Fujisawa et al. 1998). In this sense. 



2004). Actually, the e fluctuation explains the decoher- electron-phonon coupling may be reduced by optimizing 
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the size of QD structures. In addition, polar semicon- 
ductors, such as GaAs, exhibit a piezoelectric type of 
elect ron-phonon coupling, which is significant for low- 
energy excitations (< 0.1 meV for GaAs). Non-polar 
semiconductors, such as Si or carbon-based molecules, 
may be preferable for reducing the phonon contribution 
to the decoherence. 



III. DECOHERENCE DUE TO 1// NOISE 

During the last decade we witnessed of an extraordi- 
nary progress in quantum devices engineering reaching a 
high level of isolation from the local electromagnetic envi- 
ronment. Under these conditions, the material-inherent 
sources of noise play a crucial role. While the microscopic 
noise sources may have different physical origin, as elu- 
cidated in Section [nj their noise spectral densities show 
similar l//-type behavior at low frequencies. Material- 
inherent fluctuations with 1// spectrum represent the 
main limiting factor to quantum coherent behavior of the 
present generation of nanodevices. This fact has stimu- 
lated a great effort in understanding and predicting de- 
coherence due to 1// noise and to the closely related RT 
noise. 

There are two characteristic features which make any 
prediction of decoherence due to 1// noise quite compli- 
cated. Firstly, stochastic processes with 1// spectrum 
are long-time correlated. The spectral density of the 
noise increases with decreasing frequency down to the 
lowest experimentally accessible frequencies. The mea- 
surement frequency band is limited either by frequency 
filters or simply by the finite duration of each realization 
of the random process, set by the measurement time, t^. 
In particular, during tm some of the excitations respon- 
sible for the noise may not reach the equilibrium. For 
this reason 1// noise is considered as a typical nonequi- 
librium phenomenon for whic h fluctuat i on-dis sipation re- 
lations may not hold ([Galperin et at. 2QQ4|)F] The ex- 



istence of relaxation times longer than any finite mea- 
surement time tm also poses the question of stationarity 



^ In some magnetic systems, for instance in spin glasses, exper- 
iments evidenced magnetic noise with 1// spectrum satisfying, 
within experimental accuracies, fluctuation-dissipation relations. 
Despite being non-ergodic, in these systems magnetic noise is 
with satisfactory accuracy a thermal and equilibrium one. The 
mechanism of magnetic fluctuations is, however, not yet clear. 
Combined studies of magnetic and resistance noise, of response 
to ac magnetic fleld and decay of magnetization in spin glasses 
and amorphous ferromagnets revealed the role of disorder and, 
speciflcally, the role of the spectrum of metastable states, sepa- 
rated by barriers, which are created by disorder, in the formation 
of the noise spectrum. The problem of kinetics of spin glasses, 
determined both by many-body competing interactions and dis- 
order, is extremely complicated and beyond the scope of this 
review. In connection with 1// noise it is addressed in ( |Kogan[ 
19961. 



of l//-type noise, a problem which has attracted much 
attention (|Kogan , 1996). For a stationary process re- 
peated measurements yield the same power spectrum 
Sx{(jo) within experimental accuracy. In some systems 
variations of the spectrum or its "wandering" have been 



reported, see (Kogan 1996 Weissman, 1993). These ob- 



servations are more reasonably attributed to incomplete 
ergodicity of the fluctuating system rather than to non- 
stationarity. In those cases the long relaxation times may 
correspond to rare transitions of the system, overcoming 
high energy barriers, from one "valley" in phase space to 
another one in which the spectrum of relaxation times 
r < tm is different. Several experiments tried to re- 



veal a possible nonstationarity of 1// noise (see (Kogan 



1996) and references therein), but no clear manifestation 
of nonstationarity has been found. 

Secondly, in general 1// noise cannot be assumed to be 
a Gaussian random process. Even if the probability den- 
sity function of many 1// processes resembles a Gaussian 
form (which is necessary, but not sufficient to guarantee 
Gaussianity), clear evidences demonstrating deviations 
from a Gaussian random behavior have been reported, 
see for a review (Weissman, 1988), Sec. IIIA, and (iKo-1 



gan 



1996), Sec. 8.2.2. The explanation of the variety of 



observed behaviors stems from the fact that the mecha- 
nisms of 1// noise may be different in various physical 
systems, implying that their statistical properties may 
strongly differ. 

As a consequence, the two standard approximations 
allowing making simple predictions on the evolution of 
open quantum systems, the Markovian approximation 
and the modelization of the environment as a bath of 
harmonic oscillators, cannot be straightforwardly applied 
when the power spectrum of the noise is of 1// typejj In 
the context of quantum computation, the implication of 
long-time correlations of the stochastic processes is that 
the effects of 1// noise on the system evolution depend 
on the specific quantum operation performed and/or on 
the measurement protocol. Some protocols clearly reveal 
the non-Gaussian nature of the process, whereas for oth- 
ers a Gaussian approximation captures the main effects 



at least on a short time scale (Falci et al. 2005 Makhlin 



and Shnirman, 2004 Rabenstein et al. 2004) 



^ We remind that the Markovian approximation for the reduced 
dynamics of a quantum system is applicable provided that the 
noise correlation time, Tc-, and the syste m-bath coupling strength , 
-u, satisfy the condition vtc <C 1 (Cohen- Tannoudji et a/.||l992| ). 
In the case of 1// noise Tc -^ oo, excluding the possibility to 
apply any Markovian approximation for the reduced system evo- 
lution. We observe that this is not in contrast with the fact that 
1// noise can originate from a superposition of RT processes, 
each of them (even the slowest) being a discrete Markov process. 
The above inequality in fact guarantees that the quantum system 
does not perturb the environment and, therefore, no memory of 
its previous states is left. 
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In this Section we will illustrate various approaches de- 
veloped in recent years to address the problem of decoher- 
ence due to noise sources having 1// spectrum, consider- 
ing both microscopic quantum models and semi-classical 
theories. We will start addressing the decohence prob- 
lem in single-qubit gates driven by dc-pluses or by ac- 
fields, then we will consider more complex architectures 
needed to implement the set of universal gates. To cast 
these problems in a general framework, we introduce here 
the general Hamiltonian of a nanodevice including a phe- 
nomenological description of the environment. In some 
cases this general structure results from a microscopic 
description of the device, including the most relevant en- 
vironmental degrees of freedom, as illustrated in the fol- 
lowing. In general, the Hamiltonian for a nanodevice can 
be cast in the form 



^tot = ^0 



Hc(t)^Hn(t)^HR^Hi. 



(37) 



where i^o + Hcif) describes the driven closed system, 
classical noise affecting the system is included in Hn(f) 
and Hr + Hi represents quantum environment and its 
interaction with the system. This general Hamiltonian 
can be justified by the following considerations. The 
macroscopic Hamiltonian of the device ^o[q] is an op- 
erator onto a A/'-dimensional Hilbert space, H. It de- 
pends on the parameters q used to fix the bias (op- 
erating) point and allowing tunability of the device. 
The eigenstates of i/'o[q], {|^i(q)) : i = 1,...,A^}, 
form the "local basis" of the "laboratory frame" where 
^0 = Z)i=i^i(q) l</^i(q)) (</>i(q)l- External control 
is described by a time-dependent term. In a one-port 
design, the driving field A{t) couples to a single time- 
independent system operator Q, 



H,{t)=A{t)Q, 



(38) 



which is Hermitian and traceless. In general, control is 
operated via the same ports, i. e., by allowing the bias 
in Hq[c\\ to depend on time. This suggests to split q in 
a slow part, q(t), which includes static bias, and the fast 
control parameter, qc(t), as q ^ q(t) + qc(^)- Accord- 
ingly, we split 



HoHt) + qe(t)] := ifo[q(t)] + H,{t) 



(39) 



where Hc(f) describes (fast) control; in relevant situa- 
tions it can be linearized in qc(^), yielding the structure 
of Eq. (38) J The interaction with the complicated en- 
vironment of microscopic degrees of freedom in the solid 



^ Notice that while Q is an observable and does not depend on the 
local basis, its matrix representation does. The physical conse- 
quence is that the effectiveness of the fast control qc(^) in trig- 
gering transitions depends also on the slow q(t), a feature of 
artificial atoms reflecting their ease of tunability. 



state can be modeled by a phenomenological Hamilto- 
nian. Usually classical noise also acts through the same 
control ports, thus it can be modeled by adding a stochas- 
tic component (^q(t) to the drive. Again, we split slow 
and fast noise, (^q(t) -^ ^q(^) + ^.f{^)i ^^^ include the 
slow part in Hq[g^ -^ Ho[<l{t) + ^q(^)]- The same steps 
leading to Eq. (39) yield the noisy Hamiltonian 



fftot = HoHt) + 6q{t)] + He{t) + Hf 



(40) 



where Hf describes Markovian classical noise. "Quanti- 
zation" of this term, Hf{t) -^ E Q ^ Hr = Hj -\- Hr, 
yields the final system-environment Hamiltonian in the 
form of Eq. (|37|, where ^o[q(^)] + ^n(^) = Ho[(i{t) + 
(5q(t)], E operates on the environment and Hr is its 
Hamiltonian, plus possibly suitable counter-terms. 

When the nanodevice operates as a qubit, ^o[q] can 
be cast in the form Eq. ^ , where both the level splitting 
Q and the polar angle shown in Fig. [l] depend on the 
parameters q. For simplicity, we suppose that a single 
parameter is used to fix the bias point, q. We can write 



Eq. (40) as 



Hu 



h 



n[q + Sq{t)] • a + Hc{t) -^Hj^Hr. (41) 



Expanding ^[q-\- Sq{t)] about the fixed bias g, we obtain 

^tot = ^n{q)-a^Sq{t)^ •a^H,{t)^Hi^HR , (42) 

and if bias q controls only one qubit component, for in- 
stance ^z{q)^ then Eq. (42) reduces to the commonly 
used form 



Hint = 






(^x 



n,{q)^E{t)^E 



CTz 



-H,{t)- 



-Hr . 

(43) 

Here the qubit working point is parametrized by the angle 
9q which is tunable via the bias g, tan^^ = Vtx/^z{(l)^ 
the classical noise term is E(t) = 2Sq{t){dQz/dQ)^ and 
we have consistently set Q = cr^ in Hj. 

From the physical point of view, the phenomenologi- 



cal Hamiltonian (43) treats on different footings fast en- 



vironmental modes exchanging energy with the system 
and and slow modes responsible for dephasing. The fast 
modes are treated quantum mechanically and included in 
hEaz/2-\-HR while the slow ones are accounted for clas- 
sically and included in Hn{t) = hE{t)(7z/2. This term 
corresponds to the longitudinal fiuctuating part, b(t), of 
a "magnetic field" formed by a controlled part Bq plus a 
classical stochastic component, as introduced in Eq. (Isl). 
Results of measurements involve both quantum and clas- 
sical ensemble averaging. From the technical point of 
view, effects of Markovian noise, included in Hi -h Hr^ 
are studied by weak coupling quantum optical Master 



Equations (Bloch 1957 Cohen-Tannoudji et at., 1992 



Redfield 1957 Weiss, 2008), which lead to exponential 
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decay both of the diagonal (populations) and of the off- 
diagonal (coherences) elements of qubit density matrix 
elements in the Hq eigenbasis. The corresponding time 
scales, Ti and T2, are given by Eqs. (15) and (16), re- 
spectively. 

This approach fails for the classical stochastic term 
Hn{t) which describes low- frequency (e. g., 1//) noise. 
An important feature of superconducting nanodevices is 
that the Hamiltonian Ho[q\ can be tuned in a way such 
that parity symmetries are enforced where dEi/dq = 
and selection rules hold affecting the matrix elements Qij 
in the local basis. In these parity symmetry points the 
device is well protected against low frequency noise, and 



the system is said to operate at "optimal points" (Vion 



et al. 2002). In general, solid-state devices suffer from 



broadband and non-monotonic noise formally described 
by Hnit) and Hi + Hr in Eq. ^. Therefore, ap- 



proaches suitable to deal with noise acting on very differ- 
ent time scales are required. This topic will be discussed 
in Sec. IIII.B.l.bl 

To start with, we analyze the effect of 1 //-noise on the 
qubit's evolution under pure dephasing conditions. As we 
discussed in Sec. |I.B[ in some situations the Gaussian ap- 
proximation does not apply. In these cases knowledge of 
only the noise power spectrum S{uj) is not sufficient since 
noise sources with identical power spectra can have differ- 
ent decohering effects on the qubit. Therefore, it is nec- 
essary to specify the model for the noise source in more 
detail. In Sec. |III.A we illustrate dephasing by the spin- 
fluctuator model, which can be solved exactly under gen- 
eral conditions. When we departure from pure dephas- 
ing, i. e., when the noise-system interaction is not lon- 
gitudinal, no exact analytic solution for the system time 
evolution is available even for the spin-fluctuator model. 
Different exact methods have been proposed which lead 
to approximate solutions in relevant regimes. 

In Sec. |III.B] we present approximate approaches pro- 
posed to predict dephasig due to 1// noise described 
as a classical stochastic process. Approaches based on 



the adiabatic approximation (Falci et al. 2005 Ithier 



et al. 2005 ) , allow simple explanations of peculiar non- 



exponential decay reported in different experiments with 
various setups (Bylander et al.\ |2011| iChiarello et oil 



2012 



et al.\ 



Cottet et al.\ 
2003; Sank"er 



2001 

D 



2012 



Ithier et aU, ' 20051 [Martinis 
Van Harlin gen et al.\ |2004| 



Vion et a/.[ |2002[ |Yan et a/.[ |2012[ ). These approaches 



also predict the existence of operating conditions where 
leading order effects of 1// fluctuations are eliminated 
also for complex architectures, analogously to the single 
qubit "optimal point". The effect of 1// noise in solid 
state complex architectures is a subject of current inves- 



particular in the cQED scheme, see Sec. |II.C.1[ Very re- 
cently, highly sensitive superconducting circuits could be 
used as "microscopes" for probing characteristic proper- 
ties of environmental fluctuators. The reached ability to 
"directly" control these quantum TLS has opened a new 
research scenario where TLSs may play the role of nat- 
urally formed qubit s. These issues will be addressed in 
Sec. ITTLCI 

In the final part of this Section we will present current 
strategies to reduce effects of 1// noise based on tech- 
niques developed in NMR CSch lichter[ |l992|. The open 
question about the best strategy to limit 1// noise effects 
via open or closed loop control will be discussed and we 
will review the current status of the ongoing research 
along this direction. 



A. Spin-fluctuator model 

In the following we will use a simple classical model 
according to which the quantum system - qubit - is cou- 
pled to a set of two-state fluctuators. The latter ran- 
domly switch between their states due to interaction with 
a thermal bath, which can be only weakly directly cou- 
pled to the qubit. Since we are interested only in the 
low- frequency noise generated by these switches, they 
will be considered as classical. (The situations where 
quantum effects are of importance will be discussed sep- 
arately). The advantage of this approach, which is often 
referred to as the spin-fluctuator (SF) model, is that the 
system qubit+fluctuators can be described by relatively 
simple set of stochastic differential equations, which in 
many cases can be exactly solved. In particular, many 
results can just be borrowed from much earlier papers 



on magnetic resonance (Hu and Walker, 1977 Klauder 



and Anderson, 1962; Maynard et a/., 1980), on spectral 



diffusion in glasses (Black and Halperin , 1977), as well 
as works on single-molecule spectroscopy (Barkai et al.\ 



2001; "Geva 'eFal] [19961 |Moerner[ [19941 [Moerner and Qr-| 
rit[ ,1999J . 

The interaction of electrons with two-state fluctua- 
tors was previously used for evaluation of the effects of 
noise on various systems ([Galperin and Gurevich[ [1991 



Galperin et al.\ 1989[ 1994[ Hessling and Galperinj 



JG 

m 



Kogan and Nagaev} |1984b| |Kozub| |1984| [Ludviksson 



1995 



et al.\ 1984). It was recently applied to the analysis of de- 



coherence in qubit s ( [Bergli et al. 
20041 



and Galperin, 2006 



2006 2009 



Falci et al.\ 
20051 [20031 [Galperin et a/.[ |2006[ |2004| |Martin| 



Paladino et al. 



2003b 2002). Var- 



ious quantum and non-Markovian aspects of the model 



tigation (Bellomo et al.\ 2010 Brox et a/.[|2012[ D'Arrigo were addressed by Abel and Marquardt (2008); Burkard 



Paladino et al.\ |2009[ |2011[ |2010| jStorcz et al.\ |2005[ ). 
Considerable improvement in minimizing sensitivity to 
charge noise has been reached via clever engineering, in 



erari|2008||D'Arrigo and Paladino[|2012||Hu e^ a/.[|2007| (|2009D; [Coish et a/.| (|2008D; [Culcer et al 



cenzo and Loss (2005); Galperin et al. (2003); Grishin 



de Sousa et al. (2005); and Yurkevich et al. (2010) 



et al. (2005); Lutchyn et al. (2008); Paladino et al. (2002); 



(2009); DiVin- 
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1. Exact results at pure dephasing 



In Sec. LB we discussed a simple model of so-called 
pure dephasing, when the diagonal splitting e of a qubit 
represented by the Hamiltonian ^ fl uctuates in time and 
A = 0. The resulting expression (|11[) for the FID was ob- 
tained assuming that the fluctuations obey the Gaussian 
statistics. To approach speciflc features of the 1// noise 
let us take into account that such a noise can be con- 
sidered as superposition of random telegraph processes. 
To begin with let us consider decoherence due to a single 
random telegraph process. 



a. Dephasing due to a single RT Buctuator - A random 
telegraph process is defined as follows fBuckingham' 
1989| [Kirton and Uren] 1989). Consider a classical 
stochastic variable x(t), which at any time takes the 
values x(t) = ±1. It is thus suitable for describing a 
two-state fluctuator that can flnd itself in one of two 
(meta) stable states, 1 and 2, and once a while it makes a 
switch between them. The switchings are assumed to be 
uncorrelated random events with rates 71^2 and 72^1, 
which can be different. For simplicity, we will limit our- 
selves to symmetric RT process: 71^2 = 72^1 = 7/2. 
The extension to the general case can be easily made. 



see, e. g., (Itakura and Tokura 2003). The number k of 



switches that the fluctuator experiences within a time t 
follows a Poisson distribution 



Pk{t) 



2^k\ 



-7t/2 



(44) 



The number of switches, /c, determines the number of 
times the function x{t) changes its sign contributing 
(—1)^ to the correlation function, C(t) = (x(t)x(O)). 
Therefore, 



C{t) 



-ft/2 



V(-i)^M! 



-It 



k=0 



t>0. (45) 



The RT process results in a fluctuating fleld h{t) = hx{t) 
applied to the qubit. The magnitude, b = |b|, together 
with the switching rate, 7, characterizes the fluctuator. 



Using Eq. ( 10 ) we flnd power spectrum of the noise gen- 



erated by i-th fluctuator: 



(46) 



with Cj.{uj) given by Eq. (17). This expression corre- 
sponds to the high temperature limit, E <C /c^T, of 
Eq. (I28| with xi — ^2 = 26^. Thus, toghether with the 



considerations reported in Sec. |II.B[ we conclude that 
telegraph fluctuators provide a reasonable model for the 
1// noise. 

Single shot measurements and FID - Using Eq. ([5| let us 
now discuss how a single RT fluctuator effects a qubit. 



We assume that the fluctuator does not feel any feedback 
from the qubit and thus the RT function x(^) equals to 
+1 or —1 with the probability 1/2 regardless to the di- 
rection of the Bloch vector M. One can show that under 
this assumption the probability to flnd the angle cp (see 
Fig. fl]) at time t, p{(p^t)^ satisfles the second-order differ- 
ential equation ( [Bergli et a/. , 2009): 



p^jp = bd^p, 



(47) 



which is known as the telegraph equation. We can always 
choose the x-direction such that (/p = at t = 0, so 
p{(f^O) = 5{(p). The second initial condition, 

p{^,0) = ±2bd^p{^,0), 



can be derived from the integral equation for p{(p^t)^ see 
( Bergli et fl/.[|2009" ). After averaging over the fluctuator 's 
initial state, p{(p, 0) = 0. 

The FID is given as (m+) = {e'"^) = J dipp{ip,t)e''^ . 
The differential equation and initial conditions for this 
quantity can be obtained multiplying Eq. ( [47| ) by e*'^ and 
integrating over (f. In this way we obtain: 



(m+) +7(771+) = -b^{mj^) 



(48) 



with initial conditions (m+(0)) = 1, (m+(0)) = 0. The 
solution of Eq. ( 48 ) with these initial conditions is 



M^x/T^WtF- (49) 

In the context of decoherence due to discrete noise 
sources affecting superconducting qubits this expression 



was flrstly reported in (Paladino et al. 2002) where it 



has been derived as the high-temperature limiting form 
of a real-time path-integral result. Equation (49) also 



follows by direct averaging the qubit coherence over the 
stochastic bistable process h{t) = bx(t), i. e., by evalua- 
tion of the time average Z{t) = ((exp[— z J^ dt'b{t')])). 
The last two approaches do not necessarily assume a 
thermal equilibrium initial condition for the fluctuator 
at the initial time t = 0. Thus, results also depend on 
the initial population difference between the two states 
X = ±1, ^PQ' When this quantity is flxed to one of the 
classical values +1 or —1, we obtain the effect of noise 
only during the qubit quantum evolution, i. e., in an 
"ideal" single shot measurement. The result takes the 
form given by Eq. (49) with the prefactors of the two 



exponentials, /i ± 1, replaced by /i ± 1 ^ ibp^^bj^. The 
FID results by further averaging over the initial condi- 
tions, i. e., by replacing ^po with (^po) = ^Peq, where the 
thermal equilibrium value Jpeq is consistently set to zero 
in the regime E <C /c^T. Both the single shot and the 
free induction signals demonstrate qualitatively different 
behaviors for large and small values of the ratio 26/7. 



26 



26/7 = 0.8 




FIG. 21 Time dependence if the FID amplitude for 26/7 
0.8 and 5. 



This is easily seen from Eq. ( 49 ) . At 6 ^ 7 one can con- 



sider the qubit as a quantum system experiencing beat- 
ings between the states with different splittings, Bq ± 6, 
the width of these states being 7/2. In the opposite limit- 
ing case, 6 <C 7, the energy-level splitting is self-averaged 
to a certain value, the width being 6^/7. This situation 
was extensively discussed in connection with the mag- 
netic resonance and is known as the motional narrowing 
of spectral lines (Klauder and Anderson, 1962). Different 



behaviors of the FID amplitudes depending on the ratio 



26/7 are illustrated in Fig. 21 

Thus, one can discriminate between weakly coupled 
fiuctuators 26/7 <C 1 and strongly coupled fluctuators 



(in the other regimes) ( [Paladino et al. , 2QQ2[), w hich in- 
fluence the qubit in different ways. In Fig. |22[ adapted 
from ( [Paladino et al.\ 2003b), we show the decay factor 
V{t) = — ln[|(m+)(t)|] for different values of 26/7 nor- 
malized in such a way that the Gaussian decay factor 
T^{t) = {(p'^)/2, see Eqs. ^ and ([ll}, is the same for 
all curves. We report both the single-shot and the FID 
signals. Substantial deviations are clearly observed, ex- 
cept in the presence of a weakly coupled fluctuator. In 
particular, a fluctuator with 26/7 > 1 induces a slower 
dephasing compared to an oscillator environment with 
the same Sb{oj)^ a sort of saturation effect. Recurrences 
at times comparable with 1/26 are visible in T{t). In 
addition, strongly coupled fluctuators show memory ef- 
fects. This is clearly seen considering different initial 
states for the fluctuator, corresponding to the single-shot 
and FID measurement schemes. This is already seen in 
the short-times behavior, 7^ <C 1, relevant for quantum 
operations which need to be performed before the sig- 
nal decays to a very low value. In the limit 7^ <C 1, 
r^(t) ^ 6^(1 - 5pl^)t^/2, whereas 



m 



_M^2. 



2SpoSp, 



'eq 



6 



_M^,3' 



(50) 



In a single-shot process, 5po = ±1, and T{t) oc t^, show- 
ing that a fluctuator is stiffer than a bath of oscillators. 
On the other side, for repeated measurements 5po = Spe^^ 




FIG. 22 Reduced r(t) due to fluctuators prepared in a stable 
state (6po — 1 solid lines and 5po — —1 dotted lines) and in 
a thermal mixture {5pQ = dashed lines) for the indicated 
values of ^ = 26/7. Inset: longer time behavior for stable 
state preparation. The curves are normalized in a such a way 
that the oscillator approximation for all of them coincides 

|2003b| ), 



(thick dashed line). Adapted from (Paladino et al. 



with kind permission of Springer Science-hBusiness Media 



for very short times we recover the Gaussian result. This 
is due to the fact that inhomogeneous broadening due 
to the uncontrolled preparation of the fluctuator at each 
repetition adds to the effect of decoherence during the 
time evolution and results in a faster decay of V(t). 



The difference between Eq. (49) and the result (12) 
based on the Gaussian assumption is further eluci- 
dated considering the long-time behavior. Substituting 



Eq. (46) for the noise spectrum into Eq. (12) one obtains 



l/T, 



*(G) 



bVi. 



(51) 



This result should be valid for times much longer than 
the correlation time of the noise, which is 7"^. Expand- 
ing Eq. ( [49| ) at long times we flnd that FID also de- 
cays exponentially (or the beatings decay exponentially 
at 6 > 7). However, the rate of decay is parametrically 
different from Eq. (IsTl): 



1 



w-'i-% 



(52) 



At 6 <C 7 E g. ([ 52| coincides with the Gaussian result. 
Shown in Fig.|23|are the dephasing rates I/T2 and I/T2 
given by equations (51) and (52), respectively. Again 
we see that the Gaussian approximation is valid only in 
the limit 6 <C 7. The main effect of a strongly coupled 
fluctuator is a static energy shift, the contribution to the 
qubit decoherence rate saturates at ~ 7 and at 6 > 7/2 
the Gaussian assumption overestimates the decay rate. 
Apparently, this conclusion is in contradiction with the 



discussion in Sec. LB following from the central limiting 
theorem. According to this theorem p{(f>^ t) always tends 
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FIG. 23 Comparison of the dephasing rate T2^ for a single 
random telegraph process and the corresp onding Gaussian 
approximation. Adapted from (Bergli et al. 2009 J. 



to a Gaussian distribution with time-dependent variance 
provided that the time exceeds the correlation time of 
the noise. To resolve this apparent contradiction, let us 
analyze the shape of the distribution function, p{(f,t)^ 



following from the telegraph equation (47) with initial 
conditions p{ip,0) = S{(p)^ p{(p^O) = 0. This solution is 



(l/2)e-^^/2 [S{^ + bt) + S{^ - bt)] 
+ [7/6^(t)]e-^^/' [e((^ + bt) - e((^ - bt)] 
X {h[u{t)jt/2] + u{t)IoHt)-ft/2]} . (53) 

is the modified Bessel function, u{t) = 
{(p/bt)'^^ while Q{x) = 1 at x > and at x < 



p{(f,t) 



I^{x) 



Here 

is the Heaviside step function. This distribution for var- 



ious t shown in Fig. |24] consists of two delta-functions 
and a central peak. The delta- functions represent the 
finite probability for a fluctuator to remain in the same 
state during time t. As time increases, the weight of the 
delta- functions decreases and the central peak broadens. 
At long times, this peak acquires a Gaussian shape. In- 
deed, at 7^ ::^ 1 one can use asymptotic form of the 
Bessel function, I(3{z) ~ (27r2:)~"'-/^e^, as 2: ^ oo. In 
addition, at t ^ (p/b we can also expand ^/l — {(p/bt)'^ 
and convince ourselves that the central peak is described 
by the Gaussian distribution ^ with {cp'^) = 2b'^t/^. If 
the qubit-fluctuator coupling is weak, 6 <: 7, this Gaus- 
sian part of p{(f,t) dominates the average (e*'^) and the 
Gaussian approximation is valid. Contrary, when the 
coupling is strong, b > 7/2, the average is dominated by 
delta- functions at the ends of the distribution and the de- 
coherence demonstrates a pronounced non- Gaussian be- 
havior, even at long times {t > 2/7). 

Unfortunately, we are not aware of a way to measure 
the distribution p{(p^ t) in experiments with a single qubit. 
The reason is in the difference between the qubit that can 
can be viewed as a pseudospin 1/2 and a classical Bloch 
vector M. According to Eq. Q the components M^^ My^ 
Mz of M are connected with the mean component of the 
final state of the pseudospin. Therefore, to measure the 



FIG. 24 The distribution ([53} for 26/7 = 1 and 7t/2 = 1, 3, 
and 5 (numbers at the curves). Only the part with positive 
(p is shown; the function is symmetric. The arrows represent 



the d elta-functions (not to scale). Adapted from ( [Bergli et al. 
20091 



value of the phase cp (argument of 771+) that corresponds 
to a given realization of the noise one should repeat the 
experimental shot with the same realization of the noise 
many times. This is impossible because each time the 
realization of noise is different. Therefore, the only ob- 
servable in decoherence experiments is the average (e*'^). 

b. Echo - The analysis of the echo signal is rather sim- 



ilar: one has to replace (m+( ^t)) taken from Eg . (49) by 
(ml^^(2ri2)) where, cf. with (Laikhtman 1985), 

p-7t/2 



{m^;\t)) 



2v^ 






(54) 



This result can also be obtained by direct calculation of 
the function Z(2t\r]) = ((exp [-i J^^ dt'r]{t')b{t')])) where 
r^(t') = 1 for < t' < t and r](t') = -1 for t < t' < 2t 
( Falci et fl/^|2003 ). To demonstrate non-Gaussian behav- 
ior of the echo signal let us evaluate from Eq. ( 14 ) the 
variance ('0^(2ri2)). We obtain, cf. with (Bergli et al. 
|2QQ9l ), 

(^') = (26V7') (27^2 - 3 + 4e-^"^2 - e-'^"^2) (55) 

with {rrv^ ) = e~^^ ^/^ as before. The comparison be- 
tween the two results is shown in Fig. [25] for a weakly 
(26/7 = 0.8) and a strongly (26/7 = 5) coupled fluctu- 
ator. As discussed above, the Gaussian approximation 
is accurate at 2b < 7, while at 26 > 7 the results are 
qualitatively different. In particular, the plateaus in the 
time dependence of the echo signal shown in Fig. [25] are 
beyond the Gaussian approximation. We believe that 
such a plateau was experimentally observed by Naka-] 
mura et al. (2002) (see Fig. 3 from that paper partly 
reproduced in Fig. p|. In the limit 6^7, y/j/r equa- 
tion (54) acquires a simple form: 



.(e) 



{m);\r)) = e-^^ [1 + (7/26) sin26r] 



(56) 
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FIG. 25 Solid lines - echo signal for different values of the ra- 
tio 26/7, Eq. ( 54 ). Dashed lines - calculations along the Gaus- 
sian approximation, Eq. (55). Adapted from ( Bergli et al.\ 
[20091. 



According to this expression, the plateau-like features 
(|9(m!^^)/9r| <C 1) occur at br ^ kir (where k is an in- 
teger) and their heights {rrv^^) ^ ^—^^-f/b exponentially 
decay with the number k. Measuring experimentally the 
height and the position of the first plateau, one can de- 
termine both the fluctuator coupling strength b and its 
switching rate 7. For example, the echo signal measured 



by Nakamura et al. (2002) shows a plateau-like feature 
at T12 = 3.5 ns at the height (m^^) = 0.3, which yields 
b ^ 143 MHz and 7/2 ^ 27 MHz. If the fluctuator 
is a charge trap near a gate producing a dipole electric 
field, its coupling strength is 6 = e^(a • r)/r^. Using the 
gate-CPB distance r = 0.5 /im, we obtain a reasonable 
estimate for the tunneling distance between the charge 
trap and the gate, a ~ 20 A. A more extensive discussion 



can be found in (Galperin et al. 20061). A similar analy- 



sis for arbitrary qubit working point has been reported in 



(Zhou and Joynt, 2010). The effect of RTN for arbitrary 



working point will be disussed in Sec. HI. A. 2 



c. Telegraph noise and Landau-Zener transitions - Driven 
quantum systems are exceedingly more complicated to 
study than stationary systems, and only few such prob- 
lems have been solved exactly. An important excep- 



tion is the Landau-Zener (LZ) transitions (Landau, 1932 
Stueckelberg , 1932 Zenerj Il932| . In the conventional LZ 



problem, a TLS is driven by changing an external param- 
eter in such a way that the level separation he is a linear 
function of time, e{t) = o?t. Close to the crossing point of 
the two levels, an interlevel tunneling matrix element h/S. 
lifts the degeneracy in an avoided level crossing. When 
the system is initially in the ground state, the probabil- 
ity to find it in the excited state after the transition is 
g-7rA /2a Hcucc, a fast rate drives the system to the 
excited state, while the system ends in the ground state 
when driven slowly. 



In connection with decoherence of qubits, there has 
recently been increased interest in Landau-Zener transi- 
tions in systems coupled to an environment. This prob- 
lem is both of theoretical interest and of practical impor- 



tance for qubit experiments (Sillanpaa et al., 2006). The 



noisy Landau-Zener problem has been discussed by sev- 



eral authors (Kayanuma', 1984, 1985; Nishino et a/.[|2001 
Pokrovsky and Sinitsyn' 2003 ; Pokrovsky and Sun( '2007 



Saito et a/. ,2007; Shimshoni and Gefen, 1991; Shimsh^m] 
and Stern[ |1993[ |Wubs et al.\ |2006|. Here we will dis- 
cuss the role of a classical telegraph fluctuator following 
Vest garden et al.\ (|2008 ) . 

Let us model a driven qubit by the Hamiltonian (|3| 
with e{t) = a^t and A(t) = A + vx(^) where x(^) is a RT 
process. Representing the density matrix as in Eq. (pi) 
one can study the dynamics of the Bloch vector averaged 
over the realizations of the telegraph noise. To this end, 
one has to study the solution of the Bloch equation ([5| 
in the presence of the "magnetic fields" B± = Bq ± v(t). 
Introducing the partial probabilities p-i-(M,t) to be in 
state M at time t under rotations around B+ and B_, 
respectively and defining 

M± = [ d^M[p^{M,t)±p-{M,t)] 



one arrives at the set of equations, cf. with ( Vestgarden 
et a/.l[2008) 



M+ = M+ X Bo 



V X M_ , 

M_ = M_ X Bo + V X M+ - 7M_ . (57) 



Here M+ is the average Bloch vector whose z-component 
is related to the occupancy of the qubit 's levels. To 
demonstrate the role of telegraph noise, let us consider 
a simplified problem putting A = 0. In this way we will 
study the LZ transitions induced solely by the telegraph 



process. In this case, the Eqs. (57) reduce to an integro- 
differential equation for M^ 



MS) 



M+, 



a^ti{2t-ti] 



e-Tl*ilM^(t-ti). 
(58) 



This equation can be easily analyzed for the case of fast 
noise, 7 ::^ a. A series expansion in ti leads to the follow- 
ing solution of Eq. (ISSl) with initial condition Mz{^) = 1: 



(59) 



M^(oo) = e-^^'/" 



This is the usual expression for the LZ transition proba- 
bility with replacement A -^ \f2v. This result holds for 
any noise correlated at short times (<C «~^), as shown 
byjPokrovsky and Sinitsyn (2003). Note that the result 
(59) can be obtained using the Gaussian approximation. 



The case of a "slow" fluctuator with 7 < a leads to very 
different results which, for an arbitrary ratio vj a^ require 
a numerical solution of the integro-differential equation 
(58). The result of such analysis is shown in Fig. 26 
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FIG. 26 (Color online) The Mz (oo) versus the ratio v/a for 
different switching rate 7 (shown in units of a) . Adapted from 
(Vest garden et a/. [[2008). 



adapted from (Vestgarden et al., 2008). The static case 
7 = is equivalent to the standard Landau- Zener tran- 
sition with the noise strength \/2v replacing the tunnel 
coupling A between the diabatic levels. In the adiabatic 
limit, V ^ a, we see that Mz{oo) -^ —1, which corre- 
sponds to the transition to the opposite diabatic state. 
In this case the dynamics is fully coherent, but since we 
average over the fluctuator's initial states, the Bloch vec- 
tor asymptotically lies on the z axis of the Bloch sphere. 
At finite 7, the noise also stimulates Landau-Zener tran- 
sitions, however the transition probability decreases with 
increasing switching rate. Note that the curves cross the 
line Mz{oo) = at some v/a^ which depends on 7. For 
this noise strength, the final state is at the center of the 
Bloch sphere, corresponding to full decoherence. How- 
ever, increasing the noise strength beyond this point re- 
sults in a final state with negative Mz{oo). Thus, we have 
the surprising result that under some conditions, increas- 
ing the noise strength will also increase the system pu- 
rity after the transition. The results for slow and strong 
noise cannot be obtained in the Gaussian approximation. 
These considerations point out that telegraph noises can 
facilitate LZ transitions. This processes may prevent the 
implementation of protocols for adiabatic quantum com- 
puting. 



d. Ensemble of Buctuators: effects ofl/f noise - As a re- 
sult of the above considerations, we conclude that the 
role of a fluctuator in decoherence of a qubit depends 
on the ratio between the interaction strength, 6, and the 
correlation time, 7"^, of the random telegraph process. 
Weakly coupled fluctuators, i. e., "weak" and relatively 
"fast" fluctuators for which 26/7 <C 1, can be treated as 
a Gaussian noise acting on the qubit. Contrary, the in- 
fluence of strongly co^/p/ec? fluctuators, i. e., "strong" and 
"slow" fluctuators for which 26/7 > 1, is characterized 
by quantum beatings in the qubit evolution. 

As we have seen, a set of fluctuators characterized by 
the distribution function V{j) oc I/7 of the relaxation 
rates provides a realistic model for 1 //-noise. Therefore, 



it is natural to study in which way the qubit is deco- 
hered by a sum of the contributions of many fluctuators, 
^(^) = ^i^iXi{^)' A key question here is "Which fluc- 
tuators - weak or strong - are responsible for the qubit 
decoherence?" In the first case the noise can be treated 
as Gaussian while in the second case more accurate de- 
scription is needed. 

Here we will analyze this issue using a simple model 
assuming that the dynamics of different fluctuators are 
not correlated, i. e., (Xi(^)Xj(^O) = (^ije~^*l^~^ L Under 
this assumption, the average (m+) is the product of the 
partial averages, 

(m+(t)) = [](m+,(t)) =exp f ^ln(m+,(t)) J . (60) 

Since the logarithm of a product is a self- averaging quan- 
tity, it is natural to approximate the sum of logarithms, 
^^ ln(7Ti+i(t)) by its average value, — /C^(t), where 



^m{t) = -^ln(m+i(t)) . 



(61) 



Here, the bar denotes the average over both the cou- 
pling constants, 6, of the fluctuators and their transi- 
tion rates, 7. This expression can be further simplified 
when the total number A/r of thermally excited fluctu- 
ators is very large. Then we can replace ^^ln(m+^(t)) 
by A/Tln(m+(t)). Furthermore, we can use the so-called 
Holtsmark procedu re ([Chandresekhar 1943), i. e., to re- 
place ln(r?i+(t)) by (m+) — 1 assuming that in the rele- 
vant time domain each (m+^) is close to 1. Thus, ICm{t) 



is approximately equal to (Galperin et al. 2003 [20041 



Klauder and Anderson[ p!962[ [Laikhtman 1985)^ 



lCm{t)^MT l-(m+(t)) 



/ 



dbd-iV{h,-f)[l-{m+{h,-i\t))]. (62) 



Here (m+^) is specified as (7n+(6, 7|t)), which depends on 
the parameters h and 7 according to Eq. (49). The free 
induction signal is then exp[— /C^(t)]. Analysis of the 
echo signal is rather similar: one has to replace {m-^{t)) 
taken from Eq. (|49| with Eq. ^. 



To evaluate the time dependence of the free induc- 
tion or echo signal, one has to specify the distribution 
function of the coupling constants 6, i. e., partial contri- 
butions of different fluctuators to the random "magnetic 
field" b{t). To start with, we consider the situation where 
the couplings bi are distributed with a small dispersion 
around an average value b. Under these conditions the 
total power spectrum reads 



S{uj) = 62 



/ 



dj -^ C^{uo) 



A 



Irr^ 27 ' • • UJ 



(63) 



where the amplitude A can be expressed in terms of 
the number of fluctuators per noise decade, Ud = 
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ArTln(10)/ln(7M/7m), as follows A = b^Po/4. = 
6^71(^/(2 In(lO)). The spectrum exhibits a crossover to 
cj"^ behavior at a; ~ jm- In Fig. 27 we show the 
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FIG. 27 Saturation effect of slow fluctuators for a 1// spec- 
trum and coupling distributed with {Ab)/{b) = 0.2. Labels 
indicate the number of decades included. Adapted from ( |Pal-| 
adino et a/. 1 120021). 
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r(t)/r2(t) for 1// spectrum 



2x10^ Hzand7M/(27r) = 2x10^ Hz with 
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different numbers of fluctuators per decade: (a) Ud 
nd = 4 X 10^, (c) rid = 8 X 10^, (d) and (dl) n^ = 4 x 10 
(e) and (el) n^ = 4 x 10^. Solid lines correspond to 6poj — 
±1, dashed lines corre spond to equilibrium initial conditions 
(FID). Adapted from ([Paladino erZj|2002). 



results for a sample with a number of fluctuators per 
decade Ud = 1000 and with couplings distributed around 
S/(27r) = 4.6 x 10^ Hz. Initial conditions dpoj = ±1 are 
distributed according to (dpoj) = dp^^^ the equilibrium 
value. The different role played by weakly and strongly 
coupled fluctuators is illustrated considering sets with 
7M/(27r) = 10"*^^ Hz and different 7^. In this case the de- 
phasing is given by fluctuators with 7^ > 27r x 10^ Hz ^ 
26/10. The main contribution comes from three decades 
at frequencies around ^ 2b. The overall effect of the 
strongly coupled fluctuators (7^ < 26/10) is minimal, de- 
spite their large number. We remark that, while satu- 
ration of dephasing due to a single fluctuator is physi- 
cally intuitive, it is not a priori clear whether this holds 
also for the 1// case, where a large number (~ I/7) of 
slow (strongly coupled) fluctuators is involved. The de- 
cay factor }Cm{t) can be easily compared with the Gaus- 
sian approximation ( 1 1 ) ( Paladino et at. , 2002| ) . This ap- 
proximation fails to describe fluctuators with 26/7 ^ 1. 
For instance, {^p^) at a fixed t scales with the number of 
decades and does not show saturation. On the other side 
the Gaussian approximation becomes correct if the en- 
vironment has a very large number of extremely weakly 
coupled fluctuators. This is shown in Fig. [28] where the 
power spectrum is identical for all the curves but it is 
obtained by sets of fluctuators with different Ud and h^. 
The Gaussian behavior is recovered for t ^ l/7m if ^d is 
large (all the fluctuators are weakly coupled) . If in addi- 
tion we take (dpoj) = dpe^^^ ICm{t) approaches {(j)'^) also at 
short times. Hence decoherence depends separately on Ud 
and 5, whereas in the Gaussian approximation only the 
combination ridb'^, which enters 6'(cj), matters. In other 
words, the characterization of the effect of slow sources 
of 1// noise requires knowledge of moments of the bias 
fluctuations higher than S{uj). A more accurate aver- 



aging procedure considering (Spoj) = ±1 for fluctuators 
with jitjn < 1 and {Spoj) = Spe^^ for jitm > 1, shows 
that, for long measurement times btm ^ 1, fluctuators 
with 7 < 1/tm are saturated and therefore not effective, 
whereas other fluctuators, being averaged, give for short 
enough times 



m 



duo S{uj) 



l/tr. 



smujt/2 
uj/2 



(64) 



This is illustrated in Fig. 29 The above considerations 
explain the often used approximate form ^^ obtained 
from Eq. (11) with low- frequency cut-off at l/tm which 
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FIG. 29 Different averages over Spoj for 1// spectrum repro- 
duce the effect of repeated measurements. They are obtained 
by neglecting (dotted lines) or account for (solid lines) the 
str ongly correlated dynamic s of 1// noise. The noise level 



of (Nakamura et al. 
9.2 X 10° Hz, Ud 



2002| ) is used, by setting \b\/{27i) 

10^ Hz. 



10^ 7m/(27r) = 1 Hz, 7M/(27r) 
Dashed lines are the oscillator approximation with a lower 



cutoff at cj = min{|6|, 1/tm}- Adapted from (Paladino et al. 



2003b), with kind permission of Springer Science+Business 
Media. 
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was proposed in ("Cottet et al , 2001 ) and that we already 
mentioned in Section |ILC[ Eg. (32). 

The echo signal obtained from Eq. ( [62| ) with (54) as- 
suming that 7^(6,7) (X j~^5{b — 6)0(7^ — 7) leads to 
the following time dependence when b <C 7m (Galperin 



et al 2007) 



/C^^^t) 



\n{2)At^ 
a{A/b^)bt 



t << 7m' 

7m' < ^ < ^"' 



(65) 



<Ct 



where a ~ 6. At small times, the echo signal behaves as 
in the Gaussian approximation, Eq. (14), where the t^ de- 



pendence follows from the crossover of the spectrum from 
uj~^ to cj"^. On the other side, at large times b~^ <^ t 
the SF model dramatically differs from the Gaussian re- 
sult which predicts /C[?\t) = ^ln(2)t^ The origin of 
the non-Gaussian behavior comes from the already ob- 
served fact that decoherence is dominated by the fluc- 
tuators with 7 ^ 2b. Indeed, very "slow" fluctuators 
produce slowly varying fields, which are effectively refo- 
cused in the course of the echo experiment. As to the 
"too fast" fiuctuators, their infiuence is reduced due to 
the effect of motional narrowing. Since only the fluctua- 
tors with 6 <C 7 produce Gaussian noise, the noise in this 
case is essentially non-Gaussian. Only at times t <C 6~^, 
which are too short for these most important fluctuators 
to switch, the decoherence is dominated by the faster 
fluctuators contributions, and the Gaussian approxima- 
tion turns out to be valid. Instead when b ^ 7m all 
fluctuators are strongly coupled and the long-time echo 



decay is essentially non-Gaussian (Galperin et a/., 2007) 






t-€:b-\ 
t:s>b-\ 



(66) 



Thus, the long-time behavior of the echo signal depends 
on the details of the noise model, the dependence on 
the high-frequency behavior of the spectrum being the 
flrst manifestation. Unfortunately, the measured signal 
is usually very weak at long times, making it diflicult 
to flgure out the characteristics of the noise sources in 
the speciflc setup. Indeed, echo signal data measured 
on a flux qubit in the experiment by |Yoshihara et c 
( 2006 ) can be equally well flt assuming either a Gaussian 



statistics of the noise or a non-Gaussian model. The two 
models can in principle be distinguished analyzing the 
different dependence on the average coupling b. Details 
of this analysis are given in (Galperin et a/., 2007) where 
a flt with the SF model for different working points in- 
dicates that equations ( 65 ) for the case b <C 7m give an 



overall better flt. A similar experiment on a flux qubit 



has been reported in (Kakuyanagi et a/., 2007). Data for 



the echo decay rate in that case have been interpreted as- 
suming Gaussian fluctuations of magnetic flux, and con- 
sistency with this model has been observed by changing 
the qubit working point. In general, echo procedures 



allow extraction of relevant information about the noise 



spectrum, like the noise amplitude ([Galperin et al. , 2007 



Kakuyanagi et a/.]|2007[[Yoshihara et a/.[|20Q6[Zhou and 



Joyntl 2010), the noise sensitivity deflned by Eq. (90) 



( Yoshihara et al. , 2006 ) or the average change of the flux 



in the qubit loop due to a single fluctuator flip (Galperin 
et al. , 2007V Similar possibilities can be provided by the 



real time qubit tomography (Sank et a/., 2012). 



We now consider the role of the interaction strengths 
distribution, bi. Using the same model as for TLS in 



glasses (see Sec. II. A) we consider each fluctuator as a 
TLS with partial Hamiltonian 



n 



jii) 



1 



{UiTz -^ KTx) , 



(67) 



where r^ is the set of Pauli matrices describing z-th 
TLS. The energy- level splitting for this fluctuator is 
Ei = \/Uf + A| . Fluctuators switch between their states 
due to the interaction with the environment, which is 
modeled as thermal bath. It can represent a phonon 
bath as well as, e. g., electron- hole pairs in the conduct- 
ing part of the system. Fluctuations of the environment 
affect the fluctuator through the parameters Ui or A^. 
Assuming that fluctuations of the diagonal splitting U 
are most important, we describe fluctuators as TLSs in 
glasses, see Sec. |II.A| 

For the following, it is convenient to characterize fluc- 
tuators by the parameters Ei and 6i = arctan(A^/[/^). 
The mutual distribution of these parameters can be writ- 
ten as ([Laikhtman 1985) 



P{E,0) = Po/sin(9, < 6> < 7r/2 , 



(68) 



which is equivalent to the distribution (27) of relaxation 
rates. To normalize the distribution one has to cut it 
off at small relaxation rates at a minimal value 7min or 

Jlo < 1. 



cut off the distribution (68) at 0^ 



7n 



The distributions given by Eqs. (27) and (68) lead to the 
(xl/uj noise spectra at 7min <C cj <C 70. 

The variation of the qubit's energy- level spacing can be 
cast in the Hamiltonian, which [after a rotation similar 
to that leading to Eq. (25)] acquires the form 



K 



qF 



n^6,a,ri^ 



bi = g{ri)A{ni) cos Oi. (69) 



Here n^ is the direction of the elastic or electric dipole 
moment of i-th fluctuator, and r^ is the distance between 



the qubit and the i-th fluctuator. Note that in Eq. (69) 
we neglected the term ex a^r^^. This can be justifled 
as long as the fluctuator is considered to be a classical 
system. The functions A(n^) and g{ri) are not universal. 
The coupling constants, 6^, deflned by Eq. (69) contain 
cos Oi. Therefore they are statistically correlated with 
Oi. It is convenient to introduce an uncorrelated random 
coupling parameter, Ui^ as 



Ui = g{ri)A{ni), bi =UiCosi 



(70) 
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FIG. 30 (Color online) Sketch of localized charges near an 
electrode in a Cooper Pair box. Induced image charges create 
local dipoles that interact with the qubit. 



Let us also assume for simplicity that the direction n^ of 
a fiuctuator is correlated neither with its distance from 
the qubit, r^, nor to the tunneling parameter represented 
by the variable Oi. This assumption allows us to replace 
A{ni) by its average over the angles, A= (A(|n|))n. 

Now we are ready to analyze the decoherence using 
Eq. (62). The question we will address in this Section 



is whether a special group of fluctuators responsible for 
decoherence does exist. 

An interesting feature of the problem is that the result 
strongly depends on the decay of the coupling parameter, 
^(r), with the distance r, which usually can be described 
by a power law: g{r) = g/r^. To illustrate the point, let 
us compare two cases: (i) the fluctuators are distributed 
in three-dimensional space {d = 3) and (ii) the fluctuators 
are located in the vicinity of a two-dimensional manifold, 
e. g., in the vicinity of the interface between an insulator 
and a metal {d = 2). Using the distribution (68) of the 



relaxation rates one can express the distribution V{u^O) 
as 



v{u,e) 



u(d/s)+isin6' 



, n 



9 



ad 



rf ""^ - {P^kBTY'"' 



(71) 
Here we have assumed that only the fluctuators with 
E ^ ksT are important for the decoherence since the 
fluctuators with E ^ ksT are frozen in their ground 
states. The typical distance between such fluctuators is 
tt = ad{PokBT)~^^^ where ad is a d-dependent dimen- 
sionless constant. In the following, we will, for simplicity, 
assume that 



rmin < Tt < Tmax 



(72) 



where rmin (^max) are the distances between the qubit 
and the closest (most remote) fiuctuator. In this case, 
T] oc T^/^ is the typical value of the qubit-to-fluctuator 
coupling playing role of the coupling parameter b. As 



soon as the inequality (72) is violated the decoherence 
starts to depend on either rmin or rmax, i- e., becomes 
sensitive to particular mesoscopic details of the device. 

Let us first consider the case when d = s, as it is for 
charged traps located near the gate electrode (where s = 



d = 2 due to the dipole nature of the field produced by 



a charge and its induced image, as shown in Fig. 30). In 
this case one can rewrite Eq. (62) as 



ICfit) = V 



du 



7r/2 



d(9 tan i9/(ii COS (9,70 sin^i9|t). (73) 



Here f{v^ ^\t) is equal either to 1 — (771+ ('U, ^\t)) or to 1 — 
(ttiI^ ('U,7|t)), depending on the manipulation protocol. 
Equation ^2\ together with Eqs. (49) and (54) allows 



one to analyze various limiting cases. 

As an example, let us consider the case of the two- 



pulse echo. To estimate the integral in Eq. (73) let us 



look at asymptotic behaviors of the function / following 



from Eq. (54): 



r(e) 



t^7osin6>(iicos6>)^, t <C (70 sin6>)" 



ex 



tu COS 0, 



(70 sin l9) - ^<Ct<C^i- 



<Ct. 



(74) 

Note, that these dependences differ from those given by 
Eqs. (65) and (66). Splitting the regions of integration 



over u and 6 according to the domains ([74|) of different 
asymptotic behaviors, one obtains ( [Galperin et a/.[|2004[ ): 



/C^''^(2ri2) -- r]Ti2 min{7ori2, 1} . 



(75) 



The dephasing time (defined for non-exponential decay 
as the time when /C ~ 1) for the two-pulse echo is thus 
equal to 



maxJT/ \(7/7o) ^^^} . 



(76) 



The result for 70 n 2 ^ 1 has a clear physical meaning, 



cf. with (Laikhtman 1985): the decoherence occurs only 
provided at least one of the fluctuators flips. Each flip 
provides a contribution ^ rjt to the phase while 70T12 is a 
probability for a flip during the observation time ~ ri2. 
The result for long observation times, 70T12 ^ 1, is 
less intuitive since in this domain the dephasing is non- 
Markovian, see (Laikhtman 1985) for more details. In 



this case the decoherence is dominated by a set of opti- 
mal (most harmful) fluctuators located at some distance, 
ropt(^), from the qubit. This distance is determined by 
the condition 



^(ropt)^7o(T). 



(77) 



Though derivation of this estimate is rather tedious, see 



(Galperin et a/., 2004) for details, it emerges naturally 



from the behavior of the decoherence in the limiting cases 
of strong (h ^ 7) and weak {h <C 7) coupling. For strong 
coupling, the fluctuators are slow and the qubit 's behav- 
ior is determined by quantum beatings between the states 
with E ±b. Accordingly, the decoherence rate is ^ 7. In 
the opposite case, as we have already discussed, the de- 
coherence rate is ~ b'^/j- Matching these two limiting 



cases one arrives at the estimate (77). 
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What happens is d ^ s7 If the couphng decays as r~^ 
and the fluctuators are distributed in a (i-dimensional 
space, then r^~^dr is transformed to du/v}^^^^^\ There- 
fore, V{u) oc l/u^^^^^^\ As a result, dX d < s the deco- 
herence is controhed by the optimal fluctuators located 
at the distance ropt provided they exist li d < 5, but 
the closest fluctuator has 6max ^ To, then it is the quan- 
tity 6max that determines the decoherence. M d > s 
the decoherence at large time is dominated by most re- 
mote fluctuators with r = Tmax- In the last two cases, 
/C(t) oc t^, and one can apply the results of Paladino et at. 
( 2QQ2| ), substituting for h either h^in or 6max- 

Since ropt depends on the temperature, one can ex- 
pect crossovers between different regimes as a function 
of temperature. A similar mesoscopic behavior of the de- 
coherence rates was discussed for a microwave-irradiated 



\i^^~ 11^ 



Andreev interferometer ( jLundin and Galperin 2001). It 
is worth emphasizing that the result (75) for the long- 



range interaction cannot be reproduced by the Gaussian 
approximation. Indeed, if we expand the Gaussian ap- 
proximation, e~^^ ^/^, with {ip'^) given by Eq. (55) in 



the same fashion as in Eq. (74) and then substitute to 



Eq. ( 73 ) the resulting integral over u will be divergent at 
its upper limit. Physically, this would mean the domi- 
nant role of nearest neighbors of the qubit. At the same 
time, the SF model implies that the most important fluc- 



tuators are those satisfying Eq. (77). 



The existence of selected groups of fluctuators, out of 
an ensemble of many, which are responsible for decay of 
speciflc quantities, is closely related to the self- averaging 
property of the corresponding decay laws. In (|Schriefl| 



et al. , 2006 ) a class of distribution functions of the form 



(71) with < d/s < 2 has been analyzed. In this case 
the average over the coupling constants of the Lorentzian 
functions ( 46 ) diverges at the upper limit and the noise is 



dominated by the most strongly coupled fluctuators. The 
free induction decay has been found non-self-averaging 
both at short and long times with respect to 7q"^. Non- 
self- averaging of the echo signal for 2ti2 ^ 7o^^ has also 



been demonstrated (Schriefl et al., 2006). 



e. Dephasing according to other phenomenological models - 
We remark that, beside the spin-fluctuator model, other 
stochastic processes have been considered for the descrip- 
tion of fluctuations having 1// spectral density. In this 
paragraph we mention an alternative phenomenological 
model, which has been discussed in connection with the 
problem of qubit dephasing. The motivation to consider 
these models comes from the observation of asymmetric 
telegraphic signals, with longer stays in the "down" state 



than in the "up" state, reported in tunnel junctions (Zim- 



merli et a/., 1992), in MOS tunnel diodes (jBuehler et al 



2004|), and of a spike fleld detected in a SET electrome- 
ter ( |Zorin et al. 1996). In those cases the signal exhibits 
1// spectrum at low frequencies. This suggests a descrip- 
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FIG. 31 (a) Low frequency noise as sum of RT signals. The 
switching rates for the "up" and "down" states are compa- 
rable 7+ ^ 7- and 1// noise results from a superposition 
of RT signals distributed as oc I/7. In (b) the intermittent 
limit corresponds to the limit where the noise stays in the 
down states most of the time (7+ ::^ 7_ ) . The intermittent 
noise can be approximated by a spike field with independent 
random waiting times and spikes eighths. It gives a nonsta- 
tionary 1//^ spectrum depending on the plateaus distribution 
function. Adapted from ([Schriefl et a/.|,2005aj 



tion of the 1// noise in terms of a single asymmetric RT 
signal. In (Schriefl et a/.| |2005a|b|) a phenomenological 
model for a 1//^ classical intermittent noise has been 
considered. The model can be viewed as the intermit- 
tent limit of the sum of RT signals where the duration of 
each plateau of the RT signal, Tav, is assumed to be much 
shorter than waiting times between the plateaus. In this 
limit, the noise is approximated by a spike fleld consist- 
ing of delta functions whose heights x follow a Gaussian 
distribution or, more generally, a distribution with flnite 
flrst and second moments, x and s = (x^)^/^, respec- 



et 



tively, as illustrated in Fig. 31 adapted from (Schriefl 
|2005a| ). The variance s plays the role of a cou- 



pling constant between the qubit and the stochastic pro- 
cess. A 1// spectrum for the intermittent noise is re- 
covered for a distribution of waiting times r behaving as 
r~^ at large times. Because the average waiting time 
is inflnite, no time scale characterizes the evolution of 
the noise which is non-stationary. In the pure dephasing 
regime the relative phase between the qubit states per- 
forms a continuous time random walk (CTRW) ( |Haus 
and Kehr 1987) as time goes on. Using renewal the- 
ory ( Feller[ 1962) authors flnd an exact expression for 
the Laplace transform of the dephasing factor allowing 
analytical estimates in various limiting regimes and ad- 
dress the consequences of the CTRWs nonstationarity. 
The noise is initialized at time t' = and the coupling 
with the qubit is turned on at the preparation time, tp. 
The two qubit states accumulate a random relative phase 
between tp and t -\- tp and nonstationarity manifests it- 
self in the dependence of dephasing factor, on the age 
of the noise, tp. This analysis shows that two dephasing 
regimes exist separated by a crossover coupling constant, 
Sc{tp) = [2(ravAp)ln(tp/rav)]"^/^, strongly dependent on 
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FIG. 32 (a) Qubit Bloch sphere. An isolated qubit de- 
fines the mixing angle = arctanA/e, 1-L± define 0± = 

arctanA/(e ± 2b). (b) Qubit energy bands ±\/E'^ + A^: 
th e energy split tings depend on the impurity state, Q± = 

v/(e±6)2 + A2. 



the preparation time tp. For 5 < Sc{tp) dephasing is expo- 
nential and independent on tp, whereas when s > Sc{tp) 
the dephasing time depends algebraically on tp. Since 
Sc{tp) decays to zero with the noise age, any qubit cou- 
pled to the non-stationary noise will eventually fall in the 
regime s > Sc{tp). 

We remark that in Schriefi et al. ( 2QQ5a|b ) the dephas- 
ing factor is defined as a configuration average over the 
noise, rather than a time average in a given configura- 
tion. The two averages do not coincide in general for 
nonstationary or aging phenomena. Thus one should be 
careful in comparing these results with experiments. 



2. Decoherence due to the SF model at general working point 

Here we discuss decoherence due to the spin-fluctuator 
model at general operating point where 

Htot = - i^xcr^ + n,{q)a,] + - bx{t)a, . (78) 

As a difference with the pure-dephasing regime, no ex- 
act analytic solution is available under these conditions. 
Different approaches have been introduced to study the 
qubit dynamics leading to analytical approximations in 
specific limits and/or to (exact) numerical results. 

A straightforward approach consists in solving the 
stochastic Schrodinger equation. Using the theory of 
stochastic differential equations ([Brissaud and Frisch[ 



1974) a formal solution in Laplace space can be found 
which is hoverer difficult to invert analytically. A nu- 
merical solution is instead feasible and the method can 
be extended to investigate the effect of an ensemble of 
fiuctuators with proper distribution of parameters to gen- 



erate processes having 1// power spectrum (Falci et al. 
2005|). In (jCheng et al.\ |2QQ8 ) a generalized transfer- 



matrix method has been proposed. It reduces to the al- 
gebraic problem of the diagonalization of a 6 x 6 matrix 
whose eigenvalues give the decoherence "rates" entering 
the qubit dynamics. This approach also gives the system 
evolution on time scales shorter than the fluctuator cor- 
relation time. The transfer matrix method is suitable to 



address the dynamics under instantaneous dc-pulses like 
those required in the spin-echo protocol. The method can 
in principle be extended to the case of many fiuctuators, 
though the size of the matrices soon becomes intractable. 
Another approach consists in the evaluation of the quan- 
tum dynamics of a composite system composed of the 
qubit and one (or more) fiuctuator(s), treated as quan- 
tum mechanical two- level systems ( [Faladino et al. , 2003a 



2002 ) . This requires either solving the Heisenberg equa- 



tions of motion, or a master equation in the enlarged 
Hilbert space. Impurities are traced out at the end of 
the calculation, when the high-temperature approxima- 



tion for the fiuctuator is also performed (Faladino et al. 



2003a, 2002) 



Both methods lead to analytic forms in 
a limited parameter regime and a numerical analysis is 
required in the more general case. Other approximate 
methods rely on the extension of the approach based 
on the evaluation of the probability distribution p(^, t) 
which is evaluated numerically (Bergli et al. 2006[ ). 

The main effects of a RT fiuctuator coupled to a qubit 
via an interaction term of the form ( [78| ) can be simply 
illustrated as follows. The noise term —hx{t)cFz/2 in- 
duces two effective splittings, rt± — \/{^z i ^Y + ^^5 
and correspondingly two values of the polar angle, 0± = 
diictdiiiiytx/iS^z ^b))-) as illustrated in Fig. [32| 

In the adiabatic limit 7 ~ |r^+ — VL-\ <C 1^±, and 
neglecting any qubit back-action on the fiuctuator, the 
qubit coherence takes a form similar to the pure de- 
phasing result Eq. (49). Here we report the average 



{(Jy{t)) (Faladino et al.^ 2003a) 



KW) 



-Im 



,i(f^+7^/2)t 



^A(±a)e-^ 



-It 



(79) 



where a = ^l- g'^ - 2ig5p^^ - (1 - c^){l - 5pl^), 

c = cos[(^_ - ^+)/2, g = {n^- ^-)h, 
A{a) = (a + c^ - ig')p^;l{^)po + (a + c^ + ig')p^+l{^)vi 

with g^ = g -\- iSpeq{l — c^) and p\_{^) are values at 
t = of the qubit coherences in the eigenbasis |±)6'± 
of the qubit conditional Hamiltonians corresponding to 
X = ±1. In the high-temperature regime, (^Peq ^ O7 
and a ~ \/c^^^^• Thus the relevant parameter sep- 
arating the weak- from the strong-coupling regime is 
g/c^ = (1]+ — Q-)/{c^j). A fiuctuator is weakly-coupled 
if g/c^ <C 1 and strongly coupled otherwise. This con- 
dition depends on the qubit 's operating point, and at 
^ = we get g/c^ = 26/7. Thus a fiuctuator character- 
ized by a set (6, 7) affects in a qualitatively different way 
the qubit dynamics depending on the qubit 's operating 
point. In particular, a fiuctuator turns from strongly to 
weakly coupled increasing 6 from to 7r/2, as illustrated 
From ( 79 ) it easy to see that a weakly coupled 



in Fig. 33 



fiuctuator induces an exponential decay with T2 given by 



Eq. (16), whereas a strongly coupled fiuctuator induces 
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FIG. 33 The threshold value ^i^v = (Q+ - Q-)/{2bc^) for 
a fluctuator behaving as weakly coupled depends on t he op- 
erating point {v/Q = 2b /Q). Adapted from ( Paladino et al\ 
2003b), with kind permission of Springer Science+Business 



Media. 



a saturation effect which is less effective at 6> = 7r/2 than 
at l9 = (iBergH et al\ 120061 ICheng et al\ 120081 iPal 



adino et al. 



2003a| |Zhou and Joynt| |2010D . The effect 



of strongly coupled fiuctuators is more sensitive to devi- 
ations from the optimal point, = 7r/2. The dependence 
on the working point of the decoherence and relaxation 
times for a weakly and a strongly coupled fluctuator is 
shown in Fig. 34 The decay rates at ^ = 7r/2 have also 



been obtained in (Itakura and Tokura 2003|) by solving 



the stochastic differential equation. 



a. Ensemble of fiuctuators: decoherence due to 1/f noise - 
All the above mentioned approaches become intractable 
analytically as soon as the number of fiuctuators in- 
creases. Already in the presence of few fiuctuators a 
numerical solution is required. This difficulty and the 
physical scenario behind the mechanisms of dephasing 
due to 1/f noise at different operating points, can be 
traced back to the dependence of the qubit splitting on 
the induced fluctuations, cf. with Eq. (36), which we 
write here for Q±: 



6^4 






2 1^3" ■ 



(80) 



Far from the optimal point, the leading term is linear in 
the coupling h. This is the only effect left at pure de- 
phasing, Vtx = 0, and it is also the reason way an exact 
formula is available under this condition both in the case 
of a single fluctuator and in the presence of an ensem- 
ble of fiuctuators. In this last case the phase memory 
functional for any number of fiuctuators is found by sim- 
ply multiplying the phase-memory functionals for differ- 
ent fiuctuators. At the optimal point instead, Vtz = 0, 
the first order term is quadratic. Thus the effect of the 
fluctuator is reduced with respect to the pure dephasing 
regime. The effects of different fiuctuators in this case 
do not simply add up independently and no analytic so- 
lution is available. This simple observation suggests an 
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FIG. 34 (Color online) Relaxation 1/Ti and dephasing I/T2 
rates (rescaled for visibility) as functions of 0. Triangles - 
I/T2* = cos^ 0S(0)/2 for a weakly (a) and strongly (b) coupled 
fluctuator. 7/Q: (a) - 0.5, (b) - 0.1. b/Q: (a)=0.1, (b) - 
0.3. Solid (dashed) lines correspond to a symmetric ( slight ly 
asymmetric) fluctuator. Adapted from (Cheng et al. 2008). 



important physical insight on the combined effect of var- 
ious fiuctuators at different working points. When the 
first order effect on the splitting is quadratic, even if 
the different fiuctuators in themselves are independent, 
their effect on the qubit will be influenced by the posi- 
tion of all the others. At pure dephasing, since effects 
sum up independently, all slow fiuctuators are ineffec- 
tive at times 7^ <C 1. This is no longer true when the 
leading term is quadratic, as they play a role in deter- 
mining the effect of faster fiuctuators. In fact even if 
they do not have time to switch during the experiment, 
they contribute to the average effective operating point 
the qubit is working as seen by the faster fiuctuators. 
Thus very slow fiuctuators may be of great importance 
at the optimal point. For example, in Fig. [35] it is il- 
lustrated that even a single, strongly-coupled fluctuator, 
out of the many fiuctuators forming the 1/f spectrum, 
determines a dephasing more than twice the prediction 
of the weak coupling theory, T2* , even when the fluctua- 
tor is not visible in the spectrum. Further slowing down 
the fluctuator produces a saturation of dephasing, as ex- 
pected based on the above considerations. In this regime 
effects related to the initial preparation of the fluctuator, 
or equivalently effects of the measurement protocol, are 



visible as illustrated in Fig. [36] The presence of selected 
strongly coupled fiuctuators in the ensemble leading to 
1/f spectrum may give rise to striking effects like beat- 
ings between the two frequencies rt± , see Fig.|39]and the 



discussion in Sec. III.B.l.b Relevant effects of 1/f noise 
at general working point, especially for protocols requir- 
ing repeated measurements, are captured by approximate 
approaches based on the adiabatic approximation which 
we discuss in the following Subsection. 
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FIG. 35 The Fourier transform of {(Tz{t)) for a set of weakly 
coupled fluctuators plus a single strongly coupled fluctuator 
(thick line). The separate effect of the coupled fluctuator 
[g — 8.3, thin line) and the set of weakly coupled fluctuators 
(dotted line) is shown for comparison. Inset: corresponding 
power spectra. In all cases noise level at A = Ej is fixed to 
the value S{Ej)/Ej = 3.1 8 x 1Q~^ (from typical 1// noise 



amplitude in charge qubits ([Covington et al.\ 2000[ [Nakamura 



et al. 2002| Zorin et al. |1996) extrapolated at GHz frequen- 
cies) . Adapted from ( Paladino et al. 2002 ) . 



B. Approximate approaches for decoherence due to 1// 
noise 

In the context of quantum computation, the effects 
of stochastic processes with long-time correlations, like 
those characterized by 1// spectral density, depend on 
the quantum operation performed and/or on the mea- 
surement protocol. Here we present approximate ap- 
proaches proposed to predict dephasing due to 1// noise 
and its interplay with quantum noise. An approach based 
on the adiabatic approximation ( Falci et al. , 2005[ |Ithier 



et al. 



2005 ) allows simple explanations of peculiar non- 



exponential decay reported in different experiments with 
various setups. In some protocols a Gaussian approxi- 
mation captures the main effects at least on a short time 



scale (Falci et al. 2005 Makhlin and Shnirman 2004 



[Rabenstein et al. 2004). Some other protocols, instead. 



clearly reveal the non-Gaussian nature of the noise. The 
adiabatic approximation suggests a route to identify op- 
erating conditions where leading order effects of 1// noise 
are eliminated also for complex architectures. 



1. Approaches based on the adiabatic approximation 



Our staring point is the phenomenological Hamilto- 



nian Eq. ( 43 ) where we separated the effect of quantum 



noise due to (high frequency) modes exchanging energy 
with the system and of low frequency fluctuations, 5q{t)^ 
of the bias parameter q. Environments with long-time 
memory, i. e., correlated on a time scale much longer 
than the inverse of the natural system frequencies, be- 
long to the class of adiabatic noise. Stochastic processes 
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FIG. 36 The Fourier transform {az)uj for a set of weakly cou- 
pled fluctuators plus a strongly coupled fluctuator (26/7 = 
61.25) prepared in the ground (dotted line) or in the excited 
state (thick line). Inset: corresponding power spectra (the 
thin l ine corresponds to the extra fluctuator alone). Adapted 
from (iPaladino et a/. 1 120021). 



can be treated in the adiabatic approximation provided 
their contribution to spontaneous decay is negligible, a 
necessary condition being t <C Ti oc S{Q)~^. This con- 
dition is satisfied for 1// noise at short-enough times, 
considering that S{uo) ex l/uj is substantially different 
from zero only at frequencies uo <^ Vt. For pure dephas- 
ing, ^ = 0, relaxation processes are forbidden and the 
adiabatic approximation is exact for any S{uj). In the 
adiabatic approximation the instantaneous Hamiltonian 
of a qubit, manipulated only with dc-pulses, reads 



^tot = ^ (cos Oq Gz + Sm OqCFx)^ G^ 



(81) 



where we remind that E(t) 
instantaneous splitting is 



25q{t){dVLz/dq), and the 



n{q,5q{t)) = J n'^{q) ^ E'^{t) ^ 2E{t)n{q) cos Oq . (82) 



The az axis forms a time-dependent angle, Oq{t) = 
diTctdin[Q sinOq/ {ft cos 6q + E{t))]^ with a^. In the qubit 
eigenbasis {E = 0) the adiabatic Hamiltonian (81) reads 



H 



hQ{q^ Sq) 



cos [Oq{t) 



\crz' 



■S\n[eq{t) -Oq]cF^ 



(83) 

The effect of an adiabatic stochastic field E(t) is to pro- 
duce fluctuations both of the qubit splitting, Vt[q^5q{t)]^ 
and of the qubit eigenstates, or equivalently of the "di- 
rection" of the qubit Hamiltonian, expressed by the an- 
gle Oqit). Since we are interested to situations where 
\E{t)\ <C r^, both Vt[q.,5q{t)] and Oq{t) can be expanded 
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in Taylor series about ^(g, 0) and Oq respectively 

/90 1 /^^O 

n{q,5q) « n{q,0) + ^{Sq) + ~{Sq? + ■ 



(84) 
(85) 



where {5q) = Sq{t) and all derivatives are evaluated at 
Sq = 0. The adiabatic Hamiltonian (81) therefore can be 
cast in the form 

H=^ {n{q,0)cT,^ ^ Sn\\{t)(j,^ ^ 5n^{t)(Ti_) (86) 

where, from Eq.(|83|, Sft\\ includes the derivatives of 
ft[q^ Sq{t)] cos [6 q {t) — 6 q] and drt± the derivatives of 
Q[q^ 5q{t)] sin [0 q{t) — Oq]. The Pauli matrix (j± in the 
case of Eq. (|83| is ax' , in general it can be a combina- 
tion of ax' and ay' . 

The effect on the qubit dynamics of adiabatic trans- 
verse fluctuations weakly depends on time, on the other 
side, longitudinal components are responsible for phase 
errors, which accumulate in time. Thus adiabatic trans- 
verse noise has possibly some effect at very short times, 
but the phase damping channel eventually prevails. The 
relevance of these effects quantitatively depends on the 
amplitude of the noise at low frequencies. It has been 



demonstrated in (Falci et at. 2005) (see Fig. 37) that, for 



the typical noise figures of superconducting devices and 
at least for short times scales relevant for quantum com- 
puting, the leading effect of adiabatic noise is defocusing 
originated by the fluctuating splitting during the repeti- 
tions of the measurement runs. This effect is captured by 
neglecting the transverse terms in Eq. (86), i. e., making 



the longitudinal approximation oi the Hamiltonian (81) 



n ^ az 



hsnu{t) 



-a.' 



(87) 



In addition, consistently with Qq(i) ~ Oq, the spHtting 



fluctuations are further approximated as (Falci et al 



2005 Ithier et al. , 2005 



S%{t) 



on 



i^l) 



id'^n 



(Sqf 



Considering the explicit dependence on 6q of ft{q,Sq)^ 
the first terms of the expansion read, cf. with Eq. (36), 



Q{q, Sq{t)) ^ Q{q) + E{t) cos Oq 



1 {E{t) sin Oq)'^ 

2 n{q) 



. (89) 



For longitudinal noise, ^^ = 0, we recover the exact lin- 
ear dependence on the noise of the instantaneous split- 
ting. For transverse noise instead, Oq = 7r/2, the first 
non- vanishing term of the expansion is quadratic. It is 
common to refer to this regime as the "quadratic cou- 
pling" (Ithier et al. 2005), or "quadratic longitudinal 



coupling" (Makhlin and Shnirman, 2004) condition. 



We remark that Eq. (87) also applies when both com- 



ponents, Qx and 1^2, fluctuate. For instance, this is the 
case of flux qubits where both flux and critical current 
fluctuate with 1// spectrum and the physical fluctuating 
quantity Sq is a function of the magnetic flux or of the 
critical current. The partial derivatives can be expressed 
in terms of noise sensitivities as follows 



dn _ on dVt^ dVt dVtx _ ^z dn^ Qx dVtx 
dq dftz dq dVtx dq ft dq ft dq 



(90) 



where the noise sensitivities dflz/dq^ dflx/dq can be in- 
ferred from spectroscopy measurements as in the experi- 
ment ([Bylander et a/.[|201lD. 



A formal expression for the qubit dynamics in the adi- 



abatic approximation, introduced in (Falci et al., 2005) 



has been discussed in details and extended to more com- 
plex gates in (jPaladino et al. 2009). Here we report 



the results in the adiabatic and longitudinal approxima- 
tions. In this regime populations of the qubit density 
matrix in the eigenbasis do not evolve, whereas the off- 
diagonal elements are obtained by averaging over all the 
realizations of the stochastic process E(t) expressed by 
the path-integral 



.w 



.(0) 



= jnE{s 



\P[E{s) 



-i ,QdsQ.^^{q,6q(s)) 



.(91) 



Here P[E{s)\ contains information both on the stochastic 
processes and on details of the specific protocol. It is 
convenient to split it as follows 

P[E{s)] = F[E{s)]p[E{s)], 

where p[E{s)] is the probability of the realization E{s). 
The filter function F[E{s)] describes the specific oper- 
ation. For most of present day experiments on solid- 
state qubits F[E{s)] = 1. For an open- loop feedback 
protocol, which allows initial control of some collective 
variable of the environment, say Eq = 0, we should put 
F[E{s)] (X S{Eo). The different decay of coherent oscil- 
lations in each protocol in the presence of adiabatic noise 
originates from the specific filter function which needs to 
be specified at this stage. 

A critical issue is the identification of p[E{s)] for the 
specific noise sources, as those displaying 1// power 
spectrum. If we sample the stochastic process at times 
tk = kAt, with At = t/n^ we can identify 

p[E{s)]= lim p,+i(E„t;...;^i,ti;^o,0) (92) 

where Pn-\-i{' - ■) is di n -\- 1 joint probability. In general, 
this is a formidable task. However, a systematic method 
can be found to select only the relevant statistical infor- 
mation on the stochastic process out of the full charac- 



terization included in p[E{s)] (Falci et a/., 2005 Paladino 



et al. 2009) 



The signal decay in FID is obtained by performing in 



Eq. (91) the Static Path Approximation (SPA), which 
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consists in neglecting the time dependence in the path, 
E{s) = Eq and taking F[E] = I. In the SPA the prob- 
lem reduces to ordinary integrations with pi(£^o,0) = 
p{Eo). The qubit coherences can be written as poi{t) = 
Poi(O) exp[— ir^t — i^{t)] with the average phase shift 



$(t) ^iln(fdEop{Eo) e'W^'+^S+^''^ocos9\ 



(93) 



Clearly, Eq. ( 93 ) describes the effect of a distribution of 



stray energy shifts h[Q{q^Sq) — Q{q,0)] and corresponds 
to the rigid lattice breadth contribution to inhomoge- 
neous broadening. In experiments with solid state de- 
vices this approximation describes the measurement pro- 
cedure consisting in signal acquisition and averaging over 
a large number N of repetitions of the protocol, for an 
overall time tm (minutes in actual experiments). Due 
to slow fluctuations of the environment calibration, the 
initial value, ft cos + Eq^ fluctuates during the repeti- 
tions blurring the average signal, independently on the 
measurement being single-shot or not. 

The probability p{Eq) describes the distribution of the 
random variable obtained by sampling the stochastic pro- 
cess E{t) at the initial time of each repetition, i. e., at 
times tk = k tm/N^ /c = 0, . . . , A/" — 1. If Eq results from 
many independent random variables of a multi-mode en- 
vironment, the central limit theorem applies and p{Eq) 
is a Gaussian distribution with standard deviation a, 



{E') 



duo S{uj) . 



where integration limits are intended as 1/t^, and the 
high-frequency cut-off of the 1// spectrum, jm- In the 
SPA the distribution standard deviation, a, is the only 
adiabatic noise characteristic parameter. If the equilib- 



rium average of the stochastic process vanishes, Eq. (93) 
reduces to 



^{t) ^iln 



dEn 



V2 



Jt^/Q?+El+2nE^i 



na^ 



■ (94) 



Expanding the square root in the above expression we 
obtain (Falci et al^ 2005} 



m = 



n{cos Oat)'^ 
n + i sin^ OaH 



In 






(95) 



The short-times decay of coherent oscillations qualita- 
tively depends on the working point. In fact, the sup- 
pression of the signal, exp[Im$(t)], turns from an expo- 
nential behavior, ex e-(cos6'cr^) /'^^ ^it ^ to Si power 
law, [1 -h (sine'^a'^t/n)'^]-^/^, at (9 ?^ 7r/2. In these limits 
Eq. (|95| reproduces the results for Gaussian 1// envi- 
ronments in the so called "quasi-static case" reported in 
(Ithier et al. 2005). In particular, at 6> = we obtain the 



short-times limit, t <C 1/7m, of the exact result of Palma 
et al. (1996), Eq. (pTJ). In fact, for very short times we 



can approximate sin^(cc;t/2)/(cc;t/2)^ ^ 1 inside the inte- 



gral Eq. (11), obtaining the exponential quadratic decay 
law at pure dephasing predicted by the Gaussian approx- 
imation. At ^ = 7r/2 the short and intermediate times 
result of Makhlin and Shnirman (2004) for a Gaussian 



noise and "quadratic coupling" is reproduced. 

The fact that results of a diagrammatic approach with 



a quantum environment, as those of Makhlin and Shnir- 



man 



(2004), can be reproduced and generalized already 



at the simple SPA level makes the semi-classical approach 
quite promising. It shows that, at least for not too long 
times (but surely longer than times of interest for quan- 
tum state processing), the quantum nature of the envi- 
ronment may not be relevant for the class of problems, 
which can be treated in the Born-Oppenheimer approx- 
imation. Notice also that the SPA itself has surely a 
wide validity since it does not require information about 
the dynamics of the noise sources, provided they are 
slow. For Gaussian wide-band 1// noise and for times 
7m <C 1/t < 7m, the contribution of frequencies a; <C 1/t 
can be approximated by Eq. ([95]) where the noise vari- 
ance is evaluated integrating the po wer spectrum from 
7m to 1/t, that is a^ = A\ii{l/^rnt) ([Cottet et a/.l|2001 



Ithier et al\ |2005[ [Nakamura et al\ |2002[). The dia- 



grammatic approach of Makhlin and Shnirman (2004) 
for Gaussian noise and quadratic coupling also predicts 
a crossover from algebraic behavior to exponential decay 
at long times. At ^ 1^ with rate A/ 2. This behavior 
is, however, hardly detectable in experiments, where the 
long-time behavior is ruled by quantum noise (see dis- 
cussion in the following section III.B.l.b) 



Equation (91) can be systematically approximated by 



proper sampling E(t') in [0,t]. In this way also the dy- 
namics of the noise sources during the runs can be sys- 
tematically included. For the first correction, p[E{t')] can 
be approximated by the joint distribution p[£^^t; £^o, 0], 
where Et = E(t). At = 7r/2 for generic Gaussian noise 
we find 



Ht) 



In 
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where U{t) = J^{duj/a^)S{uj) (l 



sn 



is a transition 



probability, depending on the stochastic process. For 
Ornstein-Uhlenbeck processes it reduces to the result of 
Rabenstein et al] ( |2004 ). The first correction suggests 
that the SPA, in principle valid for t < I/7M, may have 
a broader validity (see Fig. 37). For 1// noise due to a 
set of bistable impurities it is valid also for t ^ I/7M, 
if 7m ^ ^- Of course, the adiabatic approximation is 
tenable if t < Ti = 2/S{n). 

As already pointed out for the SF model, the signal 
decay during echo protocols is very sensitive to the dy- 
namics of the noise sources within each pulse sequence. 
The SPA cannot capture these effects and would yield 
no decay for the echoes. On the other side, predictions 
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FIG. 37 (Color online) Simulations of an adiabatic fluctua- 
tors 1// environment at ^ = 7r/2. Relaxation studied via 
{ax) (green line) is well approximated by the weak coupling 
theory, T2(dots). Dephasing in repeated measurement damps 
the oscillations (thin black line). Part of the signal is re- 
covered if the environment is re-calibrated (thin gray line). 
Noise is produced by rid = 250 fluctuators per decade, with 



l/U 



lO^rad/s < 7i < 7m = lO^rad/s < Q = 10^"rad/s. 



The coupling v — 0.02^2 is appropriate to charge devices, and 
corresponds io S — IQnAEc/uj with A — 10~^ (Zorin et al. 
1996). The adiabatic approximation, Eq. (|91|, fully accounts 



for dephasing (red dotted-dash line). The static-path approx- 
imation (SPA), Eq. ( |95[ ), (blue solid line) and the first cor- 
rection (blue dashed line) account for the initial suppression, 
and it is valid also for times t ^ I/7M. In the inset, Ram- 
sey fringes with p arameters appropriate to the experiment 
( Vion et a/.| |2Q02| ) (thin black hues). The SPA (blue sohd 
line), Eq. ( |95| ), Is in excellent agreement with observations, 
and also predicts the correct phase shift of the Ramsey sig- 
nal (blue dots, compared with simulations for small detuning 
(5 = 5 MHz, violet line), which tends to ~ 7r/4 for large times. 
Adapted from (|Falci et al\ [2005 ) 



critically depend on the specific sampling of the stochas- 
tic process E(t')^ making it difficult to obtain reliable 
estimates. The decay of the echo signal is often approxi- 
mated assuming Gaussian quasi- static {t < 1/7m) noise 



(Ithier et a/., 2005). 



Several experiments confirmed that the effect of 1// 
noise in repeated measurements protocols can be de- 
scribed by the above simple theory ([Chiarello et al.\ 2012 
Gottet et all [200T] llthier et alA 120051 IMartinis et al 



20031 [Sank eTaTj [20121 |Van Harlingen et aL\ [20041 [Yan 
et a/., 2012). Of course, the dominant 1// noise sources 
are device and material dependent and the decay of the 
measured signal depends on the measurement procedure. 
An accurate estimate and comparison with Gaussian the- 
ory of defocusing has been done in ( Ithier et al] 2005 ) 
for the quantronium. A clear evidence of an algebraic 
decay at the optimal point due to quadratic 1// charge 
noise has been reported, see Fig. [SS} in agreement with 
the prediction of 



Falci et al. 



(2005j),Fig.[37 



In (Van Harlingen et al. , 2004) the effects of 1// noise 



in the critical current in various superconducting qubits 
incorporating Josephson junctions have been computed 
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FIG. 38 Ramsey signals at the optimal point for (jOro/27v = 
106 MHz and Au = 50 MHz, as a function of the delay At be- 
tween the two 7r/2 pulses. Top and middle panels: solid lines 
are two successive records showing the partial reproducibility 
of the experiment. Dashed lines are a fit of the envelope of 
the oscillations in the middle panel leading to T2 = 300 ns. 
The dotted line shows for comparison an exponential decay 
with the same T2. Bottom panels: zoom windows of the mid- 
dle panel. The dots represent now the experimental points 
whereas the solid line is a fit of the whole oscillation with 
A(jj/27v = 50.8 MHz. Arrows point out a few sudden jumps 
of the phase and amplitude of the oscillation, attributed to 



strongly coupled charged TLFs. Adapted from ( ,Ithier et al. 
|2005}. 



numerically showing that the envelope of the coherent os- 
cillations scales as exp [— t^/(2r?)], in qualitative agree- 
ment with the adiabatic and longitudinal approximation. 
Interestingly, the extrapolated decay time r^ depends 
both on the elapsed time of the experiment and on the 
sequence in which the measurements are taken. Two av- 
eraging methods were employed: a time-delay averaging 
(method A) which averages only high-frequency fluctu- 
ations at each time-delay point, and a time-sweep aver- 
aging (method B) which averages both high- and low- 
frequency components. Method A gives longer dephas- 
ing times than method B, since the number of decades of 
1// noise that affect the qubit dynamics in method A is 
smaller than the number of decades sampled in method 
B. For a large number of repetitions, r^ for method B 
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approaches the prediction of iMartinis et al. ( 2003 ) of the 



effect of critical current 1// noise on a phase qubit in 
the Gaussian approximation, i. e., Eq. (11) with inte- 



gration limits and Vt. The dephasing time r^ is found 
to scale as /q = ^^Sj^ (IHz), where 5'/q(1Hz) is the 
spectral density of the critical-current noise at 1 Hz, and 
A = lodQ/QdIo is computed for given parameters for 
each type of qubit that specifies the sensitivity of the 



level splitting to critical-current fluctuations (Van Har- 



lingen et a/., 2004). 



Of course, any setup also suffers from noise sources ac- 
tive at frequencies around GHz, which cannot be treated 
in the adiabatic approximation. In Sec. HI.B.l.bj we will 
discuss a convenient procedure to deal with the interplay 
of noise sources responsible for spectral fluctuations in 
different frequency ranges. 



a. 1// noise during ac-driven evolution: decay of Rabi oscil- 
lations - The manipulation of superconducting qubits is 
often achieved by varying the control parameters in a res- 
onant way at a microwave angular frequency close to the 
qubit transition frequency, Q. Rabi oscillations are rou- 



tinely measured in different labs, e. g., (Chiorescu et al 



20031 IMartinis et aU |2QQ2| [Nakamura et a/.[ |2001| |Vion| 



etal. |2002| Yu et al] 2002). The decohering effect of 1// 



noise during driven evolution is actually weaker than in 
the undriven case. The intuitive reason is that ac-field 



averages the effects of noise ( Ithier et al , 2005 ) . This 
is explained treating the ac-driven noisy system in the 



adiabatic and longitudinal approximation (Falci et al. 
[2012 ). 

A qubit acted by an external ac-field can be modeled 
by the Hamiltonian (|43]) where the control is operated 
via Hc{t) = hB cos{ujt)Q. The device is nominally bi- 
ased at q and its low- frequency fluctuations 6q{t) can be 
treated in the SPA. The problem is solved in the Rotating 
Wave approximation which keeps only control entries "ef- 
fective" in triggering transitions between different states. 
This effective part of the control also depends on the ran- 
dom variable 5q. The populations in the rotating frame 
are readily found, e. g., the population of the first excited 
state is Pi{t\Sq) = {QR/2Qfi) [1 — cos{Q fit)] where the 
flopping frequency for Rabi oscillations is Qfi{q^5q) = 
^r]^{q,Sq) ^nj^{q,Sq). Here Qniq^Sq) = BQio{q,Sq) is 
the peak Rabi frequency, Qio being the matrix element 
of Q in the instantaneous eigenbasis of Ho{q-\-Sq) and the 
detuning is r]{q^Sq) = Q{q^5q) — uo. Averages are evalu- 



Gaussian distribution of 5q with variance (Tq one obtains 



-i^{t) 



1 



1 + i n'j^iap 



exp 






.(97) 



This equation describes different regimes for the de- 
cay of Rabi oscillations, namely a Gaussian time de- 



cay 



-i^{t)\ 



"|(^/zC^q t) 



when the linear term 
in the expansion dominates, and power-law behavior 
- l/[aq{n'}ity^'^] when Q'^i -^ 0. In this regime Eq. ^ 
describes the initial suppression of the signal. The coef- 
ficients of the expansion depend on several parameters 
such as the amplitude of the control fiels and the nom- 
inal detuning r]Q = ^^(^^,0) — u. Further information, 
such as the dependence on q of the energy spectrum and 
of matrix elements Qij{q)^ comes from the characteriza- 
tion of the device. The variance can be related to the 
integrated power spectrum,, and can be extracted from 



FID or Ramsey experiments (Chiarello et al. 2012 Ithier 



et al. , 2005 ) . A similar power law decay of Rabi oscilla- 
tions has been observed for an electron spin in a QD, 
due to the interaction with a static nuclear-spin bath in 



(Koppens et al. 2007) 



Notice that even if Eq. (97) describes the same regimes 



of the SPA in the undriven case, Eq. (95), here the sit- 



uation is different. In particular, Eq. (97) quantitatively 



accounts for the fact that ac-driving greatly reduces de- 
coherence compared to undriven systems. In particular, 
at resonance, ?7o = 0, non- vanishing linear fluctuations of 
the spectrum, dVt/dq ^ 0, determine quadratic fluctua- 
tions of ^fi{q, Sq) (neglecting fluctuations of Qij). Thus 
Q'^i = and Rabi oscillations undergo power-law decay, 
whereas in the absence of drive they determine the much 
stronger Gaussian decay ^ e~2(^^/^^) ^ ^ of coherent 
oscillations. In this regime measurements of Rabi oscil- 
lations ( Bylander et a/.[|2Ql l) have been used to probe 
the environment of a flux qubit. At symmetry points, 
dVt/dq = 0, coherent oscillations decay with a power law, 
whereas Rabi oscillations are almost unaffected by low- 
frequency noise. For non- vanishing detuning decay laws 
are equal being the system driven or not. 

This picture applies to many physical situations, since 
fluctuations of Qij are usually small, corresponding to a 
fraction of VLr ^ 0, whereas 7^(g, 5q) fluctuates on the 
scale of the Bohr splitting Q. ^ Qji and may be partic- 
ularly relevant for r]o = 0. The dependence Qij{q) may 
have significant consequences in multilevel systems. 

Recently, the influence of external driving on the noise 



ated by expanding Qfi to second order as in Eq. ([84|), ^^^^^^^ ^f Usts^Ue fluctuators was investigated in (IC^ 



stantin et al. 2009 ) . The authors proposed an idea that 



nfl{q,Sq)^nfl{q,0)^n'f,{Sq)^-n'},{Sq)^^--- , (96) 
where ft'fi = dftfi/dq, ft'ji = d'^ftfi/dq'^.Si Assuming a 



external driving may saturate the fluctuators thus de- 
creasing their contribution to the dephasing. A calcu- 
lation based on the Bloch-Redfield formalism has shown 
that the saturation of some fluctuators does not lead to 
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significant decrease in decoherence. Brox et al. (2011) 



took into account the effect of driving on the dynamics 
of the fluctuators. The main result of this analysis is 
the prediction that additional low-frequency driving may 
shift the noise spectrum towards high frequencies. Since 
the dephasing is influenced mostly by the low-frequency 
tail of the noise spectrum this shift decreases decoher- 
ence. However, the predicted effect is not very strong. 



b. Broadband noise: multi-stage approach - In the last 
part of this subsection, we have to warn the reader that 
when comparing the above predictions with experiments, 
one has to keep in mind that in nanodevices noise is typ- 
ically broadband and structured. In other words, the 
noise spectrum extends to several decades, it is non- 
monotonic, sometimes a few resonances are present. We 
already mentioned that typical 1 //-noise measurements 
extend from a few Hz to ~100 Hz (Zorin et a/., 1996). 



The intrinsic high-frequency cut-off of 1// noise is hardly 
detectable. Recently, charge noise up to 10 MHz has 
been detected in a SET ( Kafanov et a/. [[2008 ). Flux and 
critical current noises with 1// spectrum consistent with 
that in the 0.01-100 Hz range were measured at consid- 
erably higher frequencies (0.2-20 MHz) ( [Bylander et al 



2011 



Yan et al. 



2012). Incoherent energy exchanges 



between system and environment, leading to relaxation 
and decoherence, occur at typical operating frequencies 
(about 10 GHz). Indirect measurements of noise spec- 
trum in this frequency range often suggest a "white" or 



Ohmic behavior ( Astafiev et a/.[[2004[|Ithier et a/.[[2005 ). 
In addition, narrow resonances at selected frequencies 
(sometimes resonant with the nano device-relevant energy 



scales) have being observed ([Cooper et al , 2004 Eroms 



et a/.[|2006[|Simmonds et aL\\2004:f . 

The various noise sources responsible for the above 
phenomenology have a qualitative different influence on 
the system evolution. This naturally leads to a classifica- 
tion of the noise sources according to the effects produced 
rather than to their specific nature. Environments with 
long-time memory belong to the class of adiabatic noise 
for which the Born-Oppenheimer approximation is appli- 
cable. This part of the noise spectrum can be classified 
as "adiabatic noise": 1// noise falls in this noise class. 
High-frequency noise is essentially responsible for spon- 
taneous decay and can be classified as "quantum noise" . 
Finally, resonances in the spectrum unveil the presence 
of discrete noise sources which severely effect the system 
performances, in particular reliability of devices. This is 
the case when classical impurities are slow enough to in- 
duce a visible bistable instability in the system intrinsic 
frequency. This part of the noise spectrum can be clas- 
sified as "strongly-coupled noise". Each noise class re- 
quires a specific approximation scheme, which is not ap- 
propriate for the other classes. The overall effect results 
from the interplay of the three classes of noise. In order 



to deal with broadband and structured noise we have to 
resort to a multi-scale theory which can be sketched as 



and Lukin 



2006) 



follows (Falci et a/., 2005 Paladino et al 



2009 Taylor 



Suppose we are interested to a reduced description of a 
n-qubit system, expressed by the reduced density matrix 
p^{t). It is formally obtained by tracing out environ- 
mental degrees of freedom from the total density matrix 
^/'Q,A,sc fj.^^ which depends on quantum (Q), adiabatic 
(A) and strongly coupled (SC) variables. The elimination 
procedure can be conveniently performed by separating 
in the interaction Hamiltonian, which we write analo- 
gously to Eq. (43), as ^- Gz (g) Ei^, various noise classes, 
e. g., by formally writing 

(j^^ ^Ei = (j^^ Ef + (j^^ Ef + (j^^ Ef^ . (98) 

Adiabatic noise is typically correlated on a time scale 
much longer than the inverse of the natural frequen- 
cies f^i, then application of the Born-Oppenheimer ap- 
proximation is equivalent to replace Ef^ with a classi- 
cal stochastic field Ef{t). This approach is valid when 
the contribution of adiabatic noise to spontaneous de- 
cay is negligible, a necessary condition being t <C Tf ex 
S"^{ft)~^. This condition is usually satisfied at short 
enough times since S"^{(ju) ocl/ou. 

This fact suggests a route to trace out different noise 
classes in the appropriate order. The total density matrix 
parametrically depends on the specific realization of the 
slow random drives E{t) = {Ef{t)} and may be writ- 
ten as W^'^'^^(t) = W^^^^[t\E{t)]. The first step is 
to trace out quantum noise. In the simplest cases this 
requires solving a master equation. In a second stage, 
the average over all the realizations of the stochastic 
processes, E{t), is performed. This leads to a reduced 
density matrix for the n-qubit system plus the strongly 
coupled degrees of freedom. These have to be traced out 
in a final stage by solving the Heisenberg equations of 
motion, or by approaches suitable to the specific micro- 
scopic Hamiltonian or interaction. In particular, this is 
the case discussed in the initial part of this Section, of 
the spin-fluctuator model at pure dephasing. The or- 
dered multi-stage elimination procedure can be formally 
written as 



p"(t) = Tr5c V[E{t)]P[E{t)]TTQ 



W^ 



SC 



{t\m)]] 



The multi-stage approach allows making predictions in 
realistic situations when the outcome of a measurement 
results from the effects of various noise classes. For in- 
stance, we can address the interplay of 1// noise with RT 
noise produced by one fluctuator which is more strongly 
coupled with the qubit, having 70 <C l/t <C 1^, but 
bo < Q. Even if the fluctuator is not resonant with the 
qubit it strongly affects the output signal. If ^0 > 1, 
it determines beats in the coherent oscillations and split 
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FIG. 39 (Color online) (a) (ay) at 6* = 7r/2, Q/(27r) = 
10^° Hz. The effect of weak adiabatic 1// noise (light gray 
line) (7 G 27r X [10^ 10^] Hz, uniform 2b = 0.002Q, n^ = 250) 
is strongly enhanced by adding a single slow (7/^ = 0.005) 
more strongly coupled {2bo/Q = 0.2) fluctuator (black line), 
which alone would give rise to beats (red line), (b) When 
the bistable fluctuator is present the Fourier transform of 
the signal may show a split-peak structure. Even if peaks 
are symmetric for the single bistable fluctuator alone (dashed 
line), 1// noise broadens them in a different way (solid line). 
Adapted from (Falci et al. 2005|. 



peaks in spectroscopy, which are signatures of a discrete 
environment. Because of the mechanism discussed in 
Sec. |III.A.2.a| the additional fluctuator makes bistable 



the working point of the qubit and amplifies defocusing 
due to 1// noise. Even if the device is initially optimally 
polarized, during tm the fluctuator may switch it to a 
different working point. The line shape of the signal will 
show two peaks, split by and differently broadened by the 
1// noise in background. The corresponding time evolu- 
tion will show damped beats, this phenomenology being 
entirely due to the non- Gaussian nature of the environ- 
ment. Figure 39, adapted from (Falci et al. 2005), shows 



results of a simulation at the optimal point, where 1// 
noise is adiabatic and weaker than the typical noise level 
in charge qubit s, this picture applies to smaller devices. 
The fact that even a single impurity on a 1// background 
causes a substantial suppression of the signal poses the 
problem of reliability of charge based devices. We remark 
that the reported beatings do not arise from a resonant 
coupling between the qubit and the fluctuator. This sit- 
uation will be address later in this Section. 

Finally, we mention that a commonly used simplified 
version of the multi-stage approach consists in simply fac- 
torizing the effects of different noise classes, in particular 
of adiabatic and quantum noise. For some measurements 
protocols, when the responsible noise classes act on suf- 
ficiently different time scales and the spectra are regular 
around the relevant frequency ranges, this approximation 



leads to reasonable prediction of the signal decay ( Ithier 
et a/.l 120051 iMartinis et a/.l 120031). 



plex algorithms and operations. To this end, scalable 
qubit coupling and control schemes capable of realizing 
gate errors small enough to achieve fault tolerance are 
required. Fluctuations with 1// spectrum represent a 
limiting factor also for the controlled generation of entan- 
gled states (two-qubit gates) and for the reliable preserva- 
tion two- (multi-) qubit quantum correlations (entangle- 
ment memory). In addition to fluctuations experienced 
by each single qubit, solid state coupled qubits, being 
usually built on-chip, may suffer from correlated noise 
sources acing simultaneously on both sub-units. Already 
in 1996, measurements on SET circuits revealed 1// be- 
havior of voltage power spectra on t wo transistor s and of 
the cross-spectrum power density (Zorin et a/., ,1996, ) p] 



In Josephson charge qubits, fluctuator s in the insulating 
substrate might influence several qubits fabricated on the 
same chip, whereas fluctuating traps concentrated inside 
the oxide layer of the tunnel junctions are expected to 
act independently on the two qubits, due to screening by 
the junction electrodes. Background charge fluctuations 
could also lead to significant gate errors and/or decoher- 
ence in semiconductor-based electron spin- qubits through 



inter-qubit exchange coupling (Hu and Das Sarma 2006). 



Recent investigations aimed at identifying operating 
conditions or control schemes allowing protection from 
1// fluctuations in complex architectures. On one side, 
passive protection strategies, like "optimal tuning" of 
nanodevices extending the single-qubit optimal point 



concept (D'Arrigo and Paladino, 2012 Paladino et al 



2011, 2010) and the identification of symmetry pro- 



tected subspaces (Brox et al. 2012 Storcz et 



2005) 



have been proposed. Alternatively, or in combination 



with passive protection, resonant rf-pulses (Rigetti et al. 



2005) and pulse sequences eventually incorporating spin 



echoes ( Kerman and 01iver[ ,2008 j have been considered 
as well. 

In this section we will focus on the first strategies; 
dynamical decoupling approaches will be presented in 
Sec. |III.D[ In addition, we will review the recent the- 
oretical studies about the impact of 1 //-noise correla- 
tions on the entanglement dynamics of coupled qubits. 
We will refer to the strictly related experimental works, 
which have also provided important indications on the 
microscopic origin of the observed noise. Rather ex- 
tensive literature on decoherence of coupled qubits in 
the presence of correlated quantum noise will be omit- 
ted, as well as the variety of relevant experimental works 
demonstrating the feasibility of universal quantum gates 
and simple quantum algorithms based on superconductor 
and semiconductor technologies. It is worth mentioning 
that coupling schemes for superconductor-based qubits 



2. 1/f noise in complex architectures 



Realizing the promise of quantum computation re- 
quires implementing a universal set of quantum gates, 
as they provide the building blocks for encoding com- 



The cross-spectrum of two stochastic processes Xi(t) and X2(t) 
is defined as SxiX2(^) — (l/^r) Jq^ <itCxiX2 (^) coscjt, where 
Cx^x^it) = {Xi(t)X2(0)> - X1X2 and Xo, = {X«(t)>. 
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have been reviewed by Clarke and Wilhelm (2008); since 



the first demonstration of quantum oscillations in su- 



perconducting charge qubits (Pashkin et al. 2003), sev- 
eral benchmarking results have been reached in different 



labs (Di Carlo et al 2010 



Fedorov et 



2012 



Lucero 



et a/.(|2012||Mariantoni et a/.[ |2011| [Neeley et a/.||2010| 



Palacios-Laloy et a/.[ |2010[ |Reed et a/.[|2012D . 

The core element of an entangling two-qubit gate can 
be modeled as 

H= 5] 4"t[fc(t)]+iJcoupling[{9a(t)}], (99) 



where for each qubit H^^^ is given by Eq. (41), and we 



indicate that the interaction term i^coupiing may also de- 
pend on the control parameters {qa}- Following the same 
steps as leading to Eq. (42) and considering only classical 
bias fluctuations one can cast Eq. ([99| as J^ + 3H where 

a ij 

5H = hJ2Y. ^? (i) at + hY, Ei,, (i) al a] , (101) 



i a=l,2 



*,J 



z,j G {x^y^z}. The stochastic processes Ef(t) include 
fluctuations affecting each unit and cross-talk effects due 
to the coupling element: Ef(t) = CaiSqa + daiSq(3 with 
Cai (X d^a/dqa and d^i ex i^ijUqa}) W t^ ^)- Fluctua- 
tions of the interaction energy huij{{qa}) are included in 
Eij{t) = {dh'ij/dqi)6qi + {dvij/dq2)Sq2. Charge fluctua- 
tions affecting the exchange splitting of two electrons in 
a gate-defined double dot ( |Hu and Das Sarma^ |2006| ) or 
background charge- induced fluctuations of the coupling 



capacitance of charge qubits ( Storcz et al. , 2005 ) can be 
modeled by a term of this form. The stochastic processes 
5qo^ might originate from the same source, from different 
sources or a combination. In the case of charge qubits 
for instance, random arrangement of background charges 
in the substrate produce correlated gate-charge fluctu- 
ations to an extent depending on their precise location 



(see Fig. 40 ), whereas impurities within tunnel junction a 



are expected to induce only gate charge fluctuations dq^ 
(jZorin et al.\ |1996|) . These correlations are quantified by 



the intrinsic correlation factor /i, which for stationary and 
zero average processes with same variance, follows from 
{^qa{t)^qi3{t)) = [^a/3 + m(1 — ^a/?)]^^- The Overall degree 
of correlation between the processes Ej {t) and Ef (t) re- 
sults both from intrinsic correlations and from cross-talk 
effects. It is expressed by the correlation coefficient /j^c 
defined, for two generic stochastic processes Xa{t)^ as 



{{X,{t) - X,){X2{t) - X2)) 

V((Xi(i)-Xi)2)((X2(t)-X2)2) 



(102) 



where (...) denotes the ensemble average and X^ 

{Xa(t)). 




Metallic Gate 



FIG. 40 a) Two qubits, Qi and Q2, coupled to a fluctuator 
in the substrate between the two qubits, where the charge 
can tunnel between two sites, b) The qubits are coupled to 
a charged impurity through its image charge on the metallic 
gate, the charge can tunnel between the gate and the im- 
purity. The coupling strengths are given by ^1 and ^2. The 
configuration demonstrated in a) gives origin to anticorrelated 
noise, while con figuration b) give s origin to correlated noise. 
Adapted from ( |Brox et a/.||2012 ). 



The adiabatic approximation scheme introduced in 
paragraph |III.B.l1 can be easily extended to complex ar- 
chitectures to investigate the short-times behavior rele- 
vant for quantum computing purposes. In this approxi- 
mation the effect of stochastic processes with 1// spec- 
trum and/or cross-spectra on universal two-qubit gates 
has been studied in (Brox et al., 12012 D'Arrigo et al. 



20081 ID'Arrigo and Paladino||2012]|Paladino et a/.||2009 



2010) and on entanglement memory element in (Bellomo 



et al., 2010). In this case adiabatic noise does not induce 



the phenomenon of entanglement sudden-death, but it 
may reduce the amount of entanglement initially stored 
faster than quantum noise for noise figures typical for 
charge-phase qubits. An extension of the multi-stage ap- 
proach to complex architectures has been reported in 
(Paladino et al., 2011), where the characteristic time 



scales of entanglement decay in the presence of broad- 
band noise have been derived. Analogously to single- 
qubit gates, low frequency noise induces fluctuations of 
the device eigenergies resulting in a defocused averaged 
signal. One way to reduce inhomogeneous broadening 
effects is to "optimally tune" multi-qubit systems. In 
(Paladino et al. 2010) a general route to reduce in- 



homogeneities due to 1// noise by exploiting tunabily 
of nanodevices has been proposed. The basic idea is 
very simple: in the adiabatic and longitudinal approx- 
imation the system evolution is related to instantaneous 
eigen- frequencies, cj^[E(t)], which depend on the noise 
realization E(t). For H given by Eq. (100), E(t) = 



{£^f (t), Eijit)}. The leading effect of low-frequency fluc- 
tuations in repeated measurements is given within the 
SPA. The frequencies uoimi^) are random variables, with 
standard deviation E^^ = ^/{SoJf^) — {SuJim)'^^ where 
^^im = ^^m(E) — ojim- "Optimal tuning" consists in 
fixing control parameters to values which minimize the 
variance E^^ of the frequencies (^/^(E). This naturally 
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FIG. 41 (Color online) Dispersion in the SWAP subspace 
ScjJ2i (blue) and in the orthogonal subspace dujso (thick gray) 
as a function of transverse fluctuations on qubit 1 xi = 6qi 
(x2 = Sq2 = 0) for resonant qubits with Vzzj^ — 0.01. 
The exact splitting (blue) is compared with a second order 
expansion (dot ted) and the single qu bit dispersion (circles). 
Adapted from ( [Paladino et al. 2010| . 



results in a enhancement of the decay time of the cor- 
responding coherence due to inhomogeneous broadening. 
The short-times decay of the coherence in the SPA is in 
fact given by 



Ke-^'"'"(^^*)i~\/r^Ts^;;^, 



(103) 



resulting in reduced defocusing for minimal variance T^ira. 
For a single-qubit gate, the optimal tuning recipe reduces 
to operating at the optimal point: if Cc;/^(E) is mono- 
tonic, then S^^ ^ Yja\d^imld(la\ cr^^, and the vari- 
ance attains a minimum for vanishing differential disper- 
sion. When bands are non-monotonic in the control pa- 
rameters, minimization of defocusing necessarily requires 
their tuning to values depending on the noise variances. 
For a multi-qubit gate, the optimal choice may be spe- 
cific to the relevant coherence for the considered opera- 
tion. By operating at an optimal coupling, considerable 
improvement of the efficiency of a Vi — SWAP gate re- 
alized via a capacitive coupling of two quantronia has 
been proved, even in the presence of moderate amplitude 

Optimiza- 



charge noise (Paladino et al., 2011 2010) 



tion against 1// flux and critical current noise of an en- 
tangling two-transmon gate has also been demonstrated 



(D'Arrigo and Paladino 2012). Other passive protection 



strategies are based on the use of a "coupler" element 
mediating a controllable interaction between qubits. In 



( [Ker man and Oliver , 2008 ) a qubit mediates the control- 
lable interaction between data qubits. By Monte Carlo 
simulation, the feasibility of a set of universal gate op- 
erations with (9(10~^) error probabilities in the presence 
of experimentally measured levels of 1// flux noise has 
been demonstrated. 

The effect of correlated or partially correlated low- 
frequency noise acting on two qubit gates has been stud- 



iedin (Brox etal, 2012 D'Arrigo et al 2008 Faoro and 



Hekking 2010 Hu et al. , 2007 Storcz et al. , 2005). Due 



the efforts have primarily resulted in numerical survey- 
ing of various situations. The natural question to ask 
is whether correlations between noise sources increase or 
suppress dephasing of the coupled systems compared to 
uncorrelated noise. The answer depends first of all on the 
symmetry of the system-environment interaction, i. e., on 
the existence of decoherence-free subspaces (DFS), possi- 
bly also one-dimensional, and on the initial system state. 
In (Brox et al. 2012), introducing a generalized Bloch- 



sphere method combined with the SPA, analytical ex- 
pressions for the dephasing rates of the two-qubit system 
as a function of the degree of correlation ii^, have been 
derived. For resonant qubits with a a^a^ -coupling in the 
presence of transverse noise, two one-dimensional DFS 
are found, |00) — |11) which does not decay in the pres- 
ence of correlated noise, but which is sensitive to anti- 
correlations (see Fig. 40), and |00) + |11) showing the 
opposite behavior. In the absence of perfect symme- 
try (for instance if qubits are not resonant), the above 
symmetric states are not eigentstates of the Hamiltonian 
and as a consequence are also less sensitive to noise cor- 
relations. In other words, it is both the symmetry of 
the initial state and how much this state overlaps with 
an eigenstate of the Hamiltonian in the absence of noise 
that determines the decoherence rate. This analysis sug- 
gests that for each setup, the most convenient subspace 
for two-qubit encoding should be based on preliminary 
investigation on the nature of noise correlations. For in- 
stance, in the absence of correlations, the SWAP sub- 
space generated by (|01) ± |10))/>/2, is more resilient to 
transverse low-frequency fluctuations with respect to the 
orthogonal subspace, this fact being ultimately due to 
the dependence of the corresponding eigenvalues on the 
deviations ^qa, as illustrated in Fig. |4l] Note that the 



SWAP subspace is expected to be also more stable for 
the single qubit. An analogous conclusion was reached 



in (De et al. 2011) for a pair of qubits coupled via the 
exchange interaction. In (D'Arrigo et a/. [[2008 Hu et al. 



to the complexity of the Hilbert space of coupled qubits. 



2007) a phenomenological model for 1// correlated noise 
affecting a two-qubit gate in a fixed coupling scheme has 
been considered. The effect of noise correlations on en- 
tanglement generation in the SWAP subspace sensitively 
depends on the ratio o jv^z between the amplitude of 
the low frequency noise and the qubits coupling strength 
( D'Arrigo et a/.[[2008 ). For small amplitude noise, cor- 
relations increase dephasing at the relevant short times 
scales (smaller than the dephasing time). On the other 
hand, under strong amplitude noise, an increasing degree 
of correlations between noise sources acting on the two 
qubits always leads to reduced dephasing. The reason for 
this behavior originates from the non-monotonic depen- 
dence of the SWAP splitting variance on the correlation 
coefficient ii^ A numerical analysis has shown that the 
above features hold true for adiabatic 1// noise extend- 
ing up to frequencies 10^ s~^, which are about two orders 
of magnitude smaller that the qubit Bohr frequencies. At 
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longer times, the entanglement decay time (defined as the 
time where the signal is reduced by a factor e~^) weakly 
increases with fic (Hu et al., 2007). 

Recent experiments on flux qubits quantified 1// flux 
noise and flux-noise correlations providing relevant indi- 
cations on its microscopic origin. In (Yoshihara et al. 



relaxation times. The observed specific correlation be- 
tween low-frequency flux noise and inductance fluctua- 
tion suggests that the flux noise is related to the nonequi- 
librium dynamics of the spin system, possibly described 
by spin glass models (Chen and Yu 2010) or fractal 



2010) flux noise correlations have been studied in a 



spin clusters, which appear naturally in a random system 
of spins with wide distribution of spin-spin interactions 



system of coupled qubits sharing parts of their loops, (Kechedzhi et al., 2011). 



whereas in (Gustavsson et al, 2011) a single, two-loop 



qubit was used to investigate flux noise correlations be- 
tween different parts within a single qubit. In both ex- 
periments the qubit dephasing rate was measured at dif- 
ferent bias points. A comparison of the data with the 
rate prediction in the Gaussian approximation based on 
the assumption of l//-type behavior both of spectra and 
cross-spectrum, provided indications of the noise ampli- 
tudes and on the sign of noise correlations. In the first 
experiment it was found that the flux fluctuations origi- 
nating from the shared branch lead to correlations in the 



noise of the two qubits. In (Gustavsson et al, 2011) flux 



fluctuations in the two loops are found to be anticorre- 
lated. Both experiments provided strong indication that 
the dominant contribution to the noise comes from /o- 



ca/ fluctuations, in agreement with (Bialczak et al\ 2007 



Faoro and loSel [20081 [Koch et al\ [2007| [Martinis et al 



2002). In particular, in the setup of Gustavsson et 



(2011) a global fluctuating magnetic field would have 



given positive correlations, which were not observed. Re- 
sults of both experiments are found to be consistent with 
a model where flux noise is generated by local magnetic 
dipoles (randomly oriented unpaired spins) distributed 
on the metal surfaces. A similar conclusion on the local 



origin of flux noise was drawn in (Lanting et al., 2010) 



from measurements of macroscopic resonant tunneling 
(MRT) between the lowest energy states of a pair of mag- 
netically coupled rf-SQUID flux qubits. In this experi- 
ment, the MRT rate peak widths indicate that each qubit 
is coupled to a local environment whose fluctuations are 
uncorrelated with that of the other qubit. Indications of 
magnetic flux noise of local origin in two phase qubits 
separated by 500 /im on the same chip has been recently 



reported in ( Sank et al. , [2012 



In (Sendelbach et aEj 2009[) the cross spectrum of in- 



ductance and flux fluctuations in a dc- SQUID has been 
measured. In this experiment, the imaginary part of 
the SQUID inductance and the quasi- static flux thread- 
ing the SQUID loop were monitored simultaneously as 
a function of time. From the two time series, the nor- 
malized cross spectral density has been computed. The 
inductance and flux fluctuations were found to be highly 
correlated at low temperature, indicating a common un- 
derlying physical mechanism. The high degree of cor- 
relation provided evidence for a small number of domi- 
nant fluctuators. The data were interpreted in terms of 
the reconfiguration of clusters of surface spins, with cor- 
related fluctuations of effective magnetic moments and 



C. Quantum coherent impurities 

The model of two-level tunneling systems formulated 
by [Anderson et a/.[ ( |1972] ) and [Phillips[ ( [1972[ ), illustrated 
in Section [II[ has been extensively tested experimentally 
by ensemble measurements performed on samples having 
a large TLS density, such as structural glasses. Results 
reported in this Section demonstrated that ensembles of 
TLSs, sparsely present in the disordered oxide barrier of 
Josephson junctions or in the insulating substrates, in- 
duce fluctuations with 1// spectrum which are a ma- 
jor source of decoherence in superconducting nanocir- 
cuits. However, the effects unambiguously proving quan- 
tum mechanical behavior of an individual fluctuator in- 
teracting with a qubit were not observed. Only recently 
highly sensitive superconducting circuits could be used as 
"microscopes" for probing spectral, spatial and coupling 
properties of selected TLSs. Understanding the origin of 
these spurious TLSs, their coherent quantum behavior, 
and their connection to 1// noise is important for any 
low- temperature application of Josephson junctions and 
it is a challenge, which will be crucial to the future of 
superconducting quantum devices. 

The first observations indicative of a considerable in- 
teraction of a superconducting circuit with a strongly 
anharmonic quantum system were reported on a large- 
area Josephson junction (^ 10 /im^) phase qubit at Boul- 
der ([Cooper et al.\ [20041 [Simmonds et al.\ [20041). Mi- 



crowave spectroscopy revealed the presence of small unin- 
tended avoided crossings in the transition spectrum sug- 
gestive of the interaction between the device and indi- 
vidual coherent TLSs resonantly coupled with the qubit. 
In these experiments, a small number of spurious res- 
onators with a distribution of splitting size, the largest 
being ~ 25 MHz and an approximate density of one ma- 
jor TLS per ~ 60 MHz, were observed. Magnitude and 
frequency of the TLS considerably changed after thermal 
cycling to room temperature, whereas cycling to 4 K pro- 
duced no apparent effect. Moreover, different qubits in 
the same experimental setup displayed their own unique 
"fingerprint" of TLSs frequencies and splitting strengths. 
Qubit Rabi oscillations driven resonantly with a TLS 
showed considerably reduced visibility with respect to 
off-resonant driving. It has been observed that similar 
spectroscopic observations may also result from macro- 
scopic resonant tunneling in the extremely asymmetric 
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double- well potential of the phase qubits ( Johnson et al. 



2005). The TLS and MRT mechanisms could be dis- 



tinguished measuring the low frequency voltage noise in 
a Josephson junction in the dissipative (running phase) 



regime ( [Martin et al.\\2^^b ). 

Since these first experiments, similar avoided crossings 
in spectroscopy have been observed in different supercon- 
ducting circuits. In phase qubits they have been reported 



in (Bushev et al.\ 2010 Hoskinson et al\ 2009 Lisenfeld 



in ([Deppe et al., 2008 Lupascu et al. 



et a/.||2010a|b||Martinis et a/.||2005r|Nee 



Palomaki et al. , 2010 Shalibo et al. 



2010 



ey e^ a/.||20Q8l 



; in flux qubits 



2009 Plourde 



et al.\ 2005); in a Cooper-Pair-Box (ultrasmall Joseph- 
son junction with nominal area 120 x 120 nm^) - in (Kim 



et al. , 2008); in the quantronium - in (Ithier et al. 2005) 
and in the transmon - in ( Schreier et al. 2008| ). 

The close analogies among these observations, despite 
of the differences in the qubit setups, junctions size and 
materials, confirm that microscopic degrees of freedom 
located in tunnel barrier of Josephson junctions, usu- 
ally made of a 2- to 3-nm-thick layer of disordered ox- 
ide (usually AlOaj, x ^ 1), are at least one common 
cause of these effects. These microscopic degrees of free- 
dom are strongly anharmonic systems and observations 
are fully consistent with coherent TLS behavior. A fur- 
ther confirmation comes from multi-photon spectroscopy 
in phase (Bushev et al., 2010 Lisenfeld et al., 2010b| 



Palomaki et al., 2010 Sun et al. 2010) and flux qubits 
]Lupa§cu et al. 20091) where the hybridized states of the 



combined qubit-TLS systems have been probed under 
strong microwave driving. For instance, an additional 
spectroscopic line in the middle of the qubit-TLS an- 
ticrossing corresponding to a two-photon transition be- 
tween the ground state and the two excitations state 
of the qubit-TLS system has been observed in (|Bushev] 



etal. 2010 Lupascu e^ a/. 2009 Sun e^ a/. , 2010). More- 



over, in some of these experiments, spontaneous changes 
of the resonator's frequency were observed. The insta- 
bility was observed during many hours while the device 



Through the effective, qubit mediated, coupling between 
the TLS and an externally applied resonant electromag- 
netic field "direct" control the quantum state of individ- 



ual TLSs been demonstrated in (Lisenfeld et al., 2010a). 



In this experiment the qubit always remained detuned 
during TLS operations, merely acting as a detector to 
measure its resulting state. A characterization of the TLS 
coherence properties was possible via detection of TLS 
Rabi oscillations, relaxation dynamics, Ramsey fringes 
and spin echo. Measurements at different temperatures 
shown stable TLS resonance frequencies and qubit 's cou- 
pling strenghts. Energy relaxation time is found to de- 
crease quadratically with temperature, whereas the TLSs 
dephasing times had a different behavior, only one of the 
measured TLSs being close to 2Ti. In (Shalibo et al. 



2010) relaxation and dephasing times of a large ensemble 
of TLSs in a small area (~ 1 ijlw?) phase qubit were mea- 
sured (82 different TLSs obtained from 8 different cooling 
cycles of the same sample). Decay times ranged almost 
3 order of magnitudes, from 12 ns to more than 6000 
ns, whereas coherence times varied between 30 ns and 
150 ns. The average Ti followed a power-law dependence 
on the qubit-TLS coupling strength, whereas the average 
dephasing time was maximal for intermediate coupling. 
Authors suggest that both time scales naturally result 
from TLSs dipole phonon radiation and anticorrelated 
dependence of the TLS tunneling amplitude and bias en- 
ergy on low- frequency environmental fluctuations. Non- 
monotonous dependences of the qubit 's decay time on 
the qubit-TLS coupling and temperature were also pre- 



dicted in ( Paladino et al. , 2008 ) for a qubit longitudinally 



coupled to a coherent TLS. In general, different exper- 
iments shown that some TLSs exhibit coherence times 
much longer than those of the superconducting qubits 
dLisenfeld et al.\ |2010a[ [Neeley et al\ |2008, .Palomaki] 



et al.v 



2010) 'n This remarkable fact, together with the 



ability to directly control selected TLSs, shed a new light 
on these microscopic systems. Indeed, it has been pro- 
posed that TLSs in the barrier of a Josephson junction 



was cold dSimmonds et al.\ |2004D, whereas in (|Lupascu| can themselves act as naturally form ed qu bits (|Tian and 



et al. 2009) it was observed in some samples during few 
tens of minutes data acquisition time. Other samples 
were instead stable over the few months duration of the 
experiment. The instability supports the idea that the 
coupled TLSs are of microscopic origin. The qualita- 
tive trend is that small-area qubits show fewer splittings 
than do large area qubits, although larger splittings are 
observed in the smaller junctions ( [Martinis et al. , 2005). 
Time-resolved experiments on phase qubits have 
demonstrated that an individual TLS can be manipu- 
lated using the qubit as a tool to both fully control and 
read out its state. The trajectory (i. e., time record) of 
the switching current of a phase qubit revealed "quantum 
jumps" between macroscopic quantum states of the qubit 
coupled to a TLS in the Josephson tunnel junction, thus 



Jacobs, 2009 Zagoskin et al. 2006). In (Neeley et al. 



2008 ) the first quantum memory operation on a TLS in a 



phase qubit has been demonstrated. An arbitrary quan- 
tum state was transferred to a TLS, stored there for some 



time and then retrieved. In (Sun et al., 2010) creation 



and coherent manipulation of quantum states of a tripar- 
tite system formed by a phase qubit coupled to two TLS 
has been demonstrated. In this experiment, the avoided 
crossing due to the qubit-TLS interaction acted as a tun- 
able quantum beam splitter of wave functions, which was 
used to precisely control the quantum states of the sys- 



providing a way to detect the TLS state (Yu et al. 2008). 



^ TLSs' decay times following from Ramsey fringes are of the order 
of few hundreds nanoseconds, and maximal relaxation times are 
about one microsecond. 
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tern and to demonstrate Landau-Zener-Stiickelberg in- 
terference. Although TLSs were suitable for these initial 
proof-of-principle demonstrations, their use in a quantum 
computer still seems unlikely because of their intrinsically 
random nature and limited coherence times. A relevant 
step further along this direction has been done recently 



natures of this coupling were observed in most of the ex- 
periments which have been consistently explained consid- 
ering purely transverse dipolar interaction ([Bushev et al. 



2010 Lisenfeld e^ a/. 2010b Lupascu e^ a/.l 2009). Other 



by Grabovskij et al. (2012) who report an experiment in 



multilevel spectroscopy experiments did not uniquely pin 
down the coupling mechanism as well. The similar fea- 
tures observed in phase qubits and in the flux qubit ex- 



which the energy of coherent TLSs coupled resonantly 
to a phase qubit is tuned. When varying a static strain 
field in situ and performing microwave spectroscopy of 
the junction, they observed continuously changing en- 
ergies of individual coherent TLSs. Moreover, obtained 
results over 41 individual TLS between 11 and 13.5 GHz 
are explained readily by the tunneling model and, there- 
fore, provide firm evidence of the hypothesis that atomic 
TLSs are the cause of avoided level crossings in the spec- 
tra of Josephson junction qubits. Mechanical strain of- 
fers a handle to control the properties of coherent TLSs, 
which is crucial for gaining knowledge about their phys- 
ical nature. 

Alternative theoretical models of TLSs have been pro- 
posed to explain the avoided level crossings observed in 
qubit spectroscopy data in phase and fiux qubits. It has 
been suggested that the state of the TLS modulates the 
transparency of the junction and therefore its critical 
current, Ic (Constantin and Yu, 2007 Faoro and loffe 



periment (Lupascu et al, 2009) suggest that strongly 



coupled TLS have the same origin in fiux and phase 
qubits, even though the degrees of freedom manipulated 
are different. A charge coupling model is also supported 



by spectroscopic observations in a Cooper-pair box ( Kim 



et a/., 2008 



A distribution of avoided splitting sizes 
consistent with the qubit coupling to charged ions tun- 
neling between random locations in the tunnel junction 
oxide and not directly interacting with each other has 



been reported in (,Palomaki et al. 2010). In (Tian and 



Simmonds , 2007 ) a possible way to resolve the underlying 



coupling mechanism of TLSs to phase qubits through the 
use of a magnetic field applied along the plane of the tun- 
nel barrier inside the junction was proposed. More gen- 
erally, one or two non-interacting qubits may be conve- 
niently used as a probe of a coherent environmen t (,Jeske] 



2006 Ku and Yu, 2005 Simmonds et al.^ 2004). In this 



case two-level defects could be formed by Andreev bound 



states ( Faoro et al. , 2005 de Sousa et al. 2009 ) or Kondo 



et a/.||2012||Oxtoby et flm2009[ [Paladino et a/.]|2008D . 

The controllable interaction between a qubit and a mi- 
croscopic coherent TLS led to a number of interesting 
features also in the qubit time evolution. One aspect 
is the reduced visibility of qubit Rabi oscillations driven 
resonantly with a TLS first observed by [Simmonds et 



ternatively the TLS may couple to the electric field inside 
the junction, which is consistent with the TLS carrying 
a dipole moment located in the aluminum oxide tunnel 
barriers fMartin et 'd] |2005| [Martinis et 'al\ |2005| ). Re- 
cently |A£a£wal_e£aL"( 2013 ) analyzed the interaction with 
phonons of individual electrons tunneling between two 
local minima of the potential well structure due to the 
electron Coulomb interaction with the nearest atoms in 
the insulator. They concluded that the resulting strong 
polaronic effects dramatically change the TLS properties 
providing quantitative understanding of the TLS relax- 
ation and dephasing observed in Josephson junctions. In 
particular, the strain effects observed by [Grabovskij et\ 
(2012) are quantitatively interpreted. 



impurities (Faoro and loffe, 2007 Faoro et al. 2008). Al- (2004). This problem has been investigated using differ- 



ent approaches and under various driving conditions and 



TLS decoherence mechanisms in (Ashhab et al. 



'2005V'Sun et al. [20101 



2006 



Galperin e^ "olj [20051 [Ku and Yu[ [2005[ [Meier and Loss 



The main conclusion of Meier and 



[Loss] ([2005| is that fiuctuators are the dominant source of 
visibility reduction at Rabi frequencies small compared 
to the qubit-TLS coupling strength, while leakage out of 
the qubit computational subspace becomes increasingly 
important for the large Rabi frequencies of experiments 



with phase qubits. In (Galperin et al. 2005 ) the quantum 



Finally, a TLS may modulate the magnetic flux thread- 



ing the superconducting loop ( Bluhm et al. , 2009 Sendel- 



bach et al., 2008). In (Cole et al., 2010) a direct com- 



parison between these models and high precision spec- 
troscopy data on a phase qubit has been performed. Ex- 
perimental data indicate a small or nonexistent longitu- 
dinal qubit-TLS coupling relative to the transverse term. 
In phase and fiux qubits fiuctuations of the critical cur- 
rent or magnetic fiux generate both transverse and longi- 
tudinal components, whereas the coupling to the electric 
field within the junction is purely transverse. Although 
longitudinal coupling cannot be ruled out, no evident sig- 



dynamics of the four-level system subject to an arbitrar- 
ily strong driving ac-field has been investigated includ- 
ing both phase and energy relaxation of the TLS in a 
phenomenological way. It was demonstrated that if the 
fiuctuator is close to resonance with the qubit, the Rabi 
oscillations of the qubit are suppressed at short times and 
demonstrate beatings when damping is weak enough. In 
addition, it was pointed out that if the read-out signal 
depends on the state of the fiuctuator, the visibility of the 
Rabi oscillations can be substantially reduced, a possible 



scenario in (Simmonds et al. 2004). Depending on the 



relative strength of the resonant ac-driving and the qubit- 
TLS coupling, additional features in the qubit dynamics 
have been predicted as anomalous Rabi oscillations and 
two-photon processes involving transitions between the 



four-level states of the coupled qubit-TLS (Ashhab et al. 
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2006 Sun et a/., 2010). Some of these effects have been environment. 



observed in the experiments mentioned above. 

The role of coherent TLSs on qubit relaxation pro- 



cesses has been investigated in (Mueller et al. 2009). In 



this article a qubit is considered as interacting with coher- 
ent TLSs each subject to relaxation and pure dephasing 
processes (in the underdamped regime) and in resonance 
or close to resonance with the qubit. Depending on the 
distribution of the TLSs energies (uniform or with strong 
local fluctuations), the qubit Ti can either be a regu- 
lar function of the qubit splitting or display an irregu- 
lar behavior. Authors suggest that this mechanism may 
explain the smooth Ti-versus-energy curve in large- area 



junction phase qubit s (Cooper et al. 2004 Neeley et al. 



2008 Simmonds et al. , 2004 ) and the seemingly random 



dependence reported in smaller-area phase qubits and in 



flux or charge qubits ( Astafiev et al. , 2004 Ithier et al. 



2005). It is also speculated that the large splittings ob- 
served in spectroscopy of the same phase qubits may re- 
sult from many weakly coupled spectrally dense TLSs. 

A characterization of the effects of bistable coherent 
impurities in superconducting qubits has been proposed 
by Paladino et al. (2008). Introducing an effective im- 



purity description in terms of a tunable spin-boson en- 
vironment, the qubit dynamics has been investigated for 
a longitudinal qubit-TLS interaction. The asymptotic 
time limit is ruled by a dominant rate which depends 
nonmonotonically on the qubit-TLS coupling strength 
and reflects the TLS dissipative processes and temper- 
ature. At the intermediate times relevant for quantum 
computing, different rates and frequencies enter the qubit 
dynamics displaying clear signatures of non- Gaussian be- 
havior of the quantum impurity. 

Finally, the possibility to highlight the coherent inter- 
action between a superconducting circuit and a micro- 
scopic quantum TLS, in principle, allows investigating 
the important question of the applicability domain of the 
classical RTN model. In (iWold et al.\ 12012) the decoher- 



ence of a qubit coupled to either a TLS again coupled 
to an environment, or a classical fluctuator modeled by 
RTN is investigated. A model for the quantum TLS is in- 
troduced where the temperature of its environment, and 
the decoherence rate can be adjusted independently. The 
model has a well-defined classical limit at any tempera- 
ture and this corresponds to the appropriate asymmet- 
ric RT process. The difference in the qubit decoherence 
rates predicted by the two models is found to depend on 
the ratio between the qubit-TLS coupling and the de- 
coherence rate in the pointer basis of the TLS. This is 
then the relevant parameter which determines whether 
the TLS has to be treated quantum mechanically or can 
be replaced by a classical RT process. This result val- 
idates the application of the RT process model for the 
study of decoherence in qubits also when the coupling 
between the qubit and the fluctuator is strong as long 
as the fluctuator couples even more strongly to its own 



D. Dynamical decoupling and 1// noise spectroscopy 

1. Noise protection and dynamical decoupling 

In the last few years several strategies for coherence 
protection have been developed, both for quantum in- 
formation processing and in the broader perspective of 
quantum control. Optimal bias point discussed earlier 
is a passive stabilization (or error avoiding) code very 
successful in solid-state nanodevices. Dynamical Decou- 
pling (DD) relying on repeated application of pulsed 
or switched control is an active stabilization (i. e., er- 
ror correcting) scheme developed in the field of high- 
resolution NMR (Becker, 2000). DD has been proposed 



as a method to extend decoherence times in solid-state 
quantum hardware, and has been recently applied to de- 
couple spin baths in semiconducto r-based qubits ([Barthel 
eTalj [20101 JBluhm et olj [2QTT| |de Lange et al.\ |2012 ' 



2012). 



Ryan et al. [2010) and 1// noise in superconducting 



nanocircuits ([Bylander et al., 2011 Gustavsson et al. 



Coherent averaging of unwanted couplings is at the 
heart of DD. The principle is illustrated by the spin 
echo which is operated by a single 7r-pulse inducing a 
spin-flip transition. Shining a pulse X^^ (evolution op- 
erator U = ax) at half of the evolution time t, say 
{t/2,X7^,t/2}, dynamically suppresses terms ex dy^az 
in the qubit Hamiltonian. In NMR samples unwanted 
terms Hi oc SB a^ are due to static randomly distributed 
local fields. The Bloch vector dynamics for the en- 
semble of spins is defocused resulting in inhomogeneous 
broadening. The Hahn echo {X^/2, t/2, X^, t/ 2, X^/2} 
is routinely used to achieve efficient refocusing ( |Becker[ 
2000 ) . Echoes also switch off "dynamically" two-qubit J- 



couplings of the nuclear spin Hamiltonian in liquid NMR 



quantum computers (Vandersypen and Chuang 2005). 



2. Pulsed control 

Coupling to a stochastic field E{t) induces diffusion in 
the free spin precession and decoherence, which mitigates 
Hahn echoes. Carr and Purcell| ( |1954| ) (CP) recognized 
that sequences of 7r-pulses may suppress spin diffusion 
since they coherently average out E{t)^ Fig. 
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Compos- 
ite pulses are also used in NMR to stabilize a given quan- 
tum gate against errors in the control ( [Becker , 2000 Van- 
dersypen and Chuang, 2005). NMR pulse sequences be- 



yond Hahn echo have been employed in a superconduct- 



ing qubit to demonstrate both gate stabilization ( Collin 
et a/.[[2004 ) and dynamical reduction of decoherence due 
to 1// noise Qlthier et a/^|2005D. 
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From Echo to DD - DD may selectively remove a noisy 
environment with a finite correlation time Tc^ Control 
is operated via a time-dependent Hc{t) describing a se- 
quence of TT-pulses. We consider hard X^^ pulses, whose 
duration is very short tp -^ 0, at times tj (i = j, . . . , N). 
The time-evolution operator reads 

UN{t, 0) = (7^/7(t, tN-i)(Tx . . . U{t2M)(yx U{tuO) (104) 

where Uit^t') describes the noisy free evolutions between 
pulses. Notice that pulses in (JxU{tn-i^tn)crx reverse 
the sign of operators ay^z appearing in U. Therefore 
such "orthogonal" components flip in sequential steps 
of the protocol. It is convenient to use the language of 
average Hamiltonian theory ( Vandersypen an d Chuan"g{ 
[2005) and introduce the effective Hamiltonian iJAr(t), 
defined as e^^^(^)V^ := UN{t,0). Then a periodic 
train of pulses tj+i — tj = At = t/N, implementing a 
sequence named Periodic DD (PDD), with elementary 
block {At,X7^, At,X7^}, tends to average out orthogo- 
nal spin components (provided N is even) ([Viola et at. 



CPMG sequence 

^rt/2 ^n "^71 Yn Yk Yti V;t Kt Yk ^i Yji ^ull 

■ I I I I I I I I I I ■■ . . 

A/n-pulses 
CP sequence 

X'^/2 -'d -Ki '^li '^ '^^ -Ki ^ -'^n \ '^ '^n/l 

_lJ I I I I I I I I Li- 



UDD sequence 

■ I I I L 



J I \ I LLj 



Yn^nll 



FIG. 42 Timing of the CPMG, CP, and UDD pulse sequences 
for N = 10. Adapted by permission from Macmillan Publish- 
ers Ltd.: ([Bylander et a/.[|2Qllt, copyright (2011). 



Hr describes the environment alone. The structure of 
Eq. (105) is still valid, T^it) depending only on the dy- 
namics ruled by Hr. In particular, [Viola and Lloyd] 



( 1998 ) studied an environment of linearly coupled quan- 
tum oscillators, E = ^^ga{o![x + <^a)- They found that 
Eq. (106) holds true, S{uo) being related to the sym- 

. . metrized correlation function of E^ uniquely expressed 

[1999). T his emerges from the perturbative (Magnus) ex- via the spectral density J{uj) = J]^ gl d{uj-uja)^ namely 



pansion ( Vandersypen and Chuang , 2005| ) of Hn , which 
washes these terms out in the limit of continuous flipping. 
At — > 0. In the simple case of a qubit coupled to pure 
dephasing classical noise, H = |[^ + E{t)]az + Hc{t)^ 



S{Lo) = -{E{t)E{Q)- 



huj 



^(0)^(t))„ =coth(^^^) JH. 



calculations can be carried out exactly (Biercuk et al. 
[2OIIJ . Coherences decay as 

Poi(t) = Poi(0)(e-'*"(*)) = poi(0)e-r«W-'^"(*) (105) 

where the phase <I>Ar(t) = J^ dsyN{s)E{s) is obtained 
by sampling a realization of noise with the piecewise 
constant ^Ar(t) whose discontinuities reflect the effect of 
pulses at ti. The decay function T^it) obtained by noise 
averaging depends on the pulse sequence. For Gaussian 
noise with power spectrum S{uj) the averaging yields 



The ultraviolet (UV) cutoff of J{oj)^ ujc^ sets the time 
scale of fastest response of the environment. As for 
classical noise, decoherence is washed out completely for 
At ^ and greatly suppressed for uJcAt ^ 1. 

In general, open loop schemes with a finite set of pulsed 



fields allow to perform fault-tolerant control |Viola et al. 
1999), i. e., to design the dynamics of a quantum sys- 



TNit) 



duj 



S{uj)FN{ujt) 



(106) 



where the filter function FN{ujt) = \yN{^t)\'^ has been 



defined as ( ,Uhrig[|2007D : 



tem to attain a given objective. The simplest goal is the 
effective decoupling of the environment. With respect 
to other active stabilization strategies, as quantum error 
correction or closed loop (quantum feedback) schemes, 
DD has the advantage that only unitary control of a small 
and well characterized system is needed and it does not 
require additional measurement resources. Relying on 
coherent averaging, DD can suppress errors regardless 
of their amplitude. In the last few years optimization 



N 



of pulse sequences (iBiercuk et a/., 2011) has been an ac 



ioutn 1 2 



FN{uJt) = |1 + (-1)^+^6^^^ + 2 ^(-l)^'e^' 

In the absence of pulses the function Fq{uj) = 
4sin^(co't/2) reproduces the decay in a FID protocol, 
whereas in the presence of pulses it yields diffraction pat- 
terns induced by interference in the time domain (Ajoy 



et a/., 2011) and, in particular, to coherent suppression 
of F/v at low frequencies. As a consequence, T^it) de- 
creases and signal decay due to decoherence is effectively 
recovered. 



tive subject of investigation allowing substantial improve- 
ment when dealing with real open quantum systems. 



Optimized sequences and robust DD - Performances of se- 
quences strongly depend on their details such as the par- 
ity of the number of pulses or their symmetrization. For 
instance, in odd N PDD the noise during the final At 
remains uncompensated. Proper symmetrization of the 
sequences, such as CP sequence (Fig. 42) may lead to 
higher order cancellations in H]s[{t). [Meiboom and GiU 



Viola and Lloyd (1998) applied such techniques to se- 



lectively decouple a pure dephasing quantum environ- 
ment, obtained by letting ^azE{t) -^ ^azE-\- Hr, where 



(1958) proposed a refinement (CPMG sequence), which 
is usually very efficient against spin diffusion ([Becker 



2000 ) , since it also averages errors due to control field in- 



50 



homogeneities. Indeed, if pulses are implemented by res- 
onant ac fields, the component producing spin-flip fluc- 
tuates in amplitude due to the same noise responsible for 
the spin diffusion. For a given initial state of the Bloch 
vector, CP accumulates errors in X^^ at second order. 



while in CPMG Y^ errors appear at fourth order ( Borne- 



man et a/., 2010). 



Recently Uhrig (2007) found that non-equidistant 



pulses improve performances in a pure-dephasing spin- 
boson environment. In the Uhrig DD (UDD) sequence 
tj/t = sin^[7rj/(2A^ + 2)] times are such that the first N 
derivatives of the filter function vanish, \d^FN/dz^^ -n ~ 
0, j G {1, 2, . . . A/'}. This ensures that pure dephasing 
is suppressed to 0{t^) in the series expansion of E{t). 
Lee et aL\ ( |2008 ) conjectured that UDD is universal for 
generic pure dephasing model, and [Yang and Liu" (2008) 
proved that generalized UDD suppresses both pure de- 
phasing and relaxation to 0{t^). In the spirit of UDD, 
several new pulse sequences were introduced in the last 
few years, achieving optimization for a given sequence 



duration (Biercuk et al, 2009), or being nearly optimal 



for generic single-qubit decoherence (West et al. 2010a) 



or for specific environments ( |Pasini and Uhrig , 2010). 
Concatenated DD (CDD) proposed by Khodjasteh and 



Lidar (2005) is an alternative scheme based on the idea 



of recursively defined sequences, which guarantee to re- 
duce decoherence below a pulse noise level. Within this 
framework high fidelity quantum gates have been demon- 
strated numerically ( ,West et al. , 2010b). 

The quest for robust DD arises from the general prob- 
lem of the tradeoff between the control resources involved 
and efficient suppression of decoherence. Ideal DD re- 
quires available couplings allowing the synthesis of con- 
trolled evolution (iViola et al, 19991), large pulse repe- 



tition rate, and pulse hardness. Optimization can be 
used together with realistic bounded amplitude control or 



continuous always-on fields schemes (Jones et al., 2012 



Khodjasteh et al. 2011 Viola and Knill 2003) to allow a 



flexible use of resources needed to attain a given decou- 
pling error V^it). 



3. DD of 1// noise 

DD has a large potential impact in solid state coherent 
nanosystems where noise has large low-frequency compo- 
nents. Indeed, DD of Gaussian 1// noise, does not always 
require ultrafast pulse rates ( Shiokawa and Lidar , 2004 ) . 



However, sources responsible for 1// noise are often dis- 
crete producing non- Gaussian noise. For a proper treat- 



ment Eq. ( 106 ) must be generalized accordingly. Initially 



addressed for understanding charge noise in supercon- 



ducting qubits (Falci et a/., 2004 Faoro and Viola , 2004), 



DD in non-Gaussian environments is important in other 
implementations since, independently on the microscopic 
origin, critical current noise and flux noise may also be 



described as due to a collection of discrete sources, and 
also resulting in a 1//" spectrum. Recently this topic 
has attracted large interest also for solid-state quantum 
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hardware based on electron and nuclear spins ( Lee et al. 
Witzel and Das Sarma, 2007a|b). Besides deter- 



mining the decay of coherences, non-Gaussian noises are 
responsible for additional structure (splitting of spectro- 
scopic peaks and beats) observable in the qubit dynamics. 
This deteriorates the fidelity and must be washed out by 
stabilization. 



a. RT noise and quantum impurities - The simplest phys- 
ical non-Gaussian environment is a single impurity cou- 
pled to the qubit. Models of quantum impurities were 



studied in ^Falci et a/.||2004||Lutchyn et a/.[[2008[|Reben 



trost et al. , 2009 ) , whereas the classical counterpart was 



addressed in ( Bergh and Faoro 2007 Cheng et al. , 2008 



Faoro and Viola 2004 Gutmann et al. 2005). Falci 



et al. (2004) modeled the environment by an electron 



tunneling with switching rate 7 from an impurity level 
to an electronic band (Paladino et al. , [2002 ). The pa- 
rameter quantifying Gaussianity is ^ = (0+ — Q-)/j 
(Sec. |III.A.2 ). The problem is tackled by studying the 
reduced dynamics of the Qubit plus Impurity system 
(QI) for arbitrary qubit bias. The band acts a Marko 
vian environment for the QI reduced dynamics ( Paladino 



et al.\ |2003a) and can be treated exactly by a Mas- 



ter Equation. The key point is that DD has no effect 
on a Markovian environment. Thus the quantum map 
in the presence of a number N of pulses can be writ- 
ten as pQi{t) = SNlPQim = {Vexp[CAt]}^pQi{0). 
Here the superoperators C and V describe, respectively, 
the QI reduced dynamics and Xj^ pulses on the qubit, 
1/ (g) ax. The reduced qubit dynamics is obtained by 
tracing out the impurity PQ{t) = Tr/[pQ/(t)] at the end 
of the whole protocol. In this way non- Gaussianity and 
non-Markovianity of the impurity are accounted for ex- 
actly. At pure dephasing {Qx = 0) a simple analytic 
form can be found. The different physics due to a weakly 
(^ < 1) or strongly coupled g > I impurity is discussed 
in Sec. |III.A.l.a] This difference is washed out for large 
flipping rates {N ^ jt) where the environment becomes 
effectively Gaussian with universal behavior TN{t) ^ g^ 



(Fig. 43). On the other hand, for N < jt a crossover is 



clearly observed between different domains of g. Notice 
that in the intermediate regime, A^ < 7^, DD is still able 
to cancel fast noise, ^ < 1, and all the features of the 
qubit dynamics appearing when g ^ 1. A qualitatively 
similar behavior is found also for ^^^^ 7^ where the solu- 
tion requires the diagonalization of En- The new feature 
in this regime is that DD of slow fluctuators, g > 1, 
may be nonmonotonic with the flipping rate, yielding for 
N < jt decoherence acceleration which is reminiscent of 
the anti-Zeno effect. 

Notice that this model reduces to a classical RT flue- 
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FIG. 43 (Color online) Scaled decay rate of the qubit coher- 
ence, VN{t)/g^^ at fixed t — 107~^ and DD (here N enu- 
merates echo pair of pulses). N — ^ corresponds to FID. A 
Gaussian environment with the same power spectrum would 
give, for arbitrary ^, the curve labeled with g=0.1, since 
VN{t) oc g^ . Inset: VN{t) for ^ = 1.1 for different inter- 
vals between pulses At (lines with dots. At = 5, 2, 0.2) are 
compared with the FID ro(t) (thick line) and with results 
obtained by numerical so lution of the stoch astic Schrodinger 
equation. Adapted from (Falci et al. 2004). 



tuator if mutual QI back-action, described by frequency 
shifts, is dropped out. Numerical simulations in this limit 



by Gutmann et al. ( 2005 ) confirmed that decoherence is 



suppressed for large pulse rates N ^ jt; Gutmann et c 
( 2004|) addressed imperfect DD pulses and Bergli and 
Faoro (2007) also found analytic solutions and showed 



that a train of Y^ pulses avoids decoherence acceleration. 
A quantum impurity modeled by a "rotating wave" spin- 



boson model was recently studied by Rebentrost et al. 
( 2009| ) with the numerical Gradient Ascent Pulse Engi- 
neering (GRAPE) algorithm. Authors found decoher- 
ence acceleration at t/N ^ l]/27r next to the optimal 
point, and optimal pulses allowing for relaxation-limited 
gates at larger pulse rates. 



b. 1// noise - An environment composed by set of inde- 
pendent impurities can be treated along the same lines. 
At pure dephasing TN{t) is the sum of independent sin- 
gle impurity contributions, and the analytic solution can 



be found for arbitrary distribution of parameters (Falci 
[et al.\ |2004[ ). An analytic expression valid in the clas- 
sical limit was also found by Faoro and Viola (2004) 
who pointed out that relatively slow DD control rates 
(At ~ 1/7m which is only a soft cutoff of the environ- 
ment) suffice for a drastic improvement. For figures of 
noise typical for experiments, pulse rates yielding recov- 
ery are insensitive to the average coupling strength of 



the impurities ( Falci et a/.[|2004, ). The situation changes 
when the distribution includes individual more strongly 



than PDD ( [Faoro and Viola[ |2004| ). While no decoher- 
ence acceleration is found at pure dephasing, this may 
happen for low pulse rates {QAt ^10) when fi^. ^ and 
noise acquires a transverse part ( Faoro and Viola 2004 ) . 



c. Robust DD - Suppression of 1// pure dephasing lon- 
gitudinal flux noise was demonstrated in a recent exper- 



iment (Bylander et al. 2011) using a 200 pulses CPMG 



sequence yielding a 50-fold improved T2 over the baseline 
value, whereas the performance of UDD sequences was 
slightly worse. Earlier work with few-pulse sequences has 
demonstrated partial suppression of low- frequency trans- 



verse charge noise (Ithier et al. 2005). 



Referring to superconducting qubits, [Cywinski et al. 
( [2QQ8| studied CMPG, UDD and CDD for 1//^ (0.5 < 
Of < 1.5) Gaussian classical noise at pure dephasing with 
UV-cutoff ujc. For pulse rates larger than this cutoff, 
CPMG is the most effective sequences increasing T2, 
UDD keeping however higher fidelity. Instead, CDD does 
not give relevant improvement even if it outperforms 



PDD for a wide range of parameters (Khodjasteh and 



Lidar 2007). For pulse rates smaller than c^c, CPMG 



slightly outperforms all other sequences. CPMG is also 
a better approach for strongly coupled RT noise, non- 
Gaussian features being suppressed in the large pulse rate 
limit. Similar conclusions were drawn by |Lutchyn et al. 
( 2008 ) for a quantum impurity environment of Andreev 
fiuctuators. Pasini and Uhrig (2010) studied sequences 



optimized for specific power-law noise spectra and found 
that they approach CPMG for soft UV-cutoff (1// noise) 
whereas for hard UV behavior (Ohmic) UDD is the lim- 
iting solution. 

Bounded amplitude "dynamical control by modula- 



tion" was proposed by Gordon et al.\ ( 2008 ) who stud 



ied optimization for Lorentzian and 1// pure dephasing 
noise. A practical limitation of this optimal chirped mod- 
ulation is the sensitivity to the low-frequency cutoff. De- 
sign of GRAPE-optimized quantum gates in the presence 
of 1// noise and inhomogeneous dephasing was recently 
investigated by Gorman et al. (2012|). 

Notice that in general solid state nanodevices suffer 
from different noise sources with multi-axis structure of 
couplings. In these cases CDD (Khodjasteh and Li- 



dar , 2005 ) or concatenated UDD sequences ( Uhrig 2009 



West et al. 2010a) may give substantial advantages. 



4. Spectroscopy 



The possibility that DD could be used as a spectro- 



scopic tool was raised in a number of early works ( Falci 



et al. 2004 Faoro and Viola 2004 ) and has been formal 



ized using the concept of filter function (Biercuk et al. 



coupled impurities ( Galperin et a/. 2003 Paladino et a/. , 2011 Uhrig 2007). The key observation is that for a 



2003b). Time symmetric CP is found to perform better Gaussian process, the filter function in Eq. (106) can be 
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interpreted, at a fixed tirae t, as a linear filter ( Biercukl 



et al] |2Q11), transforming the input phase noise E{t) 



to the output phase ^Ar(F), yielding the decay function 
rAr(f) after noise averaging. Each implementation of 
time-dependent control samples the noise in a distinctive 
way determining the form of the filter FN{ujt). For suit- 
able sequences, we define a filter frequency upi such that 
F(ujFit) r^ 1, which roughly corresponds to the minimal 
inter-pulse At. DD is described by a filter with negative 
gain for uj < ujfi^ the steepness of the attenuation yield- 
ing a measure of the effectiveness of the given sequence. 
This analysis allows to develop a filter-function-guided 
pulse design suited to a particular noise spectrum. 

Application to spectroscopy emerges from the observa- 
tion that, in addition to the decoupling regime, there ex- 
ist spectral regions of positive gain about uj = uop^ where 
the effect of the corresponding spectral components of 
noise is amplified. This is apparent from Fig. |44] where 
the modified filter function F{ujt)/uj'^ is plotted, indicat- 
ing the dominant spectral contributions to decoherence. 
For CPMG and FDD there is a single dominant peak. 
The shift of the peak towards larger cu for increasing TV 
indicates that effects of sub-Ohmic noise are reduced by 



DD. Cywinski et al. (2008) proposed that UDD spec- 



troscopy may also give informations on higher moments 
of noise via the additional structure of the filter function 
(Fig. p. 



In ( Yuge et al. 2011 ) a method for obtaining the noise 



spectrum from experimental data has been proposed. 
The method, valid for the case of pure dephasing, is 
based on the relationship between the spectrum and a 
generalized dephasing time, evaluated from the asymp- 
totic exponential decay in the presence of a sufficiently 
large number of 7r-pulses. 



Recently [Bylander et al. (2011) have exploited the 



narrow-band filtering properties of CPMG to measure 
1// flux noise in a persistent current qubit, where flux 
noise is the main source of dephasing away from the opti- 
mal point, Qz 7^ ^' Indeed, Eq. ( |106[ ) is approximated as 
T]sf{t) (X Auj {dft/dq)'^Sq{ujF)FN{(^Ft)/ujp, where ujf = 
7r/2 At is the peak frequency, Aoo is the bandwidth of the 
filter, Sq and dVt/dq are related to the power spectrum of 
fiux noise and to the qubit sensitivity, depending on the 
fiux bias, Eq. ( [QQ] ). This method has allowed to access the 
unexplored spectral region 0.2-20 MHz where 1//^"^ noise 
was detected. Decay of Rabi oscillations provides an al- 
ternative tool for environment spectroscopy in the same 
frequency range (Bylander et al. 2011). As explained in 
Sec. [III.B.l.a] quasistatic noise is efficiently averaged out 
and the observed decay of Rabi oscillations is essentially 
exponential with a contribution from frequency compo- 
nents around Q.r behaving as Fi^ oc [^z{(l)/^]^ SqiVtR) 
Ithier et al. 2005| ). Extracting this 



(Geva et al.^ 1995 



contribution yields an independent information on the 
power spectrum at u; = VLr ^ MHz. Remarkably, the 
very same power laws have been measured at much lower 
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FIG. 44 (Color online) (a) - (c) Modified filter function 
Fn {(jot) / u'^ , for the indicated pulse sequences and N values. 
Dominant spectral contributions to the measured rAr(t) ap- 
pear as peaks in the modified filter functions, (a) Modified 
filter function for FID showing large weight for low-frequency 
noise on a semilog scale, with arbitrary units, (b) Demonstra- 
tion of even-odd parity through the modified filter func tion o f 
FDD on a semilog scale. Adapted from ( [Biercuk et a/.[|2011 ). 



frequencies (0.01 — 100 Hz) by a direct method using the 
noise sensitivity of a free- induct ion Ramsey interference 
experiment ( |Yan et al] |2012 ). The peculiarity of this 
technique is that it enables measurements of noise spec- 
tra up to frequencies limited only by achievable repetition 
rates of the measurements. 



Two-qubits - The extension of DD to entanglement pro- 
tection from 1// noise is a relevant issue, currently un- 
der investigation both theoretically and experimentally. 
The first experimental demonstration of DD protection of 
pseudo two-qubit entangled states of an electron-nucleus 
ensemble in a solid-state environment has been reported 
in [Wang et al.\ ( |201l| ). DD control pulses operated on 
the electron spin suppressed inhomogeneous dephasing 
due to the static Overhauser field induced by the hyper- 
fine interaction between electron spin and surrounding 
nuclear spins in a F:Si material. Recently, in Gustavsson 



et al. (2012) single-qubit refocusing techniques have been 



extended to enhance the lifetime of an entangled state of 
a superconducting flux qubit coupled to a coherent TLS. 
Fluctuations of the qubit splitting due to 1// flux noise 
induce low frequency fluctuations of the qubit-TLS ef- 
fective interaction. Authors demonstrated that rapidly 



53 



changing the qubit's transition frequency relative to TLS, 
a refocusing pulse on the qubit improved the coherence 
time of the entangled state. Further enhancement was 
demonstrated when applying multiple refocusing pulses. 
These results highlight the potential of DD techniques for 
improving two-qubit gate fidelities, an essential prereq- 
uisite for implementing fault-tolerant quantum comput- 
ing. Quantum optimal control theory represents an alter- 
native possibility to design high-fidelity quantum gates. 
[Montangero et a/T] ( 2QQ7| ), using the GRAPE numerical 
algorithm, demonstrated a stabilized two charge-qubits 
gate robust also to 1// noise. For realistic noise figures, 
errors of 10~^ — 10""^, crossing the fault tolerance thresh- 
old, have been reached. A high-fidelity a/SWAP has been 
studied in ( Gorman et a/. [[201 2 ) using GRAPE optimiza- 
tion in the presence of 1// noise. 



IV. CONCLUSIONS AND PERSPECTIVES 



tionalities, such as quantum gates, in a scalable architec- 
ture is a non-trivial open problem. Application of opti- 
mized pulse sequences to non-Mar kovian noise is subject 
of future investigation. We reviewed the current status 
of the ongoing research along this direction, which rep- 
resents an area for future development of the field. 

Despite of the tremendous progress in this field, there 
is still a long way to go until a practically important 
quantum computer will be realized. Many details have 
to be worked out, and at present time it is actually not 
clear which physical implementation of quantum devices 
and even which architecture will be the most advanta- 
geous. However, it is fully clear that ongoing research on 
design and studies of devices for quantum information 
processing will significantly improve our understanding 
of the quantum world and of interplay between classical 
and quantum physics. It will certainly lead to signifi- 
cant development of modern mesoscopic physics and, in 
particular, of quantum electronics. 



In this review we discussed the current state of theo- 
retical work on 1// noise in nanodevices with emphasis 
on implications for solid state quantum information pro- 
cessing. Our focus was on superconducting systems and 
we referred to implementations based on semiconductors 
only when physical analogies and/or formal similarities 
were envisaged. According to the existing literature, rel- 
evant mechanisms responsible for 1// noise in supercon- 
ducting nanocircuits have been largely identified. How- 
ever, in many solid state nanodevices this problem cannot 
be considered as totally settled and details of the interac- 
tion mechanisms remain controversial (see Section [ll|). In 
some cases, available experiments do not allow drawing 
solid conclusions and further investigation is needed. 

We have discussed the role of low-frequency noises in 
decoherence of quantum bits and gates. Various methods 
to address this problem have been presented. A relevant 
issue in connection with quantum computation in the 
solid state is decoherence control and the achievement of 
the high fidelities needed for the succesfull application 
of error correction codes. Various proposals have been 
put forward to limit the effects of 1/f noise. . Since 
such noises in solid-state devices are created by material- 
inherent sources, an obvious way to improve performance 
of quantum devices is optimizing materials used for their 
fabrication. In particular, it is very important to engi- 
neer "dielectric" part of devices. Enormous effort in this 
direction based on identification of the noise sources and 
properties has resulted in significant optimizing of exist- 
ing devices and suggesting novel ones. The main focus 
of the review is relating the device performance along 
different protocols to the properties of the noise sources. 
We hope that understanding these relations may lead to 
improvement of the quantum devices. 

The integration of control tools, like dynamical decou- 
pling sequences appropriate to 1/f noise, with other func- 
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Appendix A: ABBREVIATIONS 

BCS - Bardeen, Cooper, Schrieffer 

ODD - Concatenated Dynamical Decoupling 

CPB - Cooper Pair Box 

CP - Carr Purcell 

CPMG - Carr Purcell Meilboom Gill 

CTRW - Continuous Time Random Walk 

CVD - Chemical Vapor Deposition 

DD - Dynamical Decoupling 

DQD - Double Quantum Dot 

FID - Free Induction Decay 

GRAPE - GRadient Ascent Pulse Engineering 

LZ - Landau-Zener 

MRT - Macroscopic Resonant Tunneling 

NMR - Nuclear Magnetic Resonance 

PDD - Periodic Dynamical Decoupling 

QED - Quantum Electrodynamics 

QI - Qubit plus Impurity 

QD - Quantum Dot 

QPC - Quantum Point Contact 

RKKY - Ruderman-Kittel-Kasuya-Iosida 

RT -Random Telegraph 

RTN - Random Telegraph Noise 

SEM - Scanning Electron Microscopy 

SET - Single Electron Tunneling 

SO - Spin - Orbit 

SQUID - Superconducting Quantum Interference Device 

TLS - Two-Level System 

UDD - Uhrig Dynamical Decoupling 

UV - Ultra- Violet 



